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A multi-rate speech codec supports a 
plurality of encoding bit rate modes by adap- 
tively selecting encoding bit rate modes to 
match communication channel restrictions. In 
higher bit rate encoding modes, an accurate 
representation of speech through CELP (code 
excited linear prediction) and other associated 
modeling parameters are generated for higher 
quality decoding and reproduction. To sup- 
port lower bit rate encoding modes, a variety 
of techniques are applied many of which in- 
volve the classification of the input signal. The 
encoder utilizes gain normalization wherein 
LPC (linear predictive coding) gain provides a 
smoothing factor for combining both open and 
closed loop gains. The lower the LPC gain, 
the greater the open loop gain contribution to 
a gain normalization factor. The greater the 
LPC gain, the greater the closed loop gain con- 
tribution. For background noise, the smaller of 
the closed and open loop gains are used as the 
normalization factor. The normalization factor 
is limited by the LPC gain to prevent influenc- 
ing the coding quality. 
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Background 

1. Technical Field 

The present invention relates generally to speech encoding and decoding in voice 
communication systems; and, more particularly, it relates to various techniques used with code- 
excited linear prediction coding to obtain high quality speech reproduction through a limited bit 
rate communication channel. 
2. Related Art 

Signal modeling and parameter estimation play significant roles in communicating voice 
information with limited bandwidth constraints. To model basic speech sounds, speech signals 
are sampled as a discrete waveform to be digitally processed. In one type of signal coding 
technique called LPC (linear predictive coding), the signal value at any particular time index is 
modeled as a linear function of previous values. A subsequent signal is thus linearly predictable 
according to an earlier value. As a result, efficient signal representations can be determined by 
estimating and applying certain prediction parameters to represent the signal. 

Applying LPC techniques, a conventional source encoder operates on speech signals to 
extract modeling and parameter information for communication to a conventional source decoder 
via a communication channel. Once received, the decoder attempts to reconstruct a counterpart 
signal for playback that sounds to a human ear like the original speech. 

A certain amount of communication channel bandwidth is required to communicate the 
modeling and parameter information to the decoder. In embodiments, for example where the 
channel bandwidth is shared and real-time reconstruction is necessary, a reduction in the required 
bandwidth proves beneficial. However, using conventional modeling techniques, the quality 
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requirements in the reproduced speech limit the reduction of such bandwidth below certain 
levels. 

Further limitations and disadvantages of conventional systems will become apparent to 
one of skill in the art after reviewing the remainder of the present application with reference to 
the drawings. 
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SUMMARY OF THE INVENTION 

Various aspects of the present invention can be found in a speech encoding system using 
an analysis by synthesis coding approach on a speech signal. Therein, the speech encoding 
system has an encoder processing circuit and a plurality of codebooks that generate excitation 
vectors. The encoder processing circuit calculates open loop gain and closed loop gain. The 
encoder processing circuit selectively applies the open and closed loop gains in gain 
normalization processing. 

The selective application of the open loop gain and the closed loop gain by the encoder 
processing circuit may further involve the use of a weighting factor of linear predictive coding 
gain. The encoder processing circuit may use weighting factor to linearly combine the open and 
closed loop gains. 

In certain embodiments, the selective application of the open and closed loop gains by the 
encoder processing circuit comprises applying the lesser of the open loop gain and the closed 
loop gain to the background noise. When using the linear predictive coding gain as a weighting 
factor, the encoder processing circuit may exclude such background noise from the application of 
the weighting factor. 

Additionally, the encoder processing circuit may apply a maximum limit, a minimum 
limit, or both in gain normalization processing. 

Further aspects of the present invention can also be found in a method used by a speech 
encoding system that applies an analysis by synthesis coding approach to a speech signal. The 
method involves the identification of open and closed loop gains for combination of 
contributions therefrom to generate a gain normalization factor. 
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In some embodiments, such method may further involve the use of linear predictive 
coding gain to identify appropriate contributions of the open and closed loop gains. One specific 
way to accomplish this involves the use of the linear predictive coding gain as a weighting 
factor. 

When the speech signal comprises background noise, the encoder processing system may 
also select without combination either the open loop gain or the closed loop gain in generating 
the gain normalization factor. In some cases, this involves selection of the lesser of the open 
loop gain and the closed loop gain. 

Other aspects, advantages and novel features of the present invention will become 
apparent from the following detailed description of the invention when considered in conjunction 
with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. la is a schematic block diagram of a speech communication system illustrating the 
use of source encoding and decoding in accordance with the present invention. 

Fig. Ib is a schematic block diagram illustrating an exemplary communication device 
utilizing the source encoding and decoding functionality of Fig. la. 

Figs. 2-4 are functional block diagrams illustrating a multi-step encoding approach used 
by one embodiment of the speech encoder illustrated in Figs, la and lb. In particular, Fig. 2 is a 
functional block diagram illustrating of a first stage of operations performed by one embodiment 
of the speech encoder of Figs, la and lb. Fig. 3 is a functional block diagram of a second stage 
of operations, while Fig. 4 illustrates a third stage. 

Fig. 5 is a block diagram of one embodiment of the speech decoder shown in Figs, la and 
lb having corresponding functionality to that illustrated in Figs. 2-4. 

Fig. 6 is a block diagram of an alternate embodiment of a speech encoder that is built in 
accordance with the present invention. 

Fig. 7 is a block diagram of an embodiment of a speech decoder having corresponding 
functionality to that of the speech encoder of Fig. 6. 

Fig. 8 is a flow diagram illustrating the functionality of gain normalization such as that 
represented in the block 401 of Fig. 4 by an encoder built in accordance with the present 
invention. 

Fig. 9 is a flow diagram providing a more detailed description of one embodiment of gain 
normalization functionality of Fig. 8 
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Detailed Description 

Fig. la is a schematic block diagram of a speech communication system illustrating the 
use of source encoding and decoding in accordance with the present invention. Therein, a 
speech communication system 100 supports communication and reproduction of speech across a 
communication channel 103. Although it may comprise for example a wire, fiber or optical link, 
the communication channel 103 typically comprises, at least in part, a radio frequency link that 
often must support multiple, simultaneous speech exchanges requiring shared bandwidth 
resources such as may be found with cellular telephony embodiments. 

Although not shown, a storage device may be coupled to the communication channel 103 
to temporarily store speech information for delayed reproduction or playback, e.g., to perform 
answering machine functionality, voiced email, etc. Likewise, the communication channel 103 
might be replaced by such a storage device in a single device embodiment of the communication 
system 100 that, for example, merely records and stores speech for subsequent playback. 

In particular, a microphone 1 1 1 produces a speech signal in real time. The microphone 
1 1 1 delivers the speech signal to an A/D (analog to digital) convener 1 15. The A/D convener 
1 15 converts the speech signal to a digital form then delivers the digitized speech signal to a 
speech encoder 1 17. 

The speech encoder 1 17 encodes the digitized speech by using a selected one of a 
plurality of encoding modes. Each of the plurality of encoding modes utilizes particular 
techniques that attempt to optimize quality of resultant reproduced speech. While operating in 
any of the plurality of modes, the speech encoder 1 17 produces a series of modeling and 
parameter information (hereinafter "speech indices"), and delivers the speech indices to a 
channel encoder 1 19. 
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The channel encoder 1 19 coordinates with a channel decoder 13 1 to deliver the speech 
indices across the communication channel 103. The channel decoder 131 forwards the speech 
indices to a speech decoder 133. While operating in a mode that corresponds to that of the 
speech encoder 1 17, the speech decoder 133 attempts to recreate the original speech from the 
speech indices as accurately as possible at a speaker 137 via a D/A (digital to analog) converter 



135. 



The speech encoder 1 17 adaptively selects one of the plurality of operating modes based 
on the data rate restrictions through the communication channel 103. The communication 
channel 103 comprises a bandwidth allocation between the channel encoder 1 19 and the channel 
decoder 131. The allocation is established, for example, by telephone switching networks 
wherein many such channels are allocated and reallocated as need arises. In one such 
embodiment, either a 22.8 kbps (kilobits per second) channel bandwidth, i.e., a full rate channel, 
or a 1 1.4 kbps channel bandwidth, i.e., a half rate channel, may be allocated. 

With the full rate channel bandwidth allocation, the speech encoder 1 17 may adaptively 
select an encoding mode that supports a bit rate of 1 1 .0. 8.0. 6.65 or 5.8 kbps. The speech 
encoder 1 17 adaptively selects an either 8.0, 6.65, 5.8 or 4.5 kbps encoding bit rate mode when 
only the half rate channel has been allocated. Of course these encoding bit rates and the 
aforementioned channel allocations are only representative of the present embodiment. Other 
variations to meet the goals of alternate embodiments are contemplated. 

With either the full or half rate allocation, the speech encoder 1 17 attempts to 
communicate using the highest encoding bit rate mode that the allocated channel will support. If 
the allocated channel is or becomes noisy or otherwise restrictive to the highest or higher 
encoding bit rates, the speech encoder 1 17 adapts by selecting a lower bit rate encoding mode. 
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Similarly, when the communication channel 103 becomes more favorable, the speech encoder 
1 17 adapts by switching to a higher bit rate encoding mode. 

With lower bit rate encoding, the speech encoder 1 17 incorporates various techniques to 
generate better low bit rate speech reproduction. Many of the techniques applied are based on 
characteristics of the speech itself. For example, with lower bit rate encoding, the speech 
encoder 117 classifies noise, unvoiced speech, and voiced speech so that an appropriate 
modeling scheme corresponding to a particular classification can be selected and implemented. 
Thus, the speech encoder 1 17 adaptively selects from among a plurality of modeling schemes 
those most suited for the current speech. The speech encoder 1 17 also applies various other 
techniques to optimize the modeling as set forth in more detail below. 

Fig. lb is a schematic block diagram illustrating several variations of an exemplary 
communication device employing the functionality of Fig. la. A communication device 151 
comprises both a speech encoder and decoder for simultaneous capture and reproduction of 
speech. Typically within a single housing, the communication device 151 might, for example, 
comprise a cellular telephone, portable telephone, computing system, etc. Alternatively, with 
some modification to include for example a memory element to store encoded speech 
information the communication device 151 might comprise an answering machine, a recorder, 
voice mail system, etc. 

A microphone 155 and an A/D converter 157 coordinate to deliver a digital voice signal 
to an encoding system 159. The encoding system 159 performs speech and channel encoding 
and delivers resultant speech information to the channel. The delivered speech information may 
be destined for another communication device ( not shown) at a remote location. 
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As speech information is received, a decoding system 165 performs channel and speech 
decoding then coordinates with a D/A converter 167 and a speaker 169 to reproduce something 
that sounds like the originally captured speech. 

The encoding system 159 comprises both a speech processing circuit 185 that performs 
speech encoding, and a channel processing circuit 187 that performs channel encoding. 
Similarly, the decoding system 165 comprises a speech processing circuit 189 that performs 
speech decoding, and a channel processing circuit 191 that performs channel decoding. 

Although the speech processing circuit 185 and the channel processing circuit 187 are 
separately illustrated, they might be combined in part or in total into a single unit. For example, 
the speech processing circuit 185 and the channel processing circuitry 187 might share a single 
DSP (digital signal processor) and/or other processing circuitry. Similarly, the speech processing 
circuit 1 89 and the channel processing circuit 191 might be entirely separate or combined in pan 
or in whole. Moreover, combinations in whole or in pan might be applied to the speech 
processing circuits 185 and 189, the channel processing circuits 187 and 191, the processing 
circuits 185, 187, 189 and 191, or otherwise. 

The encoding system 159 and the decoding system 165 both utilize a memory 161. The 
speech processing circuit 185 utilizes a fixed codebook 181 and an adaptive codebook 183 of a 
speech memory 177 in the source encoding process. The channel processing circuit 187 utilizes 
a channel memory 175 to perform channel encoding. Similarly, the speech processing circuit 
189 utilizes the fixed codebook 181 and the adaptive codebook 183 in the source decoding 
process. The channel processing circuit 187 utilizes the channel memory 175 to perform channel 
decoding. 
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Although the speech memory 177 is shared as illustrated, separate copies thereof can be 
assigned for the processing circuits 185 and 189. Likewise, separate channel memory can be 
allocated to both the processing circuits 187 and 191. The memory 161 also contains software 
utilized by the processing circuits 185,187,189 and 191 to perform various functionality required 
in the source and channel encoding and decoding processes. 

Figs. 2-4 are functional block diagrams illustrating a multi-step encoding approach used 
by one embodiment of the speech encoder illustrated in Figs, la and lb. In particular, Fig. 2 is a 
functional block diagram illustrating of a first stage of operations performed by one embodiment 
of the speech encoder shown in Figs, la and lb. The speech encoder, which comprises encoder 
processing circuitry, typically operates pursuant to software instruction carrying out the 
following functionality. 

At a block 215, source encoder processing circuitry performs high pass filtering of a 
speech signal 211. The filter uses a cutoff frequency of around 80 Hz to remove, for example, 
60 Hz power line noise and other lower frequency signals. After such filtering, the source 
encoder processing circuitry applies a perceptual weighting filter as represented by a block 219. 
The perceptual weighting filter operates to emphasize the valley areas of the filtered speech 
signal. 

If the encoder processing circuitry selects operation in a pitch preprocessing (PP) mode 
as indicated at a control block 245, a pitch preprocessing operation is performed on the weighted 
speech signal at a block 225. The pitch preprocessing operation involves warping the weighted 
speech signal to match interpolated pitch values that will be generated by the decoder processing 
circuitry. When pitch preprocessing is applied, the warped speech signal is designated a first 
target signal 229. If pitch preprocessing is not selected the control block 245, the weighted 
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speech signal passes through the block 225 without pitch preprocessing and is designated the 
first target signal 229. 

As represented by a block 255, the encoder processing circuitry applies a process wherein 
a contribution from an adaptive codebook 257 is selected along with a corresponding gain 257 
which minimize a first error signal 253. The first error signal 253 comprises the difference 
between the first target signal 229 and a weighted, synthesized contribution from the adaptive 
codebook 257. 

At blocks 247, 249 and 25 1, the resultant excitation vector is applied after adaptive gain 
reduction to both a synthesis and a weighting filter to generate a modeled signal that best 
matches the first target signal 229. The encoder processing circuitry uses LPC (linear predictive 
coding) analysis, as indicated by a block 239, to generate filter parameters for the synthesis and 
weighting filters. The weighting filters 219 and 25 1 are equivalent in functionality. 

Next, the encoder processing circuitry designates the first error signal 253 as a second 
target signal for matching using contributions from a fixed codebook 261. The encoder 
processing circuitry searches through at least one of the plurality of subcodebooks within the 
fixed codebook 261 in an attempt to select a most appropriate contribution while generally 
attempting to match the second target signal. 

More specifically, the encoder processing circuitry selects an excitation vector, its 
corresponding subcodebook and gain based on a variety of factors. For example, the encoding 
bit rate, the degree of minimization, and characteristics of the speech itself as represented by a 
block 279 are considered by the encoder processing circuitry at control block 275. Although 
many other factors may be considered, exemplary characteristics include speech classification, 
noise level, sharpness, periodicity, etc. Thus, by considering other such factors, a first 
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subcodebook with its best excitation vector may be selected rather than a second subcodebook s 
best excitation vector even though the second subcodebook's better minimizes the second target 
signal 265. 

Fig. 3 is a functional block diagram depicting of a second stage of operations performed by 
the embodiment of the speech encoder illustrated in Fig. 2. In the second stage, the speech 
encoding circuitry simultaneously uses both the adaptive the fixed codebook vectors found in the 
first stage of operations to minimize a third error signal 311. 

The speech encoding circuitry searches for optimum gain values for the previously 
identified excitation vectors ( in the first stage) from both the adaptive and fixed codebooks 257 
and 261. As indicated by blocks 307 and 309, the speech encoding circuitry identifies the 
optimum gain by generating a synthesized and weighted signal, i.e., via a block 301 and 303. that 
best matches the first target signal 229 (which minimizes the third error signal 3 1 1 ). Of course if 
processing capabilities permit, the first and second stages could be combined wherein joint 
optimization of both gain and adaptive and fixed codebook rector selection could be used. 

Fig. 4 is a functional block diagram depicting of a third stage of operations performed by 
the embodiment of the speech encoder illustrated in Figs. 2 and 3. The encoder processing 
circuitry applies gain normalization, smoothing and quantization, as represented by blocks 401 , 
403 and 405, respectively, to the jointly optimized gains identified in the second stage of encoder 
processing. Again, the adaptive and fixed codebook vectors used are those identified in the first 
stage processing. 

With normalization, smoothing and quantization functionally applied, the encoder 
processing circuitry has completed the modeling process. Therefore, the modeling parameters 
identified are communicated to the decoder. In particular, the encoder processing circuitry 
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delivers an index to the selected adaptive codebook vector to the channel encoder via a 
multiplexor 419. Similarly, the encoder processing circuitry delivers the index to the selected 
fixed codebook vector, resultant gains, synthesis filter parameters, etc., to the muliplexor 419. 
The multiplexor 419 generates a bit stream 421 of such information for delivery to the channel 
encoder for communication to the channel and speech decoder of receiving device. 

Fig. 5 is a block diagram of an embodiment illustrating functionality of speech decoder 
having corresponding functionality to that illustrated in Figs. 2-4. As with the speech encoder, 
the speech decoder, which comprises decoder processing circuitry, typically operates pursuant to 
software instruction carrying out the following functionality. 

A demultiplexer 51 1 receives a bit stream 513 of speech modeling indices from an often 
remote encoder via a channel decoder. As previously discussed, the encoder selected each index 
value during the multi-stage encoding process described above in reference to Figs. 2-4. The 
decoder processing circuitry utilizes indices, for example, to select excitation vectors from an 
adaptive codebook 515 and a fixed codebook 519, set the adaptive and fixed codebook gains at a 
block 52 1 , and set the parameters for a synthesis filter 53 1 . 

With such parameters and vectors selected or set, the decoder processing circuitry 
generates a reproduced speech signal 539. In particular, the codebooks 515 and 519 generate 
excitation vectors identified by the indices from the demultiplexer 51 1 . The decoder processing 
circuitry applies the indexed gains at the block 521 to the vectors which are summed. At a block 
527, the decoder processing circuitry modifies the gains to emphasize the contribution of vector 
from the adaptive codebook 515. At a block 529, adaptive tilt compensation is applied to the 
combined vectors with a goal of flattening the excitation spectrum. The decoder processing 
circuitry performs synthesis filtering at the block 531 using the flattened excitation signal. 
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Finally, to generate the reproduced speech signal 539, post filtering is applied at a block 535 
deemphasizing the valley areas of the reproduced speech signal 539 to reduce the effect of 
distortion. 

In the exemplary cellular telephony embodiment of the present invention, the A/D 
converter 1 15 (Fig. la) will generally involve analog to uniform digital PCM including: 1) an 
input level adjustment device; 2) an input anti-aliasing filter; 3) a sample-hold device sampling at 
8 kHz; and 4) analog to uniform digital conversion to 13-bit representation. 

Similarly, the D/A converter 135 will generally involve uniform digital PCM to analog 
including: 1) conversion from 13-bit/8 kHz uniform PCM to analog; 2) a hold device; 3) 
reconstruction filter including x/sin(x) correction; and 4) an output level adjustment device. 

In terminal equipment, the A/D function may be achieved by direct conversion to 13-bit 
uniform PCM format, or by conversion to 8-bit/A-law compounded format. For the D/A 
operation, the inverse operations take place. 

The encoder 1 17 receives data samples with a resolution of 13 bits left justified in a 
16-bit word. The three least significant bits are set to zero. The decoder 133 outputs data in the 
same format. Outside the speech codec, further processing can be applied to accommodate 
traffic data having a different representation. 

A specific embodiment of an AMR (adaptive multi-rate) codec with the operational 
functionality illustrated in Figs. 2-5 uses five source codecs with bit-rates 1 1.0, 8.0, 6.65, 5.8 and 
4.55 kbps. Four of the highest source coding bit-rates are used in the full rate channel and the 
four lowest bit-rates in the half rate channel. 

All five source codecs within the AMR codec are generally based on a code-excited 
linear predictive (CELP) coding model. A 10th order linear prediction (LP), or short-term, 
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synthesis filter, e.g., used at the blocks 249, 267, 301, 407 and 531 (of Figs. 2-5), is used which 
is given by: 

H(z) = J— = — j . (i) 

where 2 jt i = 1 m, are the (quantized) linear prediction (LP) parameters. 

A long-term filter, i.e., the pitch synthesis filter, is implemented using the either an 
adaptive codebook approach or a pitch pre-processing approach. The pitch synthesis filter is 
given by: 
1 1 

Wz) \-g p z- Tt (2) 
where T is the pitch delay and g p is the pitch gain. 

With reference to Fig. 2, the excitation signal at the input of the short-term LP synthesis 
filter at the block 249 is constructed by adding two excitation vectors from the adaptive and the 
fixed codebooks 257 and 261, respectively. The speech is synthesized by feeding the two 
properly chosen vectors from these codebooks through the short-term synthesis filter at the block 
249 and 267, respectively. 

The optimum excitation sequence in a codebook is chosen using an analysis-by-synthesis 
search procedure in which the error between the original and synthesized speech is minimized 
according to a perceptually weighted distortion measure. The perceptual weighting filter, e.g., at 
the blocks 251 and 268, used in the analysis-by-synthesis search technique is given by: 

A(z/y 2 ) (3) 
where A(z) is the unquantized LP filter and 0 < y 2 < Yi £ 1 are the perceptual weighting 
factors. The values y, = [0.9, 0.94] and y 2 = 0.6 are used. The weighting filter, e.g., at the 
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blocks 251 and 268. uses the unquantized LP parameters while the formant synthesis filter, e.g.. 
at the blocks 249 and 267, uses the quantized LP parameters. Both the unquantized and 
quantized LP parameters are generated at the block 239. 

The present encoder embodiment operates on 20 ms (millisecond) speech frames 
corresponding to 160 samples at the sampling frequency of 8000 samples per second. At each 
160 speech samples, the speech signal is analyzed to extract the parameters of the CELP model, 
i.e., the LP filter coefficients, adaptive and fixed codebook indices and gains. These parameters 
are encoded and transmitted. At the decoder, these parameters are decoded and speech is 
synthesized by filtering the reconstructed excitation signal through the LP synthesis filter. 

More specifically, LP analysis at the block 239 is performed twice per frame but only a 
single set of LP parameters is converted to line spectrum frequencies (LSF) and vector quantized 
using predictive multi-stage quantization (PMVQ). The speech frame is divided into subframes. 
Parameters from the adaptive and fixed codebooks 257 and 261 are transmitted every subframe. 
The quantized and unquantized LP parameters or their interpolated versions are used depending 
on the subframe. An open-loop pitch lag is estimated at the block 241 once or twice per frame 
for PP mode or LTP mode, respectively. 

Each subframe, at least the following operations are repeated. First, the encoder 
processing circuitry (operating pursuant to software instruction) computes x( n ) , the first target 
signal 229, by filtering the LP residual through the weighted synthesis filter W( z)H( z) with the 
initial states of the filters having been updated by filtering the error between LP residual and 
excitation. This is equivalent to an alternate approach of subtracting the zero input response of 
the weighted synthesis filter from the weighted speech signal. 
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Second, the encoder processing circuitry computes the impulse response, Hn K of the 

weighted synthesis filter. Third, in the LTP mode, closed-loop pitch analysis is performed to 

find the pitch lag and gain, using the first target signal 229, x( n ) , and impulse response, H n ) , 

by searching around the open-loop pitch lag. Fractional pitch with various sample resolutions 
are used. 

In the PP mode, the input original signal has been pitch-preprocessed to match the 
interpolated pitch contour, so no closed-loop search is needed. The LTP excitation vector is 
computed using the interpolated pitch contour and the past synthesized excitation. 

Fourth, the encoder processing circuitry generates a new target signal jc 2 ( n ) , the second 
target signal 253, by removing the adaptive codebook contribution (filtered adaptive code vector) 
from x(n) The encoder processing circuitry uses the second target signal 253 in the fixed 
codebook search to find the optimum innovation. 

Fifth, for the 1 1 .0 kbps bit rate mode, the gains of the adaptive and fixed codebook are 
scalar quantized with 4 and 5 bits respectively (with moving average prediction applied to the 
fixed codebook gain). For the other modes the gains of the adaptive and fixed codebook are 
vector quantized (with moving average prediction applied to the fixed codebook gain). 

Finally, the filter memories are updated using the deteimined excitation signal for finding 
the first target signal in the next subframe. 

The bit allocation of the AMR codec modes is shown in table 1 . For example, for each 
20 ms speech frame, 220, 160, 133 , 1 16 or 91 bits are produced, corresponding to bit rates of 
1 1.0, 8.0, 6.65, 5.8 or 4.55 kbps, respectively. 
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Table 1: Bit allocation of the AMR coding algorithm for 20 ms frame 



| CODING RATE 


1! OKBPS 




i 6.65KBPS 


1 5 80KB PS 


1 4 55KBPS 




20ms 






LPC order 


10" -order 


Quantization 




1 predictor. 
Obit/frame 


2 predictors 
1 bit/frame 




28 bit/frame 


24 bit/frame 


18 


LPC interpolation 


2 bits/frame 


2 bits/f | 0 


2 bits/f 


o 


0 


0 


Coding mode bit 


Obit 


Obit 


1 bit/frame 


Obit 


Obit 




LTP 


LTP | 


LTP 


PP 


PP 


PP 






5ms 








30 bits/frame (9696) 


8585 


8585 


0008 


0008 


0008 






20 


13 


18 








9 bits (scalar) 










■ I i i 1 




Total 


220 bits/frame 


160 


133 


133 


116 


91 



With reference to Fig. 5, the decoder processing circuitry, pursuant to software control. 



reconstructs the speech signal using the transmitted modeling indices extracted from the received 
bit stream by the demultiplexor 511. The decoder processing circuitry decodes the indices to 
obtain the coder parameters at each transmission frame. These parameters are the LSF vectors, 
the fractional pitch lags, the innovative code vectors, and the two gains. 

The LSF vectors are converted to the LP filter coefficients and interpolated to obtain LP 
filters at each subframe. At each subframe, the decoder processing circuitry constructs the 
excitation signal by: 1) identifying the adaptive and innovative code vectors from the codebooks 
5 15 and 5 19; 2) scaling the contributions by their respective gains at the block 52 1 ; 3) summing 
the scaled contributions; and 3) modifying and applying adaptive tilt compensation at the blocks 
527 and 529. The speech signal is also reconstructed on a subframe basis by filtering the 
excitation through the LP synthesis at the block 531. Finally, the speech signal is passed through 
an adaptive post filter at the block 535 to generate the reproduced speech signal 539. 

The AMR encoder will produce the speech modeling information in a unique sequence 
and format, and the AMR decoder receives the same information in the same way. The different 
parameters of the encoded speech and their individual bits have unequal importance with respect 
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to subjective quality. Before being submitted to the channel encoding function the bits are 
rearranged in the sequence of importance. 

Two pre-processing functions are applied prior to the encoding process: high-pass 
filtering and signal down-scaling. Down-scaling consists of dividing the input by a factor of 2 
to reduce the possibility of overflows in the fixed point implementation. The high-pass filtering 
at the block 215 (Fig. 2) serves as a precaution against undesired low frequency components. A 
filter with cut off frequency of 80 Hz is used, and it is given by: 

H (z) _ 0.92727435-1. 854494 lz' 1 +0.927274352" 2 
W Z 1 - 1.90594652 +0.91 14024z" 2 

Down scaling and high-pass filtering are combined by dividing the coefficients of the numerator 

of H M (z)by 2. 

Short-term prediction, or linear prediction (LP) analysis is performed twice per speech 
frame using the autocorrelation approach with 30 ms windows. Specifically, two LP analyses 
are performed twice per frame using two different windows. In the first LP analysis 
(LP_analysis_l ), a hybrid window is used which has its weight concentrated at the fourth 
subframe. The hybrid window consists of two parts. The first part is half a Hamming window, 
and the second part is a quarter of a cosine cycle. The window is given by: 





n= 215 to 239 



In the second LP analysis (LP_analysis_2), a symmetric Hamming window is used. 
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w : (n) = < 



0.54 + 0.46cos 



0.54 — 0.46 CO! 




n = 0toll9,L = 120 



n = 120to 239 



past frame current frame 
55 160 




future frame 



25 (samples) 



In either LP analysis, the autocorrelations of the windowed speech s (n\n = 0,239 are computed 
by: 



A 60 Hz bandwidth expansion is used by lag windowing, the autocorrelations using the window: 



Moreover, r(0) is multiplied by a white noise correction factor 1 .0001 which is equivalent to 
adding a noise floor at -40 dB. 



used to obtain the reflection coefficients k i and LP filter coefficients a, , i = 1 , 10 using the 
Levinson-Durbin algorithm. Furthermore, the LP filter coefficients a, are used to obtain the 
Line Spectral Frequencies (LSFs). 

The inteipolated unquantized LP parameters are obtained by interpolating the LSF 
coefficients obtained from the LP analysis. 1 and those from LP_analysis_2 as: 

q, (n) = 0.5q A (n - 1) + Q.5q 2 (n) 
q 3 (n) = 0.5^ 2 (n) + 0.5^ 4 (n) 



r(k) = ^s (n)s (n-kX k = 0,10. 




The modified autocorrelations r (0) = 1.0001r(0) and r (it) = r(*)w^ (it), k = 1,10 are 
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where q t (n) is the interpolated LSF for subframe 1 . is the LSF of subframe 2 obtained 

from LP_analysis_2 of current frame. q 3 (n) is the interpolated LSF for subframe 3. q A (n - 1) is 

the LSF (cosine domain) from LP_analysis_l of previous frame, and q 4 (n) is the LSF for 

subframe 4 obtained from LP_analysis_l of current frame. The interpolation is carried out in the 
cosine domain. 

A VAD (Voice Activity Detection) algorithm is used to classify input speech frames into 
either active voice or inactive voice frame (background noise or silence) at a block 235 (Fig. 2). 

The input speech s{n) is used to obtain a weighted speech signal s w (n) by passing s(n) 
through a filter: 



{/*)' 



That is, in a subframe of size L_SF, the weighted speech is given by: 

10 10 

s„ (n) = s(n) + £ a k y\ s(n -i)-^ a,y' j„ (n - /), n = 0. L_SF - 1 . 



A voiced/unvoiced classification and mode decision within the block 279 using the input 
speech s(n) and the residual r w (n) is derived where: 

10 

r w in) = s(n) + ^atfsfr - 1), n = 0,L_SF- 1 . 

i=i 

The classification is based on four measures: 1) speech sharpness P1_SHP^ % f normalized one 
delay correlation P2.R1; 3) normalized zero-crossing rate P3.ZC; and 4) normalizedLP residual 
energy P4_RE. 

The speech sharpness is given by: 
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L 



^abs(rjn)) 



P\_SHP = 



MaxL 



where Max is the maximum of abs(r w (n)) over the specified interval of length L . The 



where sgn is the sign function whose output is either 1 or -1 depending that the input sample is 
positive or negative. Finally, the normalized LP residual energy is given by: 

P4_RE = \-yJlpc_gain 



where Ipc _ gain = fj (1 - kf ) , where k, are the reflection coefficients obtained from LP 
analysis_ 1 . 

The voiced/unvoiced decision is derived if the following conditions are met: 

if P2_R\< 0.6 and PI _ SHP > 0.2 set mode = 2, 
if P3 _ ZC> 0.4 and PI _SHP> 0.1 8 set mode =2, 
if PA _ RE < 0.4 and PI _ SHP > 0.2 set mode « 2, 
if (P2 _ Rl < - 1 .2 + 3.2P1 _ SHP) set VUV = -3 
if{P4_RE< -0.2 1 + 1 .4286P1 _ SHP) set VUV = -3 
if (P3 _ ZC> 0.8 - 0.6P1 _ SHP) set VUV = -3 
*7(P4 _ RE < 0. 1) set VUV = -3 

Open loop pitch analysis is performed once or twice (each 10 ms) per frame depending 
on the coding rate in order to find estimates of the pitch lag at the block 241 (Fig. 2). It is based 



normalized one delay correlation and normalized zero-crossing rate are given by: 



L-l 



Jj(n)i(n + 1) 

D1 P.- «=0 




P3_ ZC = -L [i sgn[5(/)] - sgn[j(/ - 1)] l], 



10 
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on the weighted speech signal j„ (n + n m ).n = 0.1 ,79, in which n m defines the location of this 

signal on the first half frame or the last half frame. In the first step, four maxima of the 
correlation: 

79 

c k = X 5 w("m + n)s w (n m + n~k) 

n=0 

are found in the four ranges 17....33. 34....67, 68.. ..135, 136.. ..145, respectively. The retained 
maxima C^. / = 1,23,4, are nonnalized by dividing by: 

■J2,„ s U n m + «-*). «' = 1 4, respectively. 

The nonnalized maxima and corresponding delays are denoted by (/?,,*,),/= 1,2,3,4. 

In the second step, a delay, among the four candidates, is selected by maximizing the 
four nonnalized correlations. In the third step. *, is probably conected to *, (i<7) by favoring the 
lower ranges. That is, *, (i<f) is selected if*, is within fk/m-4, k/m+4].m=2.3,4.5. and if 

> k, 0.95 7 "' D, i < I, where D is 1.0, 0.85, or 0.65, depending on whether the previous frame 
is unvoiced, the previous frame is voiced and *, is in the neighborhood (specified by ± 8) of the 
previous pitch lag, or the previous two frames are voiced and *, is in the neighborhood of the 
previous two pitch lags. The final selected pitch lag is denoted by T^. 

A decision is made every frame to either operate the LTP (long-term prediction) as the 
traditional CELP approach (LTP_mode=l), or as a modified time warping approach 
(LTP_mode=0) herein referred to as PP (pitch preprocessing). For 4.55 and 5.8 kbps encoding 
bit rates. LTP.mode is set to 0 at all times. For 8.0 and 1 1 .0 kbps. LTP.mode is set to 1 all of 
the time. Whereas, for a 6.65 kbps encoding bit rate, the encoder decides whether to operate in 
the LTP or PP mode. During the PP mode, only one pitch lag is transmitted per coding frame. 
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For 6.65 kbps, the decision algorithm is as follows. First, at the block 24 1 . a prediction 
of the pitch lag pit for the current frame is determined as follows: 

if (LTPJAODE_m — 1 ) 

pit = lagl\ + 2A*(lagJ[3]-lagl\ ); 



pit = lagjl\] + 2.75*(lagJ[3]-lagJ[\J); 
where L7P_mod*_m is previous frame ZJP.mode, lag _f[l]Jag _f[3) are the past closed 
loop pitch lags for second and fourth subframes respectively, lagl is the current frame open-loop 
pitch lag at the second half of the frame, and t lagll is the previous frame open-loop pitch lag at 
the first half of the frame. 

Second, a normalized spectrum difference between the Line Spectrum Frequencies (LSF) 
of current and previous frame is computed as: 

e _ Isf = £ abs(LSF(i) -LSF_ m(i)) , 



if (abs(piulagl) < TH and abs(lagj[ 3 J-lagl) < lagl*0.2 ) 

if(Rp > 05 && pgain_past > 0.7 and ejsf < 05/30 ) LTP_ mod e = 0; 
else LTP_ mode = 1; 

where Rp is current frame normalized pitch correlation, pgain _ past is the quantized pitch gain 

from the fourth subframe of the past frame, TH = MIN(lagl*0A, 5 and TH = MAX( 2.0, TH) . 

The estimation of the precise pitch lag at the end of the frame is based on the normalized 
correlation: 



else 



L 



£ s w (n + nl) s w (n + n\ - k) 
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where j„ (n + nl), n = 0, 1 L-l, represents the last segment of the weighted speech signal 

including the look-ahead ( the look-ahead length is 25 samples), and the size L is defined 
according to the open-loop pitch lag 7^ with the corresponding normalized correlation C T : 
if(C T >0.6) 

op 

L = max{ 50. Top/ 

L « min{ 80, L } 

else 
L = 80 

In the first step, one integer lag k is selected maximizing the R k in the range 

* *l T op ~ 10, T op + 10] bounded by [17, 145]. Then, the precise pitch lag P m and the 

corresponding index I m for the current frame is searched around the integer lag, [k-1, k+ 1 ]. by 
up-sampling R k . 

The possible candidates of the precise pitch lag are obtained from the table named as 

PitLagTab8b[i], i=0J 127. In the last step, the precise pitch lag P„ = PitLagTab8b[l m ] is 

possibly modified by checking the accumulated delay due to the modification of the speech 
signal: 

'/ (Tncc >5) I m «= min{ I m + 1, 127} , and 
if&acc < -5) / m «=max{/„ -1,0}. 

The precise pitch lag could be modified again: 

f face > 10) / m <=min{/ m + l, 127}, and 
if&acc < -10) / m <=max{/ m -1,0}. 

The obtained index I„ will be sent to the decoder. 

The pitch lag contour, T c (n) , is defined using both the current lag P m and the previous 
lag P„.j: 
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'/ ( K - < 02 min{/> m , } ) 

W = + n(P m -P M .i)/L / , *=0.1 L r -1 

T f (/i) =P m ,n=Lf,...,nO 

else 

t c (") =Pm-i .n=0,l 39; 

t c («) = , n=40 170 

where L/= 160 is the frame size. 

One frame is divided into 3 subframes for the long-term preprocessing. For the first two 
subframes, the subframe size, L,, is 53, and the subframe size for searching, is 70. For the 
last subframe, L s is 54 and L„ is: 

L„ = min{70, L s + L w - 10- T^} , 

where 1*^=25 is the look-ahead and the maximum of the accumulated delay is limited to 
14. 

The target for the modification process of the weighted speech temporally memorized in 

{ S w (mO+ n), n = 0,1 L„ - 1 } is calculated by warping the past modified weighted speech 

buffer, s w ( mO + n ), n< 0, with the pitch lag contour, r c (n + m • L s ), m = 0,1,2 , 

s w (mO + n) = £ s w (mO + n - T c («) + 1) /,(/, T lc (n)), n = 0.l L„ - 1, 

<=-// 

where 7"c0i; and T,dn) are calculated by: 

T c (n) =tnmc{x c (n + m- L s )} , 
T lc (n) = z c (n)-T c (n), 

m is subframe number. I s {i,T lc {n)) is a set of interpolation coefficients, and /, is 10. Then, the 

target for matching, s, (n), n = 0,1 L sr - 1, is calculated by weighting 

S w (mO + n), n= 0,1 L sr - 1, in the time domain: 

l(n) = n s w (mO+n)/L t , n = 0,l L,-l, 

I (n)=s w (m0 + n),« = L s L„ - 1 
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The local integer shifting range fSRO. SRI J for searching for the best local delay is 
computed as the following: 

if speech is unvoiced 
SR0=-l t 
SR1=1, 

else 

SR0=round{ -4 minfLO, maxfO.O. 1-0.4 (P S h-0.2)}J), 
SRl=round( 4 minfLO, maxfO.O, 1-0.4 (P^r02)}}} t 

where P S h=max{P shh P sh2 }, P shl is the average to peak ratio (i.e., sharpness) from the target 



where nO = /r«nc{m0+ + 05} (here, m is subframe number and is the previous 
accumulated delay). 

In order to find the best local delay, x opu at the end of the current processing subframe, a 
normalized correlation vector between the original weighted speech signal and the modified 
matching target is defined as: 



signal: 




L sr max{\s w (mO+n% n = 0,1,..., L, r - 1} 



and P sh2 is the sharpness from the weighted speech signal: 



1 £\s w (n + n0+L s /2)\ 



(L sr - L s 1 2)max{|s w (n + n0 + L s 1 2% n = 0,1 L sr - L s 1 2 - 1} 
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A best local delay in the integer domain, k optf is selected by maximizing Rrfk) in the range of 



If RKkopi)<0.5, k r is set to zero. 

In order to get a more precise local delay in the range { k r -0.75+0JjJ=0J t ...15J around 



where {I/iJ)) is a set of interpolation coefficients. The optimal fractional delay index, j t 
selected by maximizing R/j). Finally, the best local delay, r opf , at the end of the current 
processing subframe, is given by, 

V= *r -0.75 + 0.1 
The local delay is then adjusted by: 



The modified weighted speech of the current subframe, memorized in 

{ s w (m0 + n), n = 0,1 L s - 1 } to update the buffer and produce the second target signal 253 for 

searching the fixed codebook 261, is generated by warping the original weighted speech 
{ s w (n) } from the original time region, 

I" 10 + *«■ . w<> + ^ c + L M + ] , 
to the modified time region, 
[mO, mO+LJ: 



k 6 [5/fO,S/?l] , which is corresponding to the real delay: 



k r , Rrfk) is interpolated to obtain the fractional correlation vector, R/j)> by: 



8 



R/U) = £ RAkop, + 'j + 0 //(i. y), y = 0,1 15 , 
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i H .(mO + n)= Xm«0 + » + W + 0 /,0\r w (n)). n ^ 0,1 L - 1, 

where Tw(n) and 7Wn; are calculated by: 

T w (n) = mmc{T^ c +n r op , /I,}, 
Tfwin^r^ + n r opf I L s -T w (n) y 

{I s {i*T m (n))) is a set of interpolation coefficients. 

After having completed the modification of the weighted speech for the current subframe, 
the modified target weighted speech buffer is updated as follows: 
s w (n) ^s w (n + L s ), n = 0,l n m -1. 

The accumulated delay at the end of the current subframe is renewed by: 

acc *acc T h opt * 

Prior to quantization the LSFs are smoothed in order to improve the perceptual quality. 
In principle, no smoothing is applied during speech and segments with rapid variations in the 
spectral envelope. During non-speech with slow variations in the spectral envelope, smoothing 
is applied to reduce unwanted spectral variations. Unwanted spectral variations could typically 
occur due to the estimation of the LPC parameters and LSF quantization. As an example, in 
stationary noise-like signals with constant spectral envelope introducing even very small 
variations in the spectral envelope is picked up easily by the human ear and perceived as an 
annoying modulation. 

The smoothing of the LSFs is done as a running mean according to: 
Wi Xn) = p(n) • Isf t ( n - I) + (1 - 0(n)) Isf.est^nl i = 1 10 
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where fe/_«r,(n)isthe /'* estimated LSF of frame n,and lsf t (n) is the i" 1 LSF for quantization 
of frame n . The parameter /J(n) controls the amount of smoothing, e.g. if fi(n) is zero no 
smoothing is applied. 

fi(n) is calculated from the VAD information (generated at the block 235) and two 
estimates of the evolution of the spectral envelope. The two estimates of the evolution are 
defined as: 

10 

ASP = (n) -4sf_est { (n - 1)) 2 

&SP m = £(/*/_«r, (n) - majsft (n - 1)) 2 

maJsf,(n) = - ma_isf,(n - 1) + (1 - 0(n)) • Isf_est,(n). / = 1 10 
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The parameter fi(n) is controlled by the following logic: 

Step 1 : 

if(Vad = 1 1 PastVad = 1 1 k, > 0.5) 
P(n) = 0.0 

^«/(^n»*.<™ (n - 1) > 0 & (ASP > 0.0015 1 ASP m > 0.0024)) 

0(n) = 0.0 
elseifiN^ ^in - 1) > 1 & ASP > 0.0025) 

Step 2 : 

//(VW = 0 & PastVad = 0) 

<») = ^mode.tm. (« ~ D + 1 

«W«*>.<n)>5) 

(«) = 5 

else 
endif 

where *, is the first reflection coefficient. 

In step 1, the encoder processing circuitry checks the VAD and the evolution of the 
spectral envelope, and performs a full or partial reset of the smoothing if required. In step 2, the 
encoder processing circuitry updates the counter, ta (/i) . and calculates the smoothing 
parameter, p(n) . The parameter 0(n) varies between 0.0 and 0.9, being 0.0 for speech, music. 
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tonal-like signals, and non-stationary background noise and ramping up towards 0.9 when 
stationary background noise occurs. 

The LSFs are quantized once per 20 ms frame using a predictive multi-stage vector 
quan- ation. A minimal spacing of 50 Hz is ensured between each two neighboring LSFs before 
quantization. A set of weights is calculated from the LSFs, given by u> = K\P(f,tf* where /, is 
the /* LSF value and />(/. ) is the LPC power spectrum at / ( K is an irrelevant multiplicative 
constant). The reciprocal of the power spectrum is obtained by (up to a multiplicative constant): 
(1 - cos(2#, iWicosilrtfi) - cos(2#,)] 2 
(1 + cos(2^)n[cos(2J5r j ) - cos(2#\)] 2 odd i 



and the power of - 0.4 is then calculated using a lookup table and cubic-spline interpolation 
between table entries. 

A vector of mean values is subtracted from the LSFs, and a vector of prediction error 
vector fe is calculated from the mean removed LSFs vector, using a full-matrix AR(2) predictor. 
A single predictor is used for the rates 5.8, 6.65, 8.0, and 1 1 .0 kbps coders, and two sets of 
prediction coefficients are tested as possible predictors for the 4.55 kbps coder. 

The vector of prediction error is quantized using a multi-stage VQ, with multi-surviving 
candidates from each stage to the next stage. The two possible sets of prediction error vectors 
generated for the 4.55 kbps coder are considered as surviving candidates for the first stage. 

The first 4 stages have 64 entries each, and the fifth and last table have 16 entries. The 
first 3 stages are used for the 4.55 kbps coder, the first 4 stages are used for the 5.8. 6.65 and 8.0 
kbps coders, and all 5 stages are used for the 1 1.0 kbps coder. The following table summarizes 
the number of bits used for the quantization of the LSFs for each rate. 
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prediction 


I s ' stage 


2 nd stage 


3 r0 stage 


4 m stage 


5 m stage 


total 


4.55 kbps 


1 


6 


6 


6 






19 


5.8 kbps 


0 


6 


6 


6 


6 




24 


6.65 kbps 


0 


6 


6 


6 


6 




24 


8.0 kbps 


0 


6 


6 


6 


6 




24 


11.0 kbps 


0 


6 


6 


6 


6 


4 


28 



The number of surviving candidates for each stage is summarized in the following table. 





prediction 


Surviving 


surviving 


surviving 


surviving 




candidates 


candidates 


candidates 


candidates 


candidates 




into the 1 st 


from the 


from the 


from the 


from the 




stage 


1 st stage 


2 nd stage 


3 rt stage 


4 th stage 


4.55 kbps 


2 


10 


6 


4 




5.8 kbps 


1 


8 


6 


4 




6.65 kbps 


1 


8 


8 


4 




8.0 kbps 


1 


8 


8 


4 




11.0 kbps 


1 


8 


6 


4 


4 



The quantization in each stage is done by minimizing the weighted distortion measure 
given by: 

The code vector with index which minimizes € k such that € kmm < e k for all k , is chosen to 
represent the prediction/quantization error (fe represents in this equation both the initial 

prediction error to the first stage and the successive quantization error from each stage to the next 
one). 

The final choice of vectors from all of the surviving candidates (and for the 4.55 kbps 
coder - also the predictor) is done at the end, after the last stage is searched, by choosing a 
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combined set of vectors (and predictor) which minimizes the total error. The contribution from 
all of the stages is summed to form the quantized prediction error vector, and the quantized 
prediction error is added to the prediction states and the mean LSFs value to generate the 
quantized LSFs vector. 

For the 4.55 kbps coder, the number of order flips of the LSFs as the result of the 
quantization if counted, and if the number of flips is more than 1 , the LSFs vector is replaced 
with 0.9 (LSFs of previous frame) + 0.1 • (mean LSFs value) . For all the rates, the quantized 
LSFs are ordered and spaced with a minimal spacing of 50 Hz. 

The interpolation of the quantized LSF is performed in the cosine domain in two ways 

depending on the LTP.mode. If the LTP_mode is 0, a linear interpolation between the 

quantized LSF set of the current frame and the quantized LSF set of the previous frame is 

performed to get the LSF set for the first, second and third subframes as: 

q y {n) = 0.75q 4 (n - 1) + 0.25q 4 («) 
q 2 («) = 0.5q 4 (n - 1) + 0.5? 4 (n) 
q, (n) = 0.259 4 (n - 1) + 0.15q 4 (n) 

where q 4 (n - 1) and q 4 (n) are the cosines of the quantized LSF sets of the previous and current 
frames, respectively, and q t (n) , q 2 {n) andq 3 (,n) are the interpolated LSF sets in cosine domain 
for the first, second and third subframes respectively. 

If the LTP_mode is 1, a search of the best interpolation path is performed in order to get 
the interpolated LSF sets. The search is based on a weighted mean absolute difference between a 
reference LSF set rl (n) and the LSF set obtained from LP analysis_2 / («) . The weights w are 
computed as follows: 
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h<0) = (1-/(0))(1-/(1) + /(0)) 
vt<9) = (l-/(9))(l-/(9) + /(8)) 
fori = lto9 

w(i) = (1 - /(/))(! - Min(l(i + 1) - - /(/ - 1))) 

where Min(a.b) returns the smallest of a and b. 

There are four different interpolation paths. For each path, a reference LSF set rq(n) in 
cosine domain is obtained as follows: 
rq(n) = a(k)q A (n) + (1 - a(k))q 4 (n - 1), k = 1 to 4 

S" = {0.4,0.5,0.6, 0.7} for each path respectively. Then the following distance measure is 
computed for each path as: 

D = |r/(n)-/(n)| 7 w 

The path leading to the minimum distance D is chosen and the corresponding reference LSF set 
rq(n) is obtained as : 

rq in) = CC op ,q 4 (n) + (1 - a opt )q 4 (n - 1) 

The interpolated LSF sets in the cosine domain are then given by: 

^(n) = 0.5^ 4 (n-l) + 0.5r^ (n) 
q 2 {n) = rq(ri) 

q y (ri) = 0.5rq(n) + 0.Sq A (n) 

The impulse response, h(n) , of the weighted synthesis filter 
H (z)W(z) = A{z //,)/[ A(z)A(z/Y 2 )) is computed each subframe. This impulse response is 
needed for the search of adaptive and fixed codebooks 257 and 261 . the impulse response 
h(n) is computed by filtering the vector of coefficients of the filter A{z/y x ) extended by zeros 
through the two filters 1 / A(z) and 1 / A(z I y 2 ) . 
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The target signal for the search of the adaptive codebook 257 is usually computed by subtracting 
the zero input response of the weighted synthesis filter H(z)W(z) from the weighted speech 
signal sjn). This operation is performed on a frame basis. An equivalent procedure for 
computing the target signal is the filtering of the LP residual signal r(n) through the 
combination of the synthesis filter 1/ A(z) and the weighting filter W(z) . 

After determining the excitation for the subframe, the initial states of these filters are 
updated by filtering the difference between the LP residual and the excitation. The LP residual is 
given by: 

10 

r(n) = s(n) + a, s(n - i), n = 0, L _ SF - 1 
.=i 

The residual signal r(n) which is needed for finding the target vector is also used in the adaptive 
codebook search to extend the past excitation buffer. This simplifies the adaptive codebook 
search procedure for delays less than the subframe size of 40 samples. 

In the present embodiment, there are two ways to produce an LTP contribution. One uses 
pitch preprocessing (PP) when the PP-mode is selected, and another is computed like the 
traditional LTP when the LTP-mode is chosen. With the PP-mode, there is no need to do the 
adaptive codebook search, and LTP excitation is directly computed according to past synthesized 
excitation because the interpolated pitch contour is set for each frame. When the AMR coder 
operates with LTP-mode, the pitch lag is constant within one subframe, and searched and coded 
on a subframe basis. 

Suppose the past synthesized excitation is memorized in / ext(MAX_LAG+n), n<0J. 
which is also called adaptive codebook. The LTP excitation codevector, temporally memorized 
in { ext(MAXjAG+n), 0<=n<L_SFJ, is calculated by interpolating the past excitation (adaptive 
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codebook) with the pitch lag contour. T c (n + mL_SF), m = 0.1.2.3. The interpolation is 
performed using an FIR filter (Hamming windowed sine functions): 

ext{MAX _LAG + n) = ^ext{MAX _ LAG + n - T c (n) + i) ■ I,(iJ lc (n)), n = 0,1 L_SF - W 

i—f, 

where Tdn) and T lc (n) are calculated by 
T c (n)=trunc{T e (n + m-L_ SF) } , 
T lc (n)=T c (n)-T c (n), 

m is subframe number. { I s {iJ lc {n)) } is a set of interpolation coefficients. /, is 10. MAXJLAG is 
145+1 1 . and L_SF=40 is the subframe size. Note that the interpolated values 
{ext(MAX_LAG+n), 0<=n<L_SF -17+11} might be used again to do the interpolation when the 
pitch lag is small. Once the interpolation is finished, the adaptive codevector \a=(vjn).n=0 to 
39} is obtained by copying the interpolated values: 
v a (n)=ext(MAX_LAG+n), 0<=n<L_SF 

Adaptive codebook searching is performed on a subframe basis. It consists of performing 
closed-loop pitch lag search, and then computing the adaptive code vector by interpolating the 
past excitation at the selected fractional pitch lag. The LTP parameters (or the adaptive 
codebook parameters) are the pitch lag (or the delay) and gain of the pitch filter. In the search 
stage, the excitation is extended by the LP residual to simplify the closed-loop search. 

For the bit rate of 1 1.0 kbps, the pitch delay is encoded with 9 bits for the 1 st and 3 rd 
subframes and the relative delay of the other subframes is encoded with 6 bits. A fractional pitch 

delay is used in the first and third subframes with resolutions: 1/6 in the range [17,93-], and 

6 

integers only in the range [95.145]. For the second and fourth subframes, a pitch resolution of 
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1/6 is always used for the rate 1 1.0 kbps in the range [7*. -S-.r, +4-], where T x is the pitch 

6 6 

lag of the previous (1 SI or 3 rd ) subframe. 

The close-loop pitch search is performed by minimizing the mean-square weighted error 
between the original and synthesized speech. This is achieved by maximizing the term: 

Yj gs {n)y k (n) 

R(k) = where T 8 s( n ) is * e ^get signal and y k (n) is the past filtered 

V n=0 

excitation at delay k (past excitation convoluted with h(n)). The convolution y k (n) is 
computed for the first delay in the search range, and for the other delays in the search range 
* = 'mm + 1 » it is updated using the recursive relation: 

y k («> s (» - 1) + «(-)*(«) . 

where «(/i),n = -(143+1 1) to 39 is the excitation buffer. 

Note that in the search stage, the samples a(n),n = 0 to 39 , are not available and are 
needed for pitch delays less than 40. To simplify the search, the LP residual is copied to u(n) to 
make the relation in the calculations valid for all delays. Once the optimum integer pitch delay is 
determined, the fractions, as defined above, around that integor are tested. The fractional pitch 
search is performed by interpolating the normalized correlation and searching for its maximum. 

Once the fractional pitch lag is determined, the adaptive codebook vector, v(n) , is 
computed by interpolating the past excitation u(n) at the given phase (fraction). The 
interpolations are performed using two FIR filters (Hamming windowed sine functions), one for 
interpolating the term in the calculations to find the fractional pitch lag and the other for 
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interpolating the past excitation as previously described. The adaptive codebook gain, g p . is 
temporally given then by: 

X 7 « («)>•(«) 



& p 39 



£y(n)y(/i) 

bounded by 0 < g p < 1 .2 , where y(n) = v(n) * A(n) is the filtered adaptive 
codebook vector (zero state response of H(z)W(z) to v(/i) ). The adaptive codebook gain could 
be modified again due to joint optimization of the gains, gain normalization and smoothing. The 
term y(n) is also referred to herein as C p (n) . 

With conventional approaches, pitch lag maximizing correlation might result in two or 
more times the correct one. Thus, with such conventional approaches, the candidate of shorter 
pitch lag is favored by weighting the correlations of different candidates with constant weighting 
coefficients. At times this approach does not correct the double or treble pitch lag because the 
weighting coefficients are not aggressive enough or could result in halving the pitch lag due to 
the strong weighting coefficients. 

In the present embodiment, these weighting coefficients become adaptive by checking if 
the present candidate is in the neighborhood of the previous pitch lags (when the previous frames 
are voiced) and if the candidate of shorter lag is in the neighborhood of the value obtained by 
dividing the longer lag (which maximizes the correlation) with an integer. 

In order to improve the perceptual quality, a speech classifier is used to direct the 
searching procedure of the fixed codebook (as indicated by the blocks 275 and 279) and to- 
control gain normalization (as indicated in the block 401 of Fig. 4). The speech classifier serves 
to improve the background noise performance for the lower rate coders, and to get a quick start- 
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up of the noise level estimation. The speech classifier distinguishes stationary noise-like 
segments from segments of speech, music, tonal-like signals, non-stationary noise, etc. 

The speech classification is performed in two steps. An initial classification 
(speech_mode) is obtained based on the modified input signal. The final classification 
(exc_mode) is obtained from the initial classification and the residual signal after the pitch 
contribution has been removed. The two outputs from the speech classification are the excitation 
mode, exc_mode, and the parameter fi sllb (n) , used to control the subframe based smoothing of the 
gains. 

The speech classification is used to direct the encoder according to the characteristics of 
the input signal and need not be transmitted to the decoder. Thus, the bit allocation, codebooks, 
and decoding remain the same regardless of the classification. The encoder emphasizes the 
perceptually important features of the input signal on a subframe basis by adapting the encoding 
in response to such features. It is important to notice that misclassification will not result in 
disastrous speech quality degradations. Thus, as opposed to the VAD 235, the speech classifier 
identified within the block 279 (Fig. 2) is designed to be somewhat more aggressive for optimal 
perceptual quality. 
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The initial classifier (speechjclassifier) has adaptive thresholds and is performed in six steps: 

1 . Adapt thresholds: 

if(updates_noise > 30 & updates jspeech > 30) 



SNR_max = 3.5 
endif 

if(SNR_max<U5) 
decijmaxjmes = 1.30 
decijnajcp = 0.70 
update jnax_mes = 1.10 
update _majcpjspeech = 0.72 

elseifiSNRjnax < 2.50) 
deci_max_mes = 1.65 
deci_ma_cp = 0.73 
update jmaxjmes = 1.30 
update _majcp_speech = 0.72 

decijmaxjmes = 1 .75 
decijnajcp = 0.77 
update jmaxjmes = 1.30 
update jrnajcpjspeech = 0.77 

2. Calculate parameters: 
Pitch correlation: 




£l(0 s(i-tag) 
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Running mean of pitch correlation: 
m<*-Cp(n) = 0.9 ■ majop(n - 1) + 0. 1 cp 

Maximum of signal amplitude in current pitch cycle: 

max(n) = max{j(i)| f i = start L_SF - 1} 

where: 

start = min{L_SF - lagjO) 
Sum of signal amplitudes in current pitch cycle: 

L_SF-\ 

mean(n)= ^\s(i)\ 

iastart 

Measure of relative maximum: 
max(n) 

max_mes = — 

ma_max_noise{n - 1) 

Maximum to long-term sum: 

maxlsum - |4 max ^ r \l 

^mean(n — k) 

Maximum in groups of 3 subframes for past 15 subframes: 
max^groupin.k) = max<msx(n - 3 • (4 - *) - j\ j = 0 2} k = 0 4 



Group-maximum to minimum of previous 4 group-maxima: 

endmax2minmax = max_group(nA) 

Tmntyiax_group(n f k\ k = 0, ... ,3} 



Slope of 5 group maxima: 

4 

slope = 0. 1 £ (k - 2) • max_group(n y k) 
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3. Classify subframe: 

if(((max_mes < decijmax_mes & ma_cp < deci_ma_cp) I (VAD = 0)) & 

(LTPJ40DE = 1 1 S.Skbit I s 1 4.55kbit I s)) 

speech_mode = 0/* classl * / 
else 

speech_mode = 1 /* classl * / 
endif 

4. Check for change in background noise level, i.e. reset required: 

Check for decrease in level: 

if (updates_noise = 31 & max_mes <= 0.3) 

if (consec.low < 15) 
consec_Iow++ 

endif 
else 

consecjow = 0 
endif 

if (consecjow = 15) 

updates.noise = 0 

lev.reset = - 1 /* low level reset */ 
endif 



Check for increase in level: 

if ((updates.noise >= 30 1 lev.reset = -1) & max_mes > 1.5 & ma_cp < 0.70 & cp < 0.85 
& kl < -0.4 & endmax2minmax < 50 & max2sum < 35 & slope > -100 & slope < 120) 
if (consec.high < 15) 

consec_high++ 
endif 

else 

consec.high = 0 
endif 

if (consecjiigh = 15 & endmax2minmax < 6 & max2sum < 5)) 

updates.noise = 30 

lev.reset = 1 /* high level reset */ 
endif 
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5. Update running mean of maximum of class 1 segments, i.e. stationary noise: 



I * 1 . condition : regular update * / 

(max_mes < update jnax_mes & majcp < 0.6 &cp< 0.65 & max_mes > 0.3) I 

/ * 2. condition : VAD continued update * / 

(consec^vadjO = 8) I 

/ * 3. condition : start - up/reset update * / 

(updates_noise < 30 & ma_cp < 0.7 & cp < 0.75 & *, < -0.4 & endmax2minmax < 5 & 
(levjreset * -1 1 (levjreset = -1 & maxjnes < 2))) 



ma_max_noise{ n ) = 0.9 • majmaxjnoise(n - 1) + 0. 1 - max(n) 

if (updates j^oise < 30) 
updates_noise + + 

levjreset — 0 
where is the first reflection coefficient. 

6. Update running mean of maximum of class 2 segments, i.e. speech, music, tonal-like signals, 
non-stationary noise, etc, continued from above: 

elseif(majcp > update _majcp_speech) 
if (updates jspeech < 80) 



) 



'speech 



= 0.95 



else 



a, 



'speech 



= 0.999 



endif 



majnax_speech(n ) = a % 



'speech 



• majnax_speech ( n - 1 ) + ( 1 - a 



'speech 



) - max(n) 



if (updates_speech £ 80) 
updates^speech + + 
endif 
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The final classifier {exc ^preselect) provides the final class. exc_mode. and the subframe based 
smoothing parameter, fi sub (n) . It has three steps; 

1 . Calculate parameters: 

Maximum amplitude of ideal excitation in current subframe: 
"*a*res2 (") = max{rw2(i )|, 1=0,..., L_SF - 1 } 



Measure of relative maximum: 
maxjmes^ = ret2V 7 



majmax^in-l) 



2. Classify subframe and calculate smoothing: 
if (speech jnode = 1 1 max^mes^ > 1.75) 
excjnode = 1 /* class 2 * / 
= 0 

Njmode_sub(n ) = -4 

excjnode = 0 1* class 1 * / 
N_mode_sub(n) = N_mode_sub(n - 1) + 1 
if(N_rnode_sub(n) > 4) 
N_mode_sub(n) = 4 

if (Njmodejsub(n) > 0) 
/Lb («) = ^ (Njmode_sub(n) - 1 J 2 

= 0 

em/// 
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3. Update running mean of maximum: 
if(max_mes m2 < 0.5) 
if(consec < 51) 

consec + + 
endif 
else 

consec = 0 
endif 

if((exc_mode = 0 & (max_mes ni2 > 0.5 1 consec > 50)) I 

(updates < 30 & mfl.cp < 0.6 &cp< 0.65)) 
ma_max(n) = 0.9 • majmax(n - 1) + 0. 1 - mox^^) 
if (updates < 30) 

updates + + 

em/if 

When this process is completed, the final subframe based classification, exc.mode, and the 
smoothing parameter, 0 SU b(n), are available. 

To enhance the quality of the search of the fixed codebook 261, the target signal, T g (n), is 
produced by temporally reducing the LTP contribution with a gain factor, G r : 
T g (n) = Tgs(n)-G r .*,.Y a (n), n=0,l_.,39 

where T gs (n) is the original target signal 253, YJn) is the filtered signal from the adaptive 

codebook, g p is the LTP gain for the selected adaptive codebook vector, and the gain factor is 

determined according to the normalized LTP gain, Rp, and the bit rate: 

if (rate <=0) /*for 4.45kbps and 5.8kbps*/ 
G r = 0.7R p +0J; 

if (rate = i) /* for 6.65kbps */ 
G r = 0.6 R p +0.4; 
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if (rate ==2 J /* for 8.0kbps V 
C r = 0.3 R p +0.7; 

if (rate==3) /*for 11. 0kbps V 
G r = 0.95; 

if (T„p>L_SF & g p >0.5 & rate<-2) 
G r <^G r (0.3'R p ' +'0.7);and 

where normalized LTP gain, R p , is defined as: 



Another factor considered at the control block 275 in conducting the fixed codebook 
search and at the block 401 (Fig. 4) during gain normalization is the noise level + ")" which is 



where E, is the energy of the current input signal including background noise, and E„ is a 
running average energy of the background noise. E„ is updated only when the input signal is 
detected to be background noise as follows: 

if (first background noise frame is true) 

E„ = 0.75E S ; 
else if (background noise frame is true) 

E„ = 0.75 £„_„ + 0.25 E t ; 

where E„_„ is the last estimation of the background noise energy. 

For each bit rate mode, the fixed codebook 261 (Fig. 2) consists of two or more 
subcodebooks which are constructed with different structure. For example, in the present 
embodiment at higher rates, all the subcodebooks only contain pulses. At lower bit rates, one of 



39 



2X(n) Y a (n) 




given by: 
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the subcodebooks is populated with Gaussian noise. For the lower bit-rates (e.g., 6.65. 5.8. 4.55 
kbps), the speech classifier forces the encoder to choose from the Gaussian subcodebook in case 
of stationary noise-like subframes, exc_mode = 0. For exc_mode = 1 all subcodebooks are 
searched using adaptive weighting. 

For the pulse subcodebooks, a fast searching approach is used to choose a subcodebook 
and select the code word for the current subframe. The same searching routine is used for all the 
bit rate modes with different input parameters. 

In particular, the long-term enhancement filter, F p (z), is used to filter through the selected 
pulse excitation. The filter is defined as F p (z) = ^ _ ^ T y where T is the integer part of 

pitch lag at the center of the current subframe, and fi is the pitch gain of previous subframe, 
bounded by [0.2, 1.0]. Prior to the codebook search, the impulsive response h(n) includes the 
filter F/z). 

For the Gaussian subcodebooks, a special structure is used in order to bring down the 
storage requirement and the computational complexity. Furthermore, no pitch enhancement is 
applied to the Gaussian subcodebooks. 

There are two kinds of pulse subcodebooks in the present AMR coder embodiment. All 
pulses have the amplitudes of +1 or - 1 . Each pulse has 0, 1 , 2, 3 or 4 bits to code the pulse 
position. The signs of some pulses are transmitted to the decoder with one bit coding one sign. 
The signs of other pulses are determined in a way related to the coded signs and their pulse 
positions. 

In the first kind of pulse subcodebook, each pulse has 3 or 4 bits to code the pulse 
position. The possible locations of individual pulses are defined by two basic non-regular tracks 
and initial phases: 
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POS(n p .i) = TRACK (m p ,i) + PHAS{n p , phas_mode) . 

where i=0./ 7 or 15 (corresponding to 3 or 4 bits to code the position), is the possible position 

index, n p = 0,...,N P -1 (N p is the total number of pulses), distinguishes different pulses. m p ~0or I. 

defines two tracks, and phase_mode-0 or 1, specifies two phase modes. 

For 3 bits to code the pulse position, the two basic tracks are: 

/ TRACK(O.i) MO, 4. 8. 12. 18. 24. 30, 36), and 
/ TRACK(l.i) }=(0. 6, 12. 18. 22. 26. 30, 34}. 

If the position of each pulse is coded with 4 bits, the basic tracks are: 

/ TRACK(O.i) MO, 2. 4. 6, 8. 10, 12, 14, 17. 20. 23. 26. 29. 32. 35. 38/, and 
/ TRACKU.i) MO, 3. 6. 9, 12. 15. 18, 21. 23. 25, 27, 29, 31. 33. 35. 37}. 

The initial phase of each pulse is fixed as: 

PHAS(n p ,0) = modulus{n p /MAXPHAS) 
PHAS(n p , 1) = PHAS(N p -\-n p ,0) 

where MAXPHAS is the maximum phase value. 

For any pulse subcodebook, at least the first sign for the first pulse, SIGN(n p ), n p =0, is 
encoded because the gain sign is embedded. Suppose N llgn is the number of pulses with encoded 
signs; that is. SIGN(n„), for n p <A^ i , n ,<=A^, is encoded while SIGN(n p ), for n p >=N ngn , is not 
encoded. Generally, all the signs can be determined in the following way: 
SlGN(n p )= - SIGNinp-DJorn^Nugn, 

due to that the pulse positions are sequentially searched from n p =0 to n p =N p -l using an iteration 
approach. If two pulses are located in the same track while only the sign of the first pulse in the 
track is encoded, the sign of the second pulse depends on its position relative to the first pulse. If 
the position of the second pulse is smaller, then it has opposite sign, otherwise it has the same 
sign as the first pulse. 
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In the second kind of pulse subcodebook, the innovation vector contains 10 signed 
pulses. Each pulse has 0, 1 , or 2 bits to code the pulse position. One subframe with the size of 
40 samples is divided into 10 small segments with the length of 4 samples. 10 pulses are 
respectively located into 10 segments. Since the position of each pulse is limited into one 
segment, the possible locations for the pulse numbered with n p are, f4n p J t {4np, 4n p +2J. or {4n^ 
4n p +L 4n p +2 t 4n p +3 }, respectively for 0, 1, or 2 bits to code the pulse position. All the signs 
for all the 10 pulses are encoded. 

The fixed codebook 261 is searched by minimizing the mean square error between the 
weighted input speech and the weighted synthesized speech. The target signal used for the LTP 
excitation is updated by subtracting the adaptive codebook contribution. That is: 
x 2 ( n)-x(n )- g p y(n), n=0„ . .39, 

where y(n)=v(n)+h(n) is the filtered adaptive codebook vector and g p is the modified 
(reduced) LTP gain. 

If c k is the code vector at index k from the fixed codebook, then the pulse codebook is 
searched by maximizing the term: 

, lt _(C t f,(d'c t f 

where d = H f x } is the correlation between the target signal x^n) and the impulse response 
h(n), H is a the lower triangular Toepliz convolution matrix with diagonal h(6) and lower 

diagonals h(l) A(39) , and <I> = H'H is the matrix of correlations of h(n) . The vector d 

(backward filtered target) and the matrix * are computed prior to the codebook search. The 
elements of the vector d are computed by: 
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39 

d(*)=£ x 2 ( ' Mi-nh n =0,. . . 39. 



and the elements of the symmetric matrix <I> are computed by: 



39 

¥iJ)=^h(n-i)h(n-jl (j>i). 



The correlation in the numerator is given by: 



where m, is the position of the i th pulse and d, is its amplitude. For the complexity reason, all 
the amplitudes { t>, } are set to +1 or -1 ; that is, 

% = SIGNUX i = n p = 0, N p - 1. 

The energy in the denominator is given by: 



To simplify the search procedure, the pulse signs are preset by using the signal b(n ), 
which is a weighted sum of the normalized d(n) vector and the normalized target signal of x 2 (n) 
in the residual domain resrfn): 



If the sign of the i th (i=n p ) pulse located at m, is encoded, it is set to the sign of signal b( n ) at 
that position, i.e., S/GWi>sign[#mJ]. 



i=0 j=i+l 
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In the present embodiment, the fixed codebook 261 has 2 or 3 subcodebooks for each of 
the encoding bit rates. Of course many more might be used in other embodiments. Even with 
several subcodebooks, however, the searching of the fixed codebook 261 is very fast using the 
following procedure. In a first searching turn, the encoder processing circuitry searches the pulse 
positions sequentially from the first pulse (n,=0) to the last pulse (n p =N p -l) by considering the 
influence of all the existing pulses. 

In a second searching turn, the encoder processing circuitry corrects each pulse position 
sequentially from the first pulse to the last pulse by checking the criterion value A k contributed 
from all the pulses for all possible locations of the current pulse. In a third turn, the functionality, 
of the second searching turn is repeated a final time. Of course further turns may be utilized if 
the added complexity is not prohibitive. 

The above searching approach proves very efficient, because only one position of one 
pulse is changed leading to changes in only one term in the criterion numerator C and few terms 
in the criterion denominator E D for each computation of the A*. As an example, suppose a pulse 
subcodebook is constructed with 4 pulses and 3 bits per pulse to encode the position. Only 96 
( 4pulsesx2 i positions per pulsex3turns= 96 ) simplified computations of the criterion A k need be 
performed. 

Moreover, to save the complexity, usually one of the subcodebooks in the fixed codebook 
261 is chosen after finishing the first searching turn. Further searching turns are done only with 
the chosen subcodebook. In other embodiments, one of the subcodebooks might be chosen only 
after the second searching turn or thereafter should processing resources so permit. 

The Gaussian codebook is structured to reduce the storage requirement and the 
computational complexity. A comb-structure with two basis vectors is used. In the comb- 
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structure, the basis vectors are orthogonal, facilitating a low complexity search. In the AMR 

coder, the first basis vector occupies the even sample positions, (0,2 38) , and the second 

basis vector occupies the odd sample positions, (1,3 39) . 

The same codebook is used for both basis vectors, and the length of the codebook vectors 
is 20 samples (half the subframe size). 

All rates (6.65, 5.8 and 4.55 kbps) use the same Gaussian codebook. The Gaussian 
codebook, CB,^, has only 10 entries, and thus the storage requirement is 10-20 = 200 16-bit 
words. From the 10 entries, as many as 32 code vectors are generated. An index, idx 6 , to one 
basis vector 22 populates the corresponding part of a code vector, , in the following way: 

c ia ^(2-(i-T) + 6) = CB Qam$ (U) i = t,t + 1,...,19 
(2 • </ + 20-r) + 5) = CB^dJ) i = 0,1 t - 1 

where the table entry, /, and the shift, r , are calculated from the index, idx s , according to: 

r = truncfydx 6 /10} 
l = idx 6 -\0i 

and 6 is 0 for the first basis vector and 1 for the second basis vector. In addition, a sign is 
applied to each basis vector. 

Basically, each entry in the Gaussian table can produce as many as 20 unique vectors, all 
with the same energy due to the circular shift. The 10 entries are all normalized to have identical 
energy of 0.5, i.e., 

19 

^CB^HJ) 1 =0.5, / = 0,1,...,9 

That means that when both basis vectors have been selected, the combined code vector, c ^ ... 
will have unity energy, and thus the final excitation vector from the Gaussian subcodebook will 
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have unity energy since no pitch enhancement is applied to candidate vectors from the Gaussian 
subcodebook. 

The search of the Gaussian codebook utilizes the structure of the codebook to facilitate a 
low complexity search. Initially, the candidates for the two basis vectors are searched 
independently based on the ideal excitation, res 2 . For each basis vector, the two best candidates, 
along with the respective signs, are found according to the mean squared error. This is 
exemplified by the equations to find the best candidate, index idx a , and its sign, s i1ix : 



where iV Gaim is the number of candidate entries for the basis vector. The remaining parameters 
are explained above. The total number of entries in the Gaussian codebook is 2 2- . The 

fine search minimizes the error between the weighted speech and the weighted synthesized 
speech considering the possible combination of candidates for the two basis vectors from the pre- 



indices k 0 and k, and the respective signs for the two basis vectors, then the final Gaussian code 
vector is selected by maximizing the term: 



over the candidate vectors, d = H r x 2 is the correlation between the target signal x 2 (n) and the 
impulse response h(n) (without the pitch enhancement), and H is a the lower triangular Toepliz 




selection. If c 



is the Gaussian code vector from the candidate vectors represented by the 
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convolution matrix with diagonal h(0) and lower diagonals h(\) A( 39 )- and 4> = H'H is the 

matrix of correlations of h(n). 

More particularly, in the present embodiment, two subcodebooks are included (or 
utilized) in the fixed codebook 261 with 31 bits in the 1 1 kbps encoding mode. In the first 



position. The signs of 6 pulses are transmitted to the decoder with 6 bits. The second 
subcodebook contains innovation vectors comprising 10 pulses. Two bits for each pulse are 
assigned to code the pulse position which is limited in one of the 10 segments. Ten bits are spent 
for 10 signs of the 10 pulses. The bit allocation for the subcodebooks used in the fixed codebook 
261 can be summarized as follows: 

Subcodebook!: 8 pulses X 3 bits/pulse + 6 signs =50 bits 
Subcodebook2: 10 pulses X 2 bits/pulse + 10 signs ~30 bits 

One of the two subcodebooks is chosen at the block 275 (Fig. 2) by favoring the second 

subcodebook using adaptive weighting applied when comparing the criterion value Fl from the 

first subcodebook to the criterion value F2 from the second subcodebook: 

if( W c Fl> F2 ), the first subcodebook is chosen, 
else, the second subcodebook is chosen, 

where the weighting, 0<W C <=7, is defined as: 



the normalized LTP gain, and P sharp is the sharpness parameter of the ideal excitation res 2 (n) 
(i.e., the "sharpness" in the block 279). 



subcodebook, the innovation vector contains 8 pulses. Each pulse has 3 bits to code the pulse 
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In the 8 kbps mode, two subcodebooks are included in the fixed codebook 26 1 with 20 
bits. In the first subcodebook, the innovation vector contains 4 pulses. Each pulse has 4 bits to 
code the pulse position. The signs of 3 pulses are transmitted to the decoder with 3 bits. The 
second subcodebook contains innovation vectors having 10 pulses. One bit for each of 9 pulses 
is assigned to code the pulse position which is limited in one of the 10 segments. Ten bits are 
spent for 10 signs of the 10 pulses. The bit allocation for the subcodebook can be summarized 2 
the following: 

Subcodebookl : 4 pulses X 4 bits/pulse + 3 signs =19 bits 

Subcodebook!: 9 pulses X 1 bits/pulse + 1 pulse X 0 bit + 10 signs =19 bits 

One of the two subcodebooks is chosen by favoring the second subcodebook using adaptive 

weighting applied when comparing the criterion value Fl from the first subcodebook to the 

criterion value F2 from the second subcodebook as in the 1 1 kbps mode. The weighting, 

0<W C <=1, is defined as: 

W, = 1.0-0.6 P„ 5R (L0-0J/? p ) min{/>^+0J, 1.0}. 

The 6.65kbps mode operates using the long-term preprocessing (PP) or the traditional 
LTP. A pulse subcodebook of 18 bits is used when in the PP-mode. A total of 13 bits are 
allocated for three subcodebooks when operating in the LTP-mode. The bit allocation for the 
subcodebooks can be summarized as follows: 
PP-mode: 

Subcodebook: 5 pulses X 3 bits/pulse + 3 signs =18 bits 
LTP-mode: 

Subcodebookl: 3 pulses X 3 bits/pulse + 3 signs =12 bits, phase jmode=l, 
SubcodeboolQ: 3 pulses X 3 bits/pulse + 2 signs =11 bits, phase_mode=0 t 
Subcodebook3: Gaussian subcodebook of 11 bits. 

One of the 3 subcodebooks is chosen by favoring the Gaussian subcodebook when searching 

with LTP-mode. Adaptive weighting is applied when comparing the criterion value from the 
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two pulse subcodebooks to the criterion value from the Gaussian subcodebook. The weighting. 
0< W c < s / , is defined as: 

W c = 1.0 - 0.9 Phsk (1.0 - 05 R p )• min [P slmp + 0.5, 1.0} . 

if (noise-like unvoiced), W C *=W C • (02R p (lO- P^^ + OA). 

The 5.8 kbps encoding mode works only with the long-term preprocessing (PP). Total 14 
bits are allocated for three subcodebooks. The bit allocation for the subcodebooks can be 
summarized as the following: 

Subcodebookl : 4 pulses X 3 bits/pulse + / signs =13 bits, phase _mode= 1 . 
Subcodebook2: 3 pulses X 3 bits/pulse + 3 signs =12 bits, phase _mode=0, 
Subcodebook3: Gaussian subcodebook of 12 bits. 

One of the 3 subcodebooks is chosen favoring the Gaussian subcodebook with aaptive weighting 
applied when comparing the criterion value from the two pulse subcodebooks to the criterion 
value from the Gaussian subcodebook. Tl.j weighting, 0<W C <=1, is defined as: 

W=l.0-P„„(l.0-0.5R p >min{P llmp +0.6,1.0} . 

if (noise - likeunvoiced), W c <=W C (03R p (1 .0 - P )+ 0.7) . 

The 4.55 kbps bit rate mode works only with the long-term preprocessing (PP). Total 10 
bits are allocated for three subcodebooks. The bit allocation for the subcodebooks can be 
summarized as the following: 

Subcodebookl: 2 pulses X 4 bits/pulse + 1 signs -9 bits. phase_mode=l 
Subcodebook2: 2 pulses- X 3 bits/pulse + 2 signs =« bits, phase mode=0, 
Subcodebook3: Gaussian subcodebook of 8 bits. 

One of the 3 subcodebooks is chosen by favoring the Gaussian subcodebook with weighting 

applied when comparing the criterion value from the two pulse subcodebooks to the criterion 

value from the Gaussian subcodebook. The weighting, 0<W C <=1. is defined as: 

W f = 1.0-1^^ (1-0-05 /? p ) min{P jtop +0.6, 1.0}, 
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if (noise - like unvoiced), W c <= W c • (0.6 R p (1.0 - P sharp ) + 0.4) . 

For 4.55, 5.8, 6.65 and 8.0 kbps bit rate encoding modes, a gain re-optimization 
procedure is performed to jointly optimize the adaptive and fixed codebook gains, g and g c 

respectively, as indicated in Fig. 3. The optimal gains are obtained from the following 
correlations given by: 



8„ = 



8c = 



_ R t R 2 — RjR* 
R S R 2 - R 3 R} 



*2 



where *, =< C p ,T t , >, R 2 =< C C ,C C >. R> =< C,,C f >,R A =< C, ,f„ > , and 
R $=<C p ,C p >■ C e ,C p , and T gs are filtered fixed codebook excitation, filtered adaptive 
codebook excitation and the target signal for the adaptive codebook search. 

For 1 1 kbps bit rate encoding, the adaptive codebook gain, g p , remains the same as that 
computed in the closeloop pitch search. The fixed codebook gain, g t , is obtained as: 



where R 6 =< C e ,T, > and T, =T t ,-g p C p . 

Original CELP algorithm is based on the concept of analysis by synthesis (waveform 
matching). At low bit rate or when coding noisy speech, the waveform matching becomes 
difficult so that the gains are up-down, frequently resulting in unnatural sounds. To compensate 
for this problem, the gains obtained in the analysis by synthesis close-loop sometimes need to be 
modified or normalized. 
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There are two basic gain normalization approaches. One is called open-loop approach 
which normalizes the energy of the synthesized excitation to the energy of the unquantized 
residual signal. Another one is close-loop approach with which the normalization is done 
considering the perceptual weighting. The gain normalization factor is a linear combination of 
the one from the close-loop approach and the one from the open-loop approach; the weighting 
coefficients used for the combination are controlled according to the LPC gain. 

The decision to do the gain normalization is made if one of the following conditions is 
met: (a) the bit rate is 8.0 or 6.65 kbps, and noise-like unvoiced speech is true; (b) the noise 
level P NS r is larger than 0.5; (c) the bit rate is 6.65 kbps, and the noise level P NS r is larger than 
0.2; and (d) the bit rate is 5.8 or 4.45kbps. 

The residual energy, £ w , and the target signal energy, E Tgs , are defined respectively as: 



L_SF-l 

E m = 5>* 2 («> 

n=0 



L_SF-\ 



Then the smoothed open-loop energy and the smoothed closed-loop energy are evaluated by: 

if (first subframe is true) 
Ol_Eg = E ns 

else 

Ol_Eg<=P sub Ol_Eg + (l-p sub )E« s 
if (first subframe is true) 

a_Eg=E Tg , 

else 

Cl.Eg^p^a.Eg+O-P^E,,, 
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where p sub is the smoothing coefficient which is determined according to the classification. After 
having the reference energy, the open-loop gain normalization factor is calculated: 



where C ol is 0.8 for the bit rate 1 1.0 kbps, for the other rates C 0 t is 0.7, and v(n) is the excitation: 
v(n) = v a (n) g p + v c (n) g c , n=0 t 1 LJSF-L 

where g p and g c are unquantized gains. Similarly, the closed-loop gain normalization factor is: 



where C ci is 0.9 for the bit rate 1 1.0 kbps, for the other rates C ct is 0.8, and y(n) is the filtered 
signal <y(n)=v(n)*h(n)): 

y(n) = yjn) g p + y c (n) g c , n=0J LJSF-L 

The final gain normalization factor, g f , is a combination of CI __g and Ol controlled in 
terms of an LPC gain parameter, Clpc* 
if (speech is true or the rate is 11 kbps) 

g/ = Clpc Ol_g + (1- Cu>c ) Cl_g 

g f = MAX(L0 t gf) 

g f = MIN(g ft 1+Cu>c) 
if (background noise is true and the rate is smaller than 11kbps) 

g/=L2MINfCU, Ol_gJ 
where Clpc is defined as: 

Cu>c = MIN{sqrt(E r e/E Tg sh 0.8J/0.8 
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Once the gain normalization factor is determined, the unquantized gains are modified: 

For 4.55 ,5.8, 6.65 and 8.0 kbps bit rate encoding, the adaptive codebook gain and the 
fixed codebook gain are vector quantized using 6 bits for rate 4.55 kbps and 7 bits for the other 
rates. The gain codebook search is done by minimizing the mean squared weighted error, Err , 
between the original and reconstructed speech signals: 

Err -K'8 P Cp'g c C c t' 

For rate 1 1.0 kbps, scalar quantization is performed to quantize both the adaptive 
codebook gain, g p , using 4 bits and the fixed codebook gain, g c , using 5 bits each. 

The fixed codebook gain, g c , is obtained by MA prediction of the energy of the scaled 
fixed codebook excitation in the following manner. Let E(n) be the mean removed energy of the 
scaled fixed codebook excitation in (dB) at subframe n be given by: 

E(n) = \0log(±g^c 2 (i))-E. 

where c(i) is the unsealed fixed codebook excitation, and E = 30 dB is the mean energy of 
scaled fixed codebook excitation. 
The predicted energy is given by: 

£(") = XM(*-') 

lot 

where [b x b 2 b,b<] = [0.680.580.340.19] are the MA prediction coefficients and R(n) is the 
quantized prediction error at subframe n . 
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The predicted energy is used to compute a predicted fixed codebook gain g c (by 
substituting E(n) by £(n)and g e by g c ). This is done as follows. First, the mean energy of the 
unsealed fixed codebook excitation is computed as: 

E,=101og(JLfy («)), 
and then the predicted gain g c is obtained as: 

g c = io (005(f< ' ,>+F - £ ' ) . 

A correction factor between the gain, g c . and the estimated one, g c , is given by: 

Y — & c / . 

It is also related to the prediction error as: 
R(n) = E(n) - E(n) = 201ogy . 

The codebook search for 4.55, 5.8, 6.65 and 8.0 kbps encoding bit rates consists of two 
steps. In the first step, a binary search of a single entry table representing the quantized 
prediction error is performed. In the second step, the index Index _ 1 of the optimum entry that is 
closest to the unquantized prediction error in mean square error sense is used to limit the search 
of the two-dimensional VQ table representing the adaptive codebook gain and the prediction 
error. Taking advantage of the particular arrangement and ordering of the VQ table, a fast search 
using few candidates around the entry pointed by Index _ I is performed. In fact, only about half 
of the VQ table entries are tested to lead to the optimum entry with Index _ 2 . Only Index _ 2 is 
transmitted. 
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For 1 1.0 kbps bit rate encoding mode, a full search of both scalar gain codebooks are 
used to quantize g p and g c . For g p , the search is performed by minimizing the error 

Err = abs{g p -g p ). Whereas for g c , the search is performed by minimizing the error 

An update of the states of the synthesis and weighting filters is needed in order to 
compute the target signal for the next subframe. After the two gains are quantized, the excitation 
signal, u(n) , in the present subframe is computed as: 

= g p v(n) + g c c{n\n = 0,39 , 
where g p and g c are the quantized adaptive and fixed codebook gains respectively, v(n) the 
adaptive codebook excitation (interpolated past excitation), and c(n) is the fixed codebook 
excitation. The state of the filters can be updated by filtering the signal r(n)-u(n) through the 
filters 1 / A(z) and W(z) for the 40-sample subframe and saving the states of the filters. This 
would normally require 3 filterings. 

A simpler approach which requires only one filtering is as follows. The local synthesized 
speech at the encoder, s(n) , is computed by filtering the excitation signal through 1/ A(z) . The 
output of the filter due to the input r(n) - u(n) is equivalent to e(n) = s(n) - s(n) , so the states of 
the synthesis filter 1/ A(z) are given by e(n) f n = 0,39. Updating the states of the filter W(z)can 
be done by filtering the error signal e(n) through this filter to find the perceptually weighted 
error e w (n). However, the signal e w (n) can be equivalently found by: 
e w (n) = T gs (n) - g p C p (n) - g c C c (n) . 

The states of the weighting filter are updated by computing e w (n) for n = 30 to 39. 
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The function of the decoder consists of decoding the transmitted parameters (dLP 
parameters, adaptive codebook vector and its gain, fixed codebook vector and its gain) and 
performing synthesis to obtain the reconstructed speech. The reconstructed speech is then 
postfiltered and upscaled. 

The decoding process is performed in the following order. First, the LP filter parameters 
are encoded. The received indices of LSF quantization are used to reconstruct the quantized LSF 
vector. Interpolation is performed to obtain 4 interpolated LSF vectors (corresponding to 4 
subframes). For each subframe, the interpolated LSF vector is convened to LP filter coefficient 
domain, a k , which is used for synthesizing the reconstructed speech in the subframe. 

For rates 4.55, 5.8 and 6.65 (during PP.mode) kbps bit rate encoding modes, the received 
pitch index is used to interpolate the pitch lag across the entire subframe. The following three 
steps are repeated for each subframe: 

1) Decoding of the gains: for bit rates of 4.55, 5.8. 6.65 and 8.0 kbps, the received index is used 
to find the quantized adaptive codebook gain, g p , from the 2-dimensionaI VQ table. The 
same index is used to get the fixed codebook gain correction factor f from the same 
quantization table. The quantized fixed codebook gain, g f . is obtained following these 
steps: 

• the predicted energy is computed £(n) = ]£&,£(« - /) ; 

•si 

• the energy of the unsealed fixed codebook excitation is calculated 
as =101og(4-Tc 2 (i)); and 
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• the predicted gain g e is obtained as g c = l0 {0 05(E(n ^ £ ' E i\ 
The quantized fixed codebook gain is given as g c = fg c . For 1 1 kbps bit rate, the 
received adaptive codebook gain index is used to readily find the quantized adaptive gain. 
g p from the quantization table. The received fixed codebook gain index gives the fixed 
codebook gain correction factor y . The calculation of the quantized fixed codebook 
gain, g c follows the same steps as the other rates. 

2) Decoding of adaptive codebook vector: for 8.0 , 1 1 .0 and 6.65 (during LTP_mode= 1 ) kbps bit 
rate encoding modes, the received pitch index (adaptive codebook index) is used to find the 
integer and fractional parts of the pitch lag. The adaptive codebook v(n) is found by 
interpolating the past excitation w(n) (at the pitch delay) using the FIR filters. 

3) Decoding of fixed codebook vector: the received codebook indices are used to extract the 
type of the codebook (pulse or Gaussian) and either the amplitudes and positions of the 
excitation pulses or the bases and signs of the Gaussian excitation. In either case, the 
reconstructed fixed codebook excitation is given as c(n) . If the integer part of the pitch lag 
is less than the subframe size 40 and the chosen excitation is pulse type, the pitch sharpening 
is applied. This translates into modifying c(n) as c(*) = c(n)+ fic(n- 7), where fi is the 
decoded pitch gain g p from the previous subframe bounded by [0.2,1.0]. 

The excitation at the input of the synthesis filter is given by 
"(*) = 8 p Hn) + g c c(n\n = 0,39 . Before the speech synthesis, a post-processing of the 
excitation elements is performed. This means that the total excitation is modified by 
emphasizing the contribution of the adaptive codebook vector: 
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[«(«) + 0.25 pg v(n), J, >0.5 
I g„ <= 0.5 

Adaptive gain control (AGC) is used to compensate for the gain difference between the 
unemphasized excitation u(n) and emphasized excitation u(n) . The gain scaling factor r/ for 
the emphasized excitation is computed by: 



n = 



£« 2 (n) 

I P >0.5 



1£« 2 («) 

10 J,<=0.5 

The gain-scaled emphasized excitation u(n) is given by: 
w'(«) = n«(w)- 

The reconstructed speech is given by: 

10 

s(n) = u (n) -^ajin - i),n = 0 to 39 , 



where a ( are the interpolated LP filter coefficients. The synthesized speech f (n) is then passed 
through an adaptive postfilter. 

Post-processing consists of two functions: adaptive postfiltering and signal up-scaling. 
The adaptive postfilter is the cascade of three filters: a formant postfilter and two tilt 
compensation filters. The postfilter is updated every subframe of 5 ms. The formant postfilter is 
given by: 
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where A(z) is the received quantized and interpolated LP inverse filter and y„ and Yj control the 
amount of the formant postfiltering. 

The first tilt compensation filter H„(z) compensates for the tilt in the formant postfilter 
H t (z) and is given by: 

where n = y,,*, is a tilt factor, with *, being the first reflection coefficient calculated on the 
truncated impulse response h f («) . of the formant postfilter *, = with- 

r * (0) 

r h (0 = £lt, (j)h, (j + 1) , ( L„ = 22) . 

The postfiltering process is performed as follows. First, the synthesized speech s(n) is 
inverse filtered through ) to produce the residual signal r(n). The signal F(n) is filtered 

by the synthesis filter \/A(zly d ) is passed to the first tilt compensation filter h„ (z) resulting in 
the postfiltered speech signal J, (n) . 

Adaptive gain control (AGC) is used to compensate for the gain difference between the 
synthesized speech signal J(n)and the postfiltered signal s,(n). The gain scaling factor y for 
the present subframe is computed by: 
"15 

v _ "=0 

The gain-scaled postfiltered signal J (n) is given by: 
s (n) = p(n)s,(n) 
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where f}(n) is updated in sample by sample basis and given by: 

where a is an AGC factor with value 0.9. Finally, up-scaling consists of multiplying the 
postfiltered speech by a factor 2 to undo the down scaling by 2 which is applied to the input 
signal. 

Figs. 6 and 7 are drawings of an alternate embodiment of a 4 kbps speech codec that also 
illustrates various aspects of the present invention. In particular, Fig. 6 is a block diagram of a 
speech encoder 601 that is built in accordance with the present invention. The speech encoder 
601 is based on the analysis-by-synthesis principle. To achieve toll quality at 4 kbps, the speech 
encoder 601 departs from the strict waveform-matching criterion of regular CELP coders and 
strives to catch the perceptual important features of the input signal. 

The speech encoder 601 operates on a frame size of 20 ms with three subframes (two of 
6.625 ms and one of 6.75 ms). A look-ahead of 15 ms is used. The one-way coding delay of the 
codec adds up to 55 ms. 

At a block 615, the spectral envelope is represented by a 10 th order LPC analysis for each 
frame. The prediction coefficients are transformed to the Line Spectrum Frequencies (LSFs) for 
quantization. The input signal is modified to better fit the coding model without loss of quality. 
This processing is denoted "signal modification" as indicated by a block 621. In order to 
improve the quality of the reconstructed signal, perceptual important features are estimated and 
emphasized during encoding. 

The excitation signal for an LPC synthesis filter 625 is build from the two traditional 
components: 1) the pitch contribution; and 2) the innovation contribution. The pitch contribution 
is provided through use of an adaptive codebook 627. An innovation codebook 629 has several 
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subcodebooks in order to provide robustness against a wide range of input signals. To each of 
the two contributions a gain is applied which, multiplied with their respective codebook vectors 
and summed, provide the excitation signal. 

The LSFs and pitch lag are coded on a frame basis, and the remaining parameters (the 
innovation codebook index, the pitch gain, and the innovation codebook gain) are coded for 
every subframe. The LSF vector is coded using predictive vector quantization. The pitch lag has 
an integer pan and a fractional pan constituting the pitch period. The quantized pitch period has 
a non-uniform resolution with higher density of quantized values at lower delays. The bit 
allocation for the parameters is shown in the following table. 



Table of Bit Allocation 



Parameter 


Bits per 20 ms 


LSFs 


21 


Pitch lag (adaptive codebook) 


8 


Gains 


12 


Innovation codebook 


3x13 = 39 


Total 


80 



When the quantization of all parameters for a frame is complete the indices are multiplexed to 
form the 80 bits for the serial bit-stream. 

Fig. 7 is a block diagram of a decoder 701 with corresponding functionality to that of the 
encoder of Fig. 6. The decoder 701 receives the 80 bits on a frame basis from a demultiplexer 
711. Upon receipt of the bits, the decoder 701 checks the sync- word for a bad frame indication, 
and decides whether the entire 80 bits should be disregarded and frame erasure concealment 
applied. If the frame is not declared a frame erasure, the 80 bits are mapped to the parameter 
indices of the codec, and the parameters are decoded from the indices using the inverse 
quantization schemes of the encoder of Fig. 6. 
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When the LSFs, pitch lag, pitch gains, innovation vectors, and gains for the innovation 
vectors are decoded, the excitation signal is reconstructed via a block 715. The output signal is 
synthesized by passing the reconstructed excitation signal through an LPC synthesis filter 721. 
To enhance the perceptual quality of the reconstructed signal both short-term and long-term post- 
processing are applied at a block 73 1 . 

Regarding the bit allocation of the 4 kbps codec (as shown in the prior table), the LSFs 
and pitch lag are quantized with 21 and 8 bits per 20 ms, respectively. Although the three 
subframes are of different size the remaining bits are allocated evenly among them. Thus, the 
innovation vector is quantized with 13 bits per subframe. This adds up to a total of 80 bits per 20 
ms, equivalent to 4 kbps. 

The estimated complexity numbers for the proposed 4 kbps codec are listed in the 
following table. All numbers are under the assumption that the codec is implemented on 
commercially available 16-bit fixed point DSPs in full duplex mode. All storage numbers are 
under the assumption of 16-bit words, and the complexity estimates are based on the floating 
point C-source code of the codec. 



Table of Complexity Estimates 



Computational complexity 


30 MIPS 


Program and data ROM 


18 kwords 


RAM 


3 kwords 



The decoder 701 comprises decode processing circuitry that generally operates pursuant 
to software control. Similarly, the encoder 601 (Fig. 6) comprises encoder processing circuitry 
also operating pursuant to software control. Such processing circuitry may coexists, at least in 
part, within a single processing unit such as a single DSP. 
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Fig. 8 is a flow diagram illustrating the functionality of gain normalization such as that 
represented in the block 401 of Fig. 4 by an encoder built in accordance with the present 
invention. In particular, at a block 807, encoding processing circuitry operating pursuant to 
software direction begins normalization processing. At a block 8 1 1 , closed loop gain (CL_G) is 
calculated. Similarly, open loop gain (OL_G) is calculated at a block 815. 

For background noise* as determined at the block 819, a gain normalization factor (G) is 
not limited or smoothed. Instead, for best perceptual quality, the encoding processing circuitry 
selects the lesser of open and closed loop gains to be the normalization factor, at a block 823. 
Thereafter, the adaptive and fixed codebook gains, represented by the blocks 259 and 263 (of 
Fig. 4), are normalized as follows: 

gp* = gp x G, 
gc* = gc x G, 

where the resultant g p ' and gc* comprise the normalized gains applied to adaptive codebook 

vectors and fixed codebook vectors, respectively. 

If the speech signal does not constitute background noise, the encoder processing circuit 
performs a smoothing algorithm at a block 835. The smoothing algorithm is one in which 
contributions from the closed and open loop gains are combined so as to prevent discontinuity 
that might otherwise accompany switching between such gains. Having smoothed the 
contributions of both the closed and open loop gains into a gain normalization factor (G), g p ' and 

gc* are calculated as previously with reference to a block 827. 

The smoothing could be linear or otherwise, and might be based on any of a variety of 
factors. In particular embodiments, the smoothing comprises a linear combination based on a 
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weighting coefficient of linear predictive coding (LPC) gain. The lower the LPC gain, the greater 
the open loop gain contribution to a gain normalization factor. The greater the LPC gain, the 
greater the closed loop gain contribution. 

Fig. 9 is a flow diagram providing a more detailed description of one embodiment of gain 
normalization functionality of Fig. 8. Upon beginning the gain normalization procedure, the 
encoder processing circuitry calculates the open and closed loop gains and the LPC gain 

(^LPC)* at a block 911. Next the encoder processing circuitry determines whether the present 

speech signal comprising background noise as indicated at a block 915. If the signal comprises 
background noise, the gain normalization factor (G) is set to the lesser of the open and close loop 
gains at a block 919. Thereafter, at a block 923, the adaptive and fixed codebook gains are 
normalized as follows: 

gp' = gp x G, 
gc* = gc x G, 

where the resultant g p * and gc' comprise the normalized gains applied to the adaptive and fixed 
codebook vectors, respectively. 

If background noise is not indicated at the block 915, the encoder processing circuitry 
checks to see if the LPC gain is between zero and one at a block 927. If so, the gain 
normalization factor is set to the value of the open loop gain, at a block 93 1 , and applied as 
described above in reference to the block 923. Otherwise, the open and closed loop gains are 
smoothly combined at a block 935. 
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Because the LPC gain changes slowly, a linear combination of the open and closed loop 
gains based on LPC gain weighting yields a slowly changing gain normalization factor. In 
particular, the generation of this gain normalization factor is given by: 

G = (<*=LPC) CLJ3 + (1 - ~u>c) OL_G. 

Thus, when the LPC gain is low for example with unvoiced speech, the open loop gain 
contribution is emphasized. Likewise, when the LPC gain is high, the closed loop gain 
contribution is emphasized. 

After calculating the gain normalization factor at the block 935, the encoder processing 
circuitry limits the gain normalization factor to a predetermined range. Specifically, if the gain 
normalization factor is not greater or equal to one as tested at a block 939, the gain normalization 
factor is set to one at a block 943. With the value of one, no adjustment is made to the adaptive 
and fixed codebook vector gains at the block 923. 

If the gain normalization factor is greater or equal to one but less than one plus the LPC 
gain (as indicated by the blocks 939 and 947), the gain normalization factor need not be limited 
and is therefore applied in the gain normalization process at the block 923. However, if the gain 
normalization factor is not less than one plus the LPC gain, the gain normalization factor is set to 
one plus the LPC gain and delivered to the gain normalization process at the block 923. Such 
limitation is used mainly because the large memory energy of the synthesis filter due to the large 
normalization factor could influence the coding quality of the following subframe when the LPC 
gain is high. 

By applying such smoothed gain normalization, especially at low bit rates for noisy 
speech, the perceptual quality of the synthesized speech is significantly improved. Of course 
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many other ways of smoothly combining the closed and open loop gains and other limits (or 
none at all) could be alternatively or additionally applied. 

Of course, many other modifications and variations are also possible. In view of the 
above detailed description of the present invention and associated drawings, such other 
modifications and variations will now become apparent to those skilled in the art It should also 
be apparent that such other modifications and variations may be effected without departing from 
the spirit and scope of the present invention. 

In addition, the following Appendix A provides a list of many of the definitions, symbols 
and abbreviations used in this application. Appendices B and C respectively provide source and 
channel bit ordering information at various encoding bit rates used in one embodiment of the 
present invention. Appendices A, B and C comprise part of the detailed description of the 
present application, and, otherwise, are hereby incorporated herein by reference in its entirety. 
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APPENDIX A 



For purposes of this application, the following symbols, definitions and abbreviations 



apply. 

adaptive codebook: 



adaptive postfiher: 



Adaptive Multi Rate codec: 



AMR handover: 

channel mode: 

channel mode adaptation: 

channel repacking: 

closed-loop pitch analysis: 



codec mode: 



codec mode adaptation: 



The adaptive codebook contains excitation vectors that are adapted 
for every subframe. The adaptive codebook is derived from the 
long term filter state. The pitch lag value can be viewed as an 
index into the adaptive codebook. 

The adaptive postfilter is applied to the output of the short term 
synthesis filter to enhance the perceptual quality of the 
reconstructed speech. In the adaptive multi-rate codec (AMR), the 
adaptive postfilter is a cascade of two filters: a formant postfilter 
and a tilt compensation filter. 

The adaptive multi-rate code (AMR) is a speech and channel codec 
capable of operating at gross bit-rates of 1 1.4 kbps ("half-rate") 
and 22.8 kbs ( 44 full-rate"). In addition, the codec may operate at 
various combinations of speech and channel coding (codec mode) 
bit-rates for each channel mode. 

Handover between the full rate and half rate channel modes to 
optimize AMR operation. 

Half-rate (HR) or full-rate (FR) operation. 

The control and selection of the (FR or HR) channel mode. 

Repacking of HR (and FR) radio channels of a given radio cell to 
achieve higher capacity within the cell. 

This is the adaptive codebook search, i.e., a process of estimating 
the pitch (lag) value from the weighted input speech and the long 
term filter state. In the closed-loop search, the lag is searched using 
error minimization loop (analysis-by-synthesis). In the adaptive 
multi rate codec, closed-loop pitch search is performed for every 
subframe. 

For a given channel mode, the bit partitioning between the speech 
and channel codecs. 

The control and selection of the codec mode bit-rates. Normally, 
implies no change to the channel mode. 
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direct form coefficients: 



fixed codebook: 



fractional lags: 

full-rate (FR): 
frame: 

gross bit-rate: 
half-rate (HR): 
in-band signaling: 

integer lags: 
interpolating filter: 

inverse filter: 
lag: 

Line Spectral Frequencies: 
Line Spectral Pair: 



One of the formats for storing the short term filter parameters. In 
the adaptive multi rate codec, all filters used to modify speech 
samples use direct form coefficients. 

The fixed codebook contains excitation vectors for speech 
synthesis filters. The contents of the codebook are non-adaptive 
(i.e., fixed). In the adaptive multi rate codec, the fixed codebook 
for a specific rate is implemented using a multi-function codebook. 

A set of lag values having sub-sample resolution. In the adaptive 
multi rate codec a sub-sample resolution between 1/6* and 1.0 of a 
sample is used. 

Full-rate channel or channel mode. 

A time interval equal to 20 ms (160 samples at an 8 kHz sampling 
rate). 

The bit-rate of the channel mode selected (22.8 kbps or 1 1.4 kbps). 
Half-rate channel or channel mode. 

Signaling for DTX, Link Control, Channel and codec mode 
modification, etc. carried within the traffic. 

A set of lag values having whole sample resolution. 

An FIR filter used to produce an estimate of sub-sample resolution 
samples, given an input sampled with integer sample resolution. 

This filter removes the short term correlation from the speech 
signal. The filter models an inverse frequency response of the 
vocal tract. 

The long term filter delay. This is typically the true pitch period, or 
its multiple or sub-multiple. 

(see Line Spectral Pair) 

Transformation of LPC parameters. Line Spectral Pairs are 
obtained by decomposing the inverse filter transfer function A(z) 
to a set of two transfer functions, one having even symmetry and 
the other haying odd symmetry. The Line Spectral Pairs (also 
called as Line Spectra] Frequencies) are the roots of these 
polynomials on the z-unit circle). 
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LP analysis window: 



LP coefficients: 

LTP Mode: 
mode: 

multi-function codebook: 
open-loop pitch search: 



out-of-band signaling: 
PP Mode: 
residual: 

short tenn synthesis filter: 
perceptual weighting filter: 



For each frame, the short term filter coefficients are computed 
using the high pass filtered speech samples within the analysis 
window. In the adaptive multi rate codec, the length of the analysis 
window is always 240 samples. For each frame, two asymmetric 
windows are used to generate two sets of LP coefficient 
coefficients which are interpolated in the LSF domain to construct 
the perceptual weighting filter. Only a single set of LP coefficients 
per frame is quantized and transmitted to the decoder to obtain the 
synthesis filter. A lookahead of 25 samples is used for both HR 
and FR. 

Linear Prediction (LP) coefficients (also referred as Linear 
Predictive Coding (LPC) coefficients) is a generic descriptive term 
for describing the short term filter coefficients. 

Codec works with traditional LTP. 

When used alone, refers to the source codec mode, i.e., to one of 
the source codecs employed in the AMR codec. (See also codec 
mode and channel mode.) 

A fixed codebook consisting of several subcodebooks constructed 
with different kinds of pulse innovation vector structures and noise 
innovation vectors, where codeword from the codebook is used to 
synthesize the excitation vectors. 

A process of estimating the near optimal pitch lag directly from the 
weighted input speech. This is done to simplify the pitch analysis 
and confine the closed-loop pitch search to a small number of lags 
around the open-loop estimated lags. In the adaptive multi rate 
codec, open-loop pitch search is performed once per frame for PP 
mode and twice per frame for LTP mode. 

Signaling on the GSM control channels to support link control. 

Codec works with pitch preprocessing. 

The output signal resulting from an inverse filtering operation. 

This filter introduces, into the excitation signal, short term 
correlation which models the impulse response of the vocal tract. 

This filter is employed in the analysis-by-synthesis search of the 
codebooks. The filter exploits the noise masking properties of the 
formants (vocal tract resonances) by weighting the error less in 
regions near the formant frequencies and more in regions away 
from them. 
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sub frame: 

vector quantization: 
zero input response: 
zero state response: 



A(z) 



A(z) 



1 



B(z) 
W(z) 

F E (z) 
T 



P 

Yn 



A time interval equal to 5-10 ms (40-80 samples at an 8 kHz 
sampling rate). 

A method of grouping several parameters into a vector and 
quantizing them simultaneously. 

The output of a filter due to past inputs, i.e. due to the present state 
of the filter, given that an input of zeros is applied. 

The output of a filter due to the present input, given that no past 
inputs have been applied, i.e., given the state information in the 
filter is all zeroes. 

The inverse filter with unquantized coefficients 
The inverse filter with quantized coefficients 

The speech synthesis filter with quantized coefficients 

The unquantized linear prediction parameters (direct form 
coefficients) 

The quantized linear prediction parameters 
The long-term synthesis filter 

The perceptual weighting filter (unquantized coefficients) 
The perceptual weighting factors 
Adaptive pre-filter 

The nearest integer pitch lag to the closed-loop fractional pitch lag 
of the subframe 

The adaptive pre-filter coefficient (the quantized pitch gain) 
The formant postfilter 

Control coefficient for the amount of the formant post-filtering 
Control coefficient for the amount of the formant post-filtering 
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H,(z) 


Tilt compensation filter 


y, 


Control coefficient for the amount of the tilt compensation filtering 


m = r,*r 


A tilt factor, with being the first reflection coefficient 


h f (n) 


The truncated impulse response of the formant postfilter 


L„ 


The length of h f (n) 




The auto-correlations of hf(n) 


A(z/y ) 


The inverse filter f numerator^ nart nf iKa fnrm ont ns\crfilter 

a itv- ujvviow iiiibi \iiuuieiaiui/ poll \Jl U1C 1UI IIUU1L pOSlIlllCr 


1 / A/? /v , i 


i ne syninesis niier v denominator) part or tne tormant postfilter 




J | "hf» rpcirliial eternal tK#» inuprco filtov A / ^ /v I 

j. tic rcaiuuoj oiguaj or me inverse riiier A( Z/j^/ 




Imoulse reSDOnse of the tilt cofnnencatinn filter 




The AGC-controlled gain scaling factor of the adaptive postfilter 


a 


The AGC factor of the adaptive postfilter 


H ht (z) 


Pre-processing high-pass filter 


W/(n), 


LP analysis windows 


7 (/) 


Length of the first part of the LP analysis window 


^ (/) 


Length of the second part of the LP analysis window w i( n ) 




Length of the first part of the LP analysis window 




Length of the second part of the LP analysis window w //( n ) 




The auto-correlations of the windowed speech 




Lag window for the auto-correlations (60 Hz bandwidth 




expansion) 


/o 


The bandwidth expansion in Hz 
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/, The sampling frequency in Hz 

r' (k) 

ac The modified (bandwidth expanded) auto-correlations 

£/d(0 The prediction error in the ith iteration of the Levinson algorithm 

The ith reflection coefficient 

Oj §) The jth direct form coefficient in the ith iteration of the Levinson 

algorithm 

F\(z) Symmetric LSF polynomial 

^ z ) Antisymmetric LSF polynomial 

F } (z) Polynomial F { (z) with root z = -1 eliminated 

F 2 (z) Polynomial F 2 (*) with root z = 1 eliminated 

* The line spectral pairs (LSFs) in the cosine domain 

q An LSF vector in the cosine domain 

q/ n) The quantized LSF vector at the ith subframe of the frame n 

m * The line spectral frequencies (LSFs) 

T m {x) A mth order Chebyshev polynomial 

f\ (0. h (0 The coefficients of the polynomials F, (z) and F 2 (z) 

f\ (0,/2(0 The coefficients of the polynomials F x '(z) and F 2 \z) 

fit) The coefficients of either F } ( z ) or F 2 (z) 

Qx) Sum polynomial of the Chebyshev polynomials 

x Cosine of angular frequency o 

Recursion coefficients for the Chebyshev polynomial evaluation 

fi The line spectral frequencies (LSFs) in Hz 
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t * = [fx h • • • /io] The vector representation of the LSFs in Hz 

z (l} (n) y z {2) (n) The mean-removed LSF vectors at frame n 

1 V n J , 2 ^ n ) The LSF prediction residual vectors at frame n 

p(n) The predicted LSF vector at frame n 

r (2) (n-\) The quantized second residual vector at the past frame 

The quantized LSF vector at quantization index k 

e lsp The LSF quantization error 

w. , i = 1 10, LSF-quantization weighting factors 

d i The distance between the line spectral frequencies f;+ } and f imm} 

H n ) The impulse response of the weighted synthesis filter 

Ok The correlation maximum of open-loop pitch analysis at delay it 

O t . , i=l t . . . 3 The correlation maxima at delays t h i=l 3 

(Af l ,/,),i = l,...,3 The normalized correlation maxima Af f and the corresponding 

delays /,-,/ = 1 3 

H(z) W(z) = 1 ^[ Z ' r ] ) s The weighted synthesis filter 
A(z)A(z/y 2 ) 

Mz/'ix) The numerator of the perceptual weighting filter 

1 / M z / Y2 ) The denominator of the perceptual weighting filter 

T \ The nearest integer to the fractional pitch lag of the previous (1st 

or 3rd) subframe 

s% (n) The windowed speech signal 

s w( n ) The weighted speech signal 

s (n) Reconstructed speech signal 
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s'(n) The gain-scaled post-filtered signal 

Sf ( n ) Post-filtered speech signal (before scaling) 

x( n ) The target signal for adaptive codebook search 

x 2 (n) x ! 2 The target signal for Fixed codebook search 

resu>(n) The LP residual signal 

c(n) The fixed codebook vector 

v( n ) The adaptive codebook vector 

y(n) = v(n)* h(n) The filtered adaptive codebook vector 

The filtered fixed codebook vector 

y k (n) The past filtered excitation 

wf n ) The excitation signal 

The fully quantized excitation signal 

The gain-scaled emphasized excitation signal 

Top The best open-loop lag 

tmin Minimum lag search value 

'max Maximum lag search value 

Correlation term to be maximized in the adaptive codebook search 

*(*), The interpolated value of R(k) for the integer delay k and fraction 

t 

At Correlation term to be maximized in the algebraic codebook search 

at index k 

Q The correlation in the numerator of A k at index k 

EDk The energy in the denominator of A k at index k 
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The correlation between the target signal x 2 (n) and the impulse 
response h(n), i.e., backward filtered target 

The lower triangular Toepliz convolution matrix with diagonal 
h(0) and lower diagonals A(l),..., h(39) 

The matrix of correlations of h(n) 

The elements of the vector d 

The elements of the symmetric matrix 0 

The innovation vector 

The correlation in the numerator of A k 

The position of the i th pulse 

The amplitude of the i th pulse 

The number of pulses in the fixed codebook excitation 

The energy in the denominator of A k 

The normalized long-term prediction residual 

The sum of the normalized d(n) vector and normalized long-term 
prediction residual res LTP (n) 

The sign signal for the algebraic codebook search 

The fixed codebook vector convolved with H n ) 

The mean-removed innovation energy (in dB) 

The mean of the innovation energy 

The predicted energy 

The MA prediction coefficients 

The quantized prediction error at subframe k 
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The mean innovation energy 

The prediction error of the fixed-codebook gain quantization 
The quantization error of the fixed-codebook gain quantization 

The states of the synthesis filter l/A(z) 

The perceptually weighted error of the analysis-by-synthesis 
search 



tj The gain scaling factor for the emphasized excitation 

g c The fixed-codebook gain 

g c The predicted fixed-codebook gain 

g c The quantized fixed codebook gain 

g p The adaptive codebook gain 

g p The quantized adaptive codebook gain 

If 'gc = 8c * 8c A correction factor between the gain g c and the estimated one g' c 

Y gc The optimum value for y gc 

y sc Gain scaling factor 

AGC Adaptive Gain Control 

AMR Adaptive Multi Rate 

CELP Code Excited Linear Prediction 

C/I Carrier-to-Interferer ratio 

DTX Discontinuous Transmission 

EFR Enhanced Full Rate 

FIR Finite Impulse Response 

FR Full Rate 



-87- 



R(n) 

e{n) 
e w ( n ) 



WO 00/11659 



PCT/US99/19136 



HR 


Half Rate 


LP 


Linear Prediction 


LPC 


Linear Predictive Coding 


LSF 


Line Spectral Frequency 


LSF 


Line Spectral Pair 


LTP 


Long Term Predictor (or Long Term Prediction) 


MA 


Moving Average 


TFO 


Tandem Free Operation 


VAD 


Voice Activity Detection 



-88- 



WO 00/11659 



PCT/US99/19136 



APPENDIX B 

Bit ordering (source coding) 



Bit orderin 


e of output bits from source encoder ( 1 1 kbit/s). 


Bits 


Description 


1-6 


Index of I "LSF stage 


7-12 


Index of 2"° LSF stage 


13-18 


Index of 3 ra LSF stage 


19-24 


Index of 4° LSF stage 


25-28 


Index of 5* LSF stage 


29-32 


Index of adaptive codebook gain. I* subframe 


33-37 


Index of fixed codebook gain, 1" subframe 


38-41 


Index of adaptive codebook gain, 2°° subframe 


42-46 


Index of fixed codebook gain. 2°° subframe 


47-50 


Index of adaptive codebook gain. 3" subframe 


51-55 


Index of fixed codebook gam. 3 m subframe 


56-59 


index of adaptive codebook gain. 4* subframe 


60-64 


Index of fixed codebook gain. 4 th subframe 


65-73 


Index of adaptive codebook. ! u subframe 


74-82 


Index of adaptive codebook. 3" subframe 


83-88 


Index of adaptive codebook (relative). 2°° subframe 


89-94 


Index of adaptive codebook (relative). 4° subframe 


95-96 


index for LSF interpolation 


97-127 


index for fixed codebook. 1** subframe 


128-158 


Index for fixed codebook. 2" subframe 


159-189 


Index for fixed codebook. 3 rd subframe 


190-220 


Index for fixed codebook. 4* subframe 



Bit orderin 


2 of output bits from source encoder (8 kbit/s). 


Bits 


Description 


1-6 


Index of 1 u LSF stage 


7-12 


Index of 2** LSF stage 


13-18 


Index of 3™ LSF stage 


19-24 


Index of 4* LSF stage 


25-31 


Index of fixed and adaptive codebook gains. 1 * subframe 


32-38 


Index of fixed and adaptive codebook gains. 2°* subframe 


39-45 


Index of fixed and adaptive codebook gains. 3 ra subframe 


46-52 


Index of fixed and adaptive codebook gains. 4° subframe 


53-60 


Index of adaptive codebook. 1" subframe 


61-68 


Index of adaptive codebook. 3" subframe 


69-73 


Index of adaptive codebook (relative). 2" 1 subframe 


74-78 


Index of adaptive codebook (relative). 4* subframe 


79-80 


Index for LSF interpolation - 


81-100 


Index for fixed codebook, 1" subframe 


101-120 


Index for fixed codebook. 2* subframe 


121-140 


Index for fixed codebook. 3 rt subframe 


141-160 


Index for fixed codebook. 4* subframe 
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Bit Ofdcnn 


g of output bits from source encoder (6.65 kbit/s). 


Bits 




1-6 




7-12 


Index of 2 - LSF stage " ' 


13-18 




19-24 


index of 4" LSF stage ™ " — 


25-31 




32-38 




39-45 




46-52 




53 


Index for mode (LTP or PP) 


LTPmode 




PPmode 


54-61 








62-69 








70-74 


Index of adaptive codebook {relative), T° subframe 






75-79 








80-81 








82-94 


Index for fixed codebook, 1" subframe 




index for fixed codebook. 1" subframe 


95-107 


Index for fixed codebook. 2*° subframe 




Index for fixed codebook. 2*° subframe 


108-120 


index for fixed codebook. 3 ra subframe 




Index for fixed codebook. 3 W subframe 


121-133 


Index for fixed codebook, 4™ subframe 




Index for fixed codebook. 4 tt subframe 



Bit ordering of output bits from source encoder (5.8 kbit/s). 



Bits 


Description 


1-6 


Index of I" LSF stage 


7-12 


Index of 2*" LSF stage 


13-18 


Index of 3™ LSF stage 


19-24 


Index of 4 m LSF stage 


25-31 


Index of fixed and adaptive codebook gains. 1" subframe 


32-38 


Index of fixed and adaptive codebook sains. 2 M subframe 


39-45 


Index of fixed and adaptive codebook sains. 3 ro subframe 


46-52 


Index of fixed and adaptive codebook gains. 4* subframe 


53-60 


Index of pitch 


61-74 


Index for fixed codebook. 1" subframe 


75-88 


Index for fixed codebook. T° subframe 


89-102 


index for fixed codebook. 3 M subframe 


93-116 


Index for fixed codebook. 4™ subframe 



Bits 


Description 


1-6 


Index of 1" LSF stage 


7-12 


Index of 2* LSF stage 


13-18 


Index of 3 rd LSF stage 


19 


Index of predictor 


20-25 


Index of fixed and adaptive codebook gains. 1 " subframe 


26-31 


Index of fixed and adaptive codebook gains. 2 m subframe 


32-37 


Index of fixed and adaptive codebook sains. 3" subframe 


38-43 


Index of fixed and adaptive codebook gains. 4™ subframe 


44-51 


Index of pitch 


52-61 


Index for fixed codebook. 1 " subframe 


62-71 


Index for fixed codebook. 2* subframe 


72-81 


Index for fixed codebook. 3" subframe 


82-91 


index for fixed codebook. 4* subframe 
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APPENDIX C 



Bit ordering (channel coding) 



Ordering of bits according to subjective importance (11 kbit/s FRTCH) 



Bits, see table XXX 


Descri pti on 


1 


tsfl-0 


2 


Isfl -I 


3 


Isfl -2 


4 


151 1 -j 


5 


lsft-4 




Isfl -5 


7 


lsf2-0 


g 


151 Z~ 1 


9 




10 


Is £2*3 


1 1 




12 


ISIZO 


65 


pncxi 1 -V 


66 


piicn i>i 


67 


piicn i -x 


68 


piicni -j 


69 


pitch 1-4 


70 


pitcn i -j 


74 


pucnj-u 


75 


piicn j- i 


76 


pttch3-2 


77 




78 


pUch3-4 


79 


PHCIL3- J 


29 


ami. ft 

gP' ° 


30 




38 


m2-0 


39 


ptv2-1 


47 


ffD3-0 


48 


Rp3-I 


56 


RP4-0 


57 


Kp4-1 


33 


Rcl-0 


34 


pcl-1 


35 


RCl-2 


42 


RC2-0 


43 


rc2-1 


44 


*c2-2 


51 


gc3-0 


52 


*c3-l 


53 


RC3-2 


60 




61 


JEC4-1 


62 


RC4-2 


71 


pitch 1-6 


72 


pitch 1-7 


73 


pitch 1-8 


80 


piteh3-6 


81 


piich3-7 


82 


pitch3-8 


83 


pitch2-0 


84 


pitch2-l 


85 


pitch2-2 


86 


pitch2-3 


87 


pitch2-4 


88 


_phch2-5 
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89 


pitch4-0 


90 


pitcb4-l 


91 


pttch4-2 


92 


pitcb4-3 


93 


pitch4-4 


94 


pitdi4-5 


13 


lsf3-0 


14 


IsO-l 


15 


UO-2 


16 


UO-3 


17 


IsfM 


18 


lsf3-5 


19 


lsf4-0 


20 


lsf4-l 


21 


bf4-2 


22 


isf4-3 


23 


tsf4-4 


24 


lsf4-5 


25 


lsf5-0 


26 


lsf5-l 


27 


lsf5-2 


28 


bf5-3 


31 


*p1-2 


32 


SPl-3 


40 


KP2-2 


41 


8P2-3 


49 


RP3-2 


50 


RP3-3 


58 


«>4-2 


59 


KP4-3 


36 


RCl-3 


45 


*c2-3 


54 


*c3-3 


63 


«c4-3 


97 


excl-0 


98 


excl-1 


99 


excl-2 


100 


excl-3 


101 


cxcl-4 


102 


cxcl-5 


103 


excl-6 


104 


excl-7 


105 


cxcl-8 


106 


excl-9 


107 


excl-lO 


108 


excMI 


109 


excl-12 


110 


excl-13 


111 


excI-14 


112 


excl-15 


113 


excl-16 


114 


excM7 


115 


excl-18 


116 


excl-19 


117 


exd-20 


118 


excl-21 


119 


cxcl-22 


120 


excl-23 


121 


excl-24 


122 


excl-25 


123 


cxcl-26 


124 


excl-27 


125 


excl-28 


128 


cxc2-0 


129 


exc2-l 
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130 


exc2-2 


131 


exc2-3 


132 


exc2-4 


133 


cxc2-5 


134 


cxc2-6 


135 


exc2-7 


136 


exc2-8 


137 


cxc2-9 


138 


exc2-10 


139 


exc2-ll 


140 


cxc2-l2 


141 


exc2-l3 


142 


exc2-!4 


143 


exc2-15 


144 


exc2-16 


145 


cxc2-l7 


146 


cxc2-18 


147 


cxc2-19 


148 


cxc2-20 


149 


exc2-21 


150 


exc2-22 


1S1 


exc2-23 


152 


cxc2-24 


153 


exc2-25 


154 


cxc2-26 


155 


cxc2-27 


156 


exc2-28 


159 


exc3-0 


160 


cxc3-l 


161 


cxc3-2 


162 


cxc3-3 


163 


cxc3-4 


164 


exc3-5 


165 


exc3-6 


166 


cxc3-7 


167 


cxc3-8 


168 


cxc3-9 


169 


cxc3-!0 


170 


exc3-ll 


171 


cxc3-12 


172 


exc3-l3 


173 


exc3-l4 


174 


cxc3-15 


175 


cxc3-16 


176 


exc3-I7 


177 


cxc3-l8 


178 


cxc3-l9 


179 


cxc3-20 


180 


exc3-21 


181 


exc3-22 


182 


exc3-23 


183 


cxc3-24 


184 


cxc3-25 


185 


exc3-26 


186 


cxc3-27 


187 


cxc3-28 


190 


exc4-0 


191 


cxc4-l 


192 


cxc4-2 


193 


CXC4-3 


194 


CXC4-4 


195 


cxc4-5 


196 


exc4-6 


197 


exc4-7 


198 


cxc4-8 
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199 


CXC4-9 


200 


cxc4-10 


201 


exc4.ll 


202 


exc4-l2 


203 


exc4-13 


204 


exc4-l4 


205 


exc4-!5 


206 


cxc4-I6 


207 


exc4~17 


208 


exc4-18 


209 


exc4-19 


210 


exc4-20 


211 


cxc4-21 


212 


cxc4-22 


213 


exc4-23 


214 


cxc4-24 


215 


cxc4-25 


216 


exc4-26 


217 


exc4-27 


218 


exc4-28 


37 


pcl-4 


46 


f?c2-4 


55 


gc3-4 


64 


JEC4-4 


126 


excl-29 


127 


excl-30 


157 


cxc2-29 


158 


cxc2-30 


188 


cxc3-29 


189 


exc3-30 


219 


exc4-29 


220 


exc4-30 


95 


interp-0 


96 


intefp-1 
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Ordering of bits according to subjective importance (8.0 kbit/s FRTCH) 



Bits, see table XXX 


ipuun 




lsfl-0 


2 


tsfl-1 


3 


Isfl -2 


4 


Isfl -3 


5 


Isfl -4 


6 


1511 "3 


7 




g 


leO 1 


9 


ISI^-Z 


10 




1 ] 


ISI2-4 




lsf2-5 




Kain 1-0 


lit 


gainl-I 




gain I -2 




gain 1-3 




gain! -4 


32 


gain2-0 


55 


gain2-l 


34 


gatn2-2 


35 


gain2-3 


36 


gain2-4 


39 


gain3-0 


40 


gain3-l 


41 


gain3-2 


42 


gain3-3 


jo 


*ain3-4 


4© 


gain4-0 


*» / 


gain4-I 


AS 


gain4-2 


AQ 


gain4-3 




gain4~4 




pitch 1 -0 


54 


pitch 1 - 1 


55 


pitch 1-2 


56 


pitch 1 -3 


57 


pitch 1-4 


58 


pitch 1 -5 


61 


pltCnJ-U 


62 


pitch3-l 


63 




64 


pitch3-3 


65 


pitch3-4 


66 


pitch3-5 


69 


pttch2-0 


70 


pitch2-l 


71 


pitch2-2 


74 


pitcn4-0 


75 


pitch4-l 


76 


pitch4-2 


13 


lsf3-0 


14 


IsO-l 


15 


lsD-2 


16 


lsf3-3 


17 


lsD-4 


18 


lsO-5 


30 


gain 1-5 


37 


gairt2-5 


44 


gain3-5 


51 


gain4-5 


59 


pitch 1-6 


67 


pitch3-6 
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72 


pitcH2-3 


77 




79 


tntt^p-0 


80 


tittcrp- 1 


31 


gun 1-6 


38 


gain2-6 


4* 
^ j 






gain4*6 






20 


lsf4-l 


21 


Uf4-2 


Z2 


lsf4-3 


23 


lsf4-4 


24 


lsf4-5 


60 


pitch 1-7 


68 


pitch3-7 


73 


pstch2-4 


78 


pitch4-4 


81 


excl-0 


82 


excl-1 


83 


exd-2 


84 


excl-3 


85 


excl-4 


86 


exd-S 


87 


excl-6 


88 


excl-7 


89 


excl-8 


90 


excl-9 


91 


excl-10 


92 


excl-11 


93 


excl-12 


94 


excl-13 


95 


cxcl-14 


96 


excl-15 


97 


excl-16 


98 


excl-17 


99 


excl-18 


100 


excl-19 


101 


exc2-0 




cxc2-l 


1 Al 


exc2-2 




exc2-3 


105 


cxc2-4 


106 


cxc2-5 


107 


exc2-6 


108 


exc2-7 


109 


exc2-8 


110 


exc2-9 


111 


exc2-10 


112 


cxc2-ll 


113 


exc2-12 


114 


cxc2-13 


115 


cxc2-14 


116 


cxc2-15 


117 


cxc2-16 


118 


cxc2-17 


119 


cxc2-18 


120 


exc2-19 


121 


cxc3-0 


122 


cxc3-l 


123 


exc3-2 


124 


cxc3-3 


125 


exc3-4 


126 


exc3-5 


127 


exc3-6 
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1) Decoding of the gains: for bit rates of 4.55, 5.8, 6.65 and 8.0 kbps, the received index is 
used to find the quantized adaptive codebook gain, g p , from the 2-dimensional VQ table. 
The same index is used to get the fixed codebook gain correction factor y from the same 
quantization table. The quantized fixed codebook gain, g c , is obtained following these 
steps: 

• the predicted energy is computed E(n) = ^ &,./?(/! - /) ; 

• the energy of the unsealed fixed codebook excitation is calculated 
as ^ ^lOlogC^X^^and 



i=0 



• the predicted gain g c is obtained as g c = 10 1 . 

The quantized fixed codebook gain is given as g c = yg c . For 1 1 kbps bit rate, the 
received adaptive codebook gain index is used to readily find the quantized adaptive 
gain, g p from the quantization table. The received fixed codebook gain index gives 

the fixed codebook gain correction factor y . The calculation of the quantized fixed 
codebook gain, g c follows the same steps as the other rates. 

2) Decoding of adaptive codebook vector: for 8.0 ,1 1 .0 and 6.65 (during LTP_mode=l) kbps 
bit rate encoding modes, the received pitch index (adaptive codebook index) is used to 
find the integer and fractional parts of the pitch lag. The adaptive codebook v(n) is 
found by interpolating the past excitation u(n) (at the pitch delay) using the FIR filters. 

3) Decoding of fixed codebook vector: the received codebook indices are used to extract the 

type of the codebook (pulse or Gaussian) and either the amplitudes and positions of the 

excitation pulses or the bases and signs of the Gaussian excitation. In either case, the 

reconstructed fixed codebook excitation is given as c(n) . If the integer part of the pitch 
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the filters 1/A(z) and W(z) for the 40-sarnple subframe and saving the states of the filters. 
This would normally require 3 filterings. 

A simpler approach which requires only one filtering is as follows. The local 
synthesized speech at the encoder, s(n) , is computed by filtering the excitation signal 
through 1/ A(z) . The output of the filter due to the input r(n) - u(ri) is equivalent to 
e(n) = s(ri) - s (n) , so the states of the synthesis filter 1 / A (z) are given by e(n) 9 n - 0,39. 
Updating the states of the filter W{z) can be done by filtering the error signal e(n) through 
this filter to find the perceptually weighted error e w (n). However, the signal e w (n) can be 
equivalently found by: 
*„(») = T m (n) - g p C p (n) - g c C c (n) . 

The states of the weighting filter are updated by computing e w (n) for n = 30 to 39. 

The function of the decoder consists of decoding the transmitted parameters (dLP 
parameters, adaptive codebook vector and its gain, fixed codebook vector and its gain) and 
performing synthesis to obtain the reconstructed speech. The reconstructed speech is then 
postfiltered and upscaled. 

The decoding process is performed in the following order. First, the LP filter 
parameters are encoded. The received indices of LSF quantization are used to reconstruct the 
quantized LSF vector. Interpolation is performed to obtain 4 interpolated LSF vectors 
(corresponding to 4 sub frames). For each subframe, the interpolated LSF vector is converted 
to LP filter coefficient domain, a k , which is used for synthesizing the reconstructed speech in 
the subframe. 

For rates 4.55, 5.8 and 6.65 (during PP_mode) kbps bit rate encoding modes, the 

received pitch index is used to interpolate the pitch lag across the entire subframe. The 

following three steps are repeated for each subframe: 
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It is also related to the prediction error as: 
R(n) = E(n)-E(n) = 20 log/ . 

The codebook search for 4.55, 5.8, 6.65 and 8.0 kbps encoding bit rates consists of 
two steps. In the first step, a binary search of a single entry table representing the quantized 
prediction error is performed. In the second step, the index Index _\ of the optimum entry 
that is closest to the unquantized prediction error in mean square error sense is used to limit 
the search of the two-dimensional VQ table representing the adaptive codebook gain and the 
prediction error. Taking advantage of the particular arrangement and ordering of the VQ 
table, a fast search using few candidates around the entry pointed by Index _1 is performed. 
In fact, only about half of the VQ table entries are tested to lead to the optimum entry with 
Index _ 2 . Only Index Jl is transmitted. 

For 1 1 .0 kbps bit rate encoding mode, a full search of both scalar gain codebooks are 
used to quantize g p and g c . For g p , the search is performed by minimizing the error 
Err = abs(g p -g p ). Whereas for g c , the search is performed by minimizing the error 

Err -\K-s P c p - gc c c f 

An update of the states of the synthesis and weighting filters is needed in order to 
compute the target signal for the next subframe. After the two gains are quantized, the 
excitation signal, u(n) , in the present subframe is computed as: 

= g p v(*) + g c c(n) y n = 0,39 , 
where g p and g c are the quantized adaptive and fixed codebook gains respectively, v(n) the 
adaptive codebook excitation (interpolated past excitation), and c(«) is the fixed codebook 
excitation. The state of the filters can be updated by filtering the signal r(n) - u(n) through 
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For rate 11.0 kbps, scalar quantization is performed to quantize both the adaptive 



codebook gain, g p , using 4 bits and the fixed codebook gain, g c , using 5 bits each. 

The fixed codebook gain, g c , is obtained by MA prediction of the energy of the 
scaled fixed codebook excitation in the following manner. Let E(ri) be the mean removed 
energy of the scaled fixed codebook excitation in (dB) at subframe n be given by: 



where c(i) is the unsealed fixed codebook excitation, and E = 30 dB is the mean energy of 
scaled fixed codebook excitation. 
The predicted energy is given by: 



where [b&b^ ]= [0.68 0.580.340.19] are the MA prediction coefficients and R(n) is the 
quantized prediction error at subframe n . 

The predicted energy is used to compute a predicted fixed codebook gain g c (by 

substituting E(n) by E(n) and g c by g c ). This is done as follows. First, the mean energy of 
the unsealed fixed codebook excitation is computed as: 



E( n) = 1 0 log(-L g lYc 2 (0) - E , 



4 



l(») = XM("-0 




and then the predicted gain g' c is obtained as: 




A correction factor between the gain, g e , and the estimated one, g' c , is given by: 
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Cl_g = MIN{ C cl 



i 



Cl_Eg L2 

L SF-l * / 

2> 2 («) *' 



ItsO 



where C cl is 0.9 for the bit rate 1 L0 kbps, for the other rates C c i is 0.8, andy(n) is the filtered 
signal (y(n)=v(n)*h(n)): 

y(n) = y a (n) g p + y c (n) g c , n=0J 9 L_SF-L 

The final gain normalization factor, g f9 is a combination of CI jg and Ol _g, controlled 
in terms of an LPC gain parameter, Cu> Ct 
if (speech is true or the rate is 11kbps) 

g/= Clpc Ol_g + (1- C LPC ) CI _g 

g f = MAX(lAgj 

gf = MIN(g f , 1+Cu>c) 
if (background noise is true and the rate is smaller than 11kbps) 

gf=L2 MIN{Cl_g, Ol_g} 
where Cu>c is defined as: 

Cu>c = MIN{sqrt(Eres/E Tg s), 0.8J/0.8 
Once the gain normalization factor is determined, the unquantized gains are modified: 

g P <=g P g f 

For 4.55 ,5.8, 6.65 and 8.0 kbps bit rate encoding, the adaptive codebook gain and the 
fixed codebook gain are vector quantized using 6 bits for rate 4.55 kbps and 7 bits for the 
other rates. The gain codebook search is done by minimizing the mean squared weighted 
error, Err , between the original and reconstructed speech signals: 
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The residual energy, E res , and the target signal energy, E Tgs , are defined respectively 

as: 

L SF-l 



n=0 



L SF-l 



Then the smoothed open-loop energy and the smoothed closed-loop energy are evaluated by: 

if {first subframe is true) 
Ol_Eg=E res 

else 

Ol_Eg<=/? Iui .Ol_Eg + <\-fiJ)E„ 

if (first subframe is true) 
Cl_Eg=E Tgt 

else 

a_Eg<=ft Mb Cl_Eg + (\-fi tub )E Tgf 

where p S ub is the smoothing coefficient which is determined according to the classification. 
After having the reference energy, the open-loop gain normalization factor is calculated: 



i 2y<») *' 

where C a i is 0.8 for the bit rate 11.0 kbps, for the other rates C ol is 0.7, and v(n) is the 
excitation: 

v(n) = v a (n) g p + vjn) g c , n=0, 1 ,L_SF-L 

where g p and g c are unquantized gains. Similarly, the closed-loop gain normalization factor 
is: 
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where R t =<C p ,T gs >, R 2 =<C C ,C C >, R 3 =<C p9 C c >,R 4 =<C e9 T„ >,and 

^5 =< C P > C p > . C f , , and 7^ are filtered fixed codebook excitation, filtered adaptive 
codebook excitation and the target signal for the adaptive codebook search. 

For 1 1 kbps bit rate encoding, the adaptive codebook gain, g p , remains the same as 
that computed in the closeloop pitch search. The fixed codebook gain, g c , is obtained as: 

where R 6 =< C c ,T g > and f g = f v -g p C p . 

Original CELP algorithm is based on the concept of analysis by synthesis (waveform 
matching). At low bit rate or when coding noisy speech, the waveform matching becomes 
difficult so that the gains are up-down, frequently resulting in unnatural sounds. To 
compensate for this problem, the gains obtained in the analysis by synthesis close-loop 
sometimes need to be modified or normalized. 

There are two basic gain normalization approaches. One is called open-loop approach 
which normalizes the energy of the synthesized excitation to the energy of the unquantized 
residual signal. Another one is close-loop approach with which the normalization is done 
considering the perceptual weighting. The gain normalization factor is a linear combination 
of the one from the close-loop approach and the one from the open-loop approach; the 
weighting coefficients used for the combination are controlled according to the LPC gain. 

-^he^iecision to do the gain normalization is made if one of the following conditions is 
met: (a) the bit rate is 8.0 or 6j55 kbps, and noise-like unvoiced speech is true; (b) the noise 
level Pjhk is larger than 0.5; (c) the bit rate is 6;65 kbps, and- the-no^e- level /Wis larger 
than 0.2; and (d) the bit rate is 5.8 or 4.45kbps. 
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One of the 3 subcodebooks is chosen favoring the Gaussian subcodebook with aaptive 
weighting applied when comparing the criterion value from the two pulse subcodebooks to 
the criterion value from the Gaussian subcodebook. The weighting, 0<W C <=1, is defined as: 

^ c =1.0-P^(1.0-0.5^>min{P 5W +0.6,1.0} , 

if {noise - likeunvoiced\ W c <=W C (03R p (1 .0 - P sharp )+ 0.7) . 

The 4.55 kbps bit rate mode works only with the long-term preprocessing (PP). Total 

10 bits are allocated for three subcodebooks. The bit allocation for the subcodebooks can be 

summarized as the following: 

Subcodebook! '; 2 pulses X 4 bits/pulse + 1 signs =9 bits, phase jnode=l 9 
Subcodebook2: 2 pulses X 3 bits/pulse + 2 signs =8 bits, phase _mode=0, 
Subcodebook3: Gaussian subcodebook of 8 bits. 

One of the 3 subcodebooks is chosen by favoring the Gaussian subcodebook with weighting 

applied when comparing the criterion value from the two pulse subcodebooks to the criterion 

value from the Gaussian subcodebook. The weighting, 0<W C <=1 9 is defined as: 

W C =1.0-12P NSR (1.0-05 R p )min{P sAarp +0.6, 1.0}, 

if (noise - like unvoiced), W c <= W c - (0.6 R p (1.0 - P sharp ) + 0.4) . 

For 4.55, 5.8, 6.65 and 8.0 kbps bit rate encoding modes, a gain re-optimization 
procedure is performed to jointly optimize the adaptive and fixed codebook gains, g p and 

g c , respectively, as indicated in Fig. 3. The optimal gains are obtained from the following 
correlations given by: 

P ^5^2 — ^3^3 

8e ' 
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criterion value F2 from the second subcodebook as in the 1 1 kbps mode. The weighting, 
0<W C <=1, is defined as: 

W> 1.0 -0.6 P NSR (LO-OSR^-mmiP^ +0.5, 1.0}. 

The 6.65kbps mode operates using the long-term preprocessing (PP) or the traditional 
LTP. A pulse subcodebook of 18 bits is used when in the PP-mode. A total of 13 bits are 
allocated for three subcodebooks when operating in the LTP-mode. The bit allocation for the 
subcodebooks can be summarized as follows: 
PP-mode: 

Subcodebook: 5 pulses X 3 bits/pulse + 3 signs =18 bits 
LTP-mode: 

Subcodebook J : 3 pulses X 3 bits/pulse + 3 signs =12 bits, phase_mode= 1, 
Subcodebook!: 3 pulses X 3 bits/pulse + 2 signs =11 bits, phase_mode=o[ 
Subcodebook3: Gaussian subcodebook of 11 bits. 

One of the 3 subcodebooks is chosen by favoring the Gaussian subcodebook when searching 

with LTP-mode. Adaptive weighting is applied when comparing the criterion value from the 

two pulse subcodebooks to the criterion value from the Gaussian subcodebook. The 

weighting, 0<W C <=1, is defined as: 

W c = 1.0-0.9/^ (™-05R p yrmn{P siu ^+0J>, 1.0} , 

// (noise - like unvoiced), W C <=W C • (02 R p (1.0 - P sharp ) + 0.8) . 

The 5.8 kbps encoding mode works only with the long-term preprocessing (PP). 

Total 14 bits are allocated for three subcodebooks. The bit allocation for the subcodebooks 

can be summarized as the following: 

Subcodebookl : 4 pulses X 3 bits/pulse + 1 signs =13 bits, phase jnode=l 9 
Subcodebook!: 3 pulses X 3 bits/pulse + 3 signs =12 bits, phase_mode=0, 
Subcodebook3: Gaussian subcodebook of 12 bits. 
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spent for 10 signs of the 10 pulses. The bit allocation for the subcodebooks used in the fixed 

codebook 261 can be summarized as follows: 

Subcodebookl : 8 pulses X 3 bits/pulse + 6 signs =30 bits 
Subcodebook2: 10 pulses X 2 bits/pulse + 10 signs =30 bits 

One of the two subcodebooks is chosen at the block 275 (Fig. 2) by favoring the 

second subcodebook using adaptive weighting applied when comparing the criterion value 

Fl from the first subcodebook to the criterion value F2 from the second subcodebook: 

if (W c * Fl > F2 ), the first subcodebook is chosen, 
else, the second subcodebook is chosen, 

where the weighting, 0<W C <=1, is defined as: 



Pnsr is the background noise to speech signal ratio (i.e., the "noise level" in the block 279), 

R p is the normalized LTP gain, and P sharp is the sharpness parameter of the ideal excitation 

res2(n) (i.e., the "sharpness" in the block 279). 

In the 8 kbps mode, two subcodebooks are included in the fixed codebook 261 with 

20 bits. In the first subcodebook, the innovation vector contains 4 pulses. Each pulse has 4 

bits to code the pulse position. The signs of 3 pulses are transmitted to the decoder with 3 

bits. The second subcodebook contains innovation vectors having 10 pulses. One bit for 

each of 9 pulses is assigned to code the pulse position which is limited in one of the 10 

segments. Ten bits are spent for 10 signs of the 10 pulses. The bit allocation for the 

subcodebook can be summarized as the following: 

Subcodebookl: 4 pulses X 4 bits/pulse + 3 signs =19 bits 

Subcodebook2: 9 pulses X 1 bits/pulse + 1 pulse X 0 bit + 10 signs =19 bits 

One of the two subcodebooks is chosen by favoring the second subcodebook using adaptive 

weighting applied when comparing the criterion value Fl from the first subcodebook to the 




1.0-03 P* 



if Pnsr < 0-5* 
(1.0 - OS R p ) • min {P sharp + 0.5, 1.0} , 
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idx s = 



= max \ 



1=0 



^ = sign|2 res 2 (2 • i + <?) • c i(U , (2 • i + S)j 

where Af Gauss is the number of candidate entries for the basis vector. The remaining 
parameters are explained above. The total number of entries in the Gaussian codebook is 
2 • 2 • AT Gauss 2 . The fine search minimizes the error between the weighted speech and the 
weighted synthesized speech considering the possible combination of candidates for the two 
basis vectors from the pre-selection. If is the Gaussian code vector from the candidate 
vectors represented by the indices k 0 and k x and the respective signs for the two basis 
vectors, then the final Gaussian code vector is selected by maximizing the term: 

A JsJLJiLaJL 

over the candidate vectors, d = H'x 2 is the correlation between the target signal x 2 (n) and 
the impulse response h(n) (without the pitch enhancement), and H is a the lower triangular 
Toepliz convolution matrix with diagonal h(6) and lower diagonals h(l) t ... t h(39), and 

<t> = H / H is the matrix of correlations of h(n) . 

More particularly, in the present embodiment, two subcodebooks are included (or 
utilized) in the fixed codebook 261 with 31 bits in the 1 1 kbps encoding mode. In the first 
subcodebook, the innovation vector contains 8 pulses. Each pulse has 3 bits to code the pulse 
position. The signs of 6 pulses are transmitted to the decoder with 6 bits. The second 
subcodebook contains innovation vectors comprising 10 pulses. Two bits for each pulse are 
assigned to code the pulse position which is limited in one of the 10 segments. Ten bits are 
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one basis vector 22 populates the corresponding part of a code vector, c idx ^ , in the following 
way: 

c idXrf (2 • (/ - r) + S) = CS Gauss (/, /) / = r, r + 1, . . . ,1 9 
c idXtf (2(/ + 20-r)-f^) = C5 Gauss (/,/) / = 0,l,...,r-l 

where the table entry, /, and the shift, r , are calculated from the index, idx s , according to: 

r = trunc^dx 5 /10} 
l = idx s -10-r 

and £ is 0 for the first basis vector and 1 for the second basis vector. In addition, a sign is 
applied to each basis vector. 

Basically, each entry in the Gaussian table can produce as many as 20 unique vectors, 
all with the same energy due to the circular shift. The 10 entries are all normalized to have 
identical energy of 0.5, i.e., 

^LcB Gauss (lJ) 2 =0.5, 7 = 0,1,.. .,9 

That means that when both basis vectors have been selected, the combined code vector, 
c idx 0 4dx l » wiI1 hay e unity energy, and thus the final excitation vector from the Gaussian 
subcodebook will have unity energy since no pitch enhancement is applied to candidate 
vectors from the Gaussian subcodebook. 

The search of the Gaussian codebook utilizes the structure of the codebook to 
facilitate a low complexity search. Initially, the candidates for the two basis vectors are 
searched independently based on the ideal excitation, res 2 . For each basis vector, the two 
best candidates, along with the respective signs, are found according to the mean squared 
error. This is exemplified by the equations to find the best candidate, index idx s , and its sign, 
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contributed from all the pulses for all possible locations of the current pulse. In a third turn, 
the functionality of the second searching turn is repeated a final time. Of course further turns 
may be utilized if the added complexity is not prohibitive. 

The above searching approach proves very efficient, because only one position of one 
pulse is changed leading to changes in only one term in the criterion numerator C and few 
terms in the criterion denominator E D for each computation of the A k . As an example, 
suppose a pulse subcodebook is constructed with 4 pulses and 3 bits per pulse to encode the 
position. Only 96 {A pulses^ positions per pulsex3turns= 96) simplified computations of 
the criterion Ak need be performed. 

Moreover, to save the complexity, usually one of the subcodebooks in the fixed 
codebook 261 is chosen after finishing the first searching turn. Further searching turns are 
done only with the chosen subcodebook. In other embodiments, one of the subcodebooks 
might be chosen only after the second searching turn or thereafter should processing 
resources so permit. 

The Gaussian codebook is structured to reduce the storage requirement and the 
computational complexity. A comb-structure with two basis vectors is used. In the comb- 
structure, the basis vectors are orthogonal, facilitating a low complexity search. In the AMR 
coder, the first basis vector occupies the even sample positions, (0,2,..., 3 8) , and the second 
basis vector occupies the odd sample positions, (1,3, ... ,39) . 

The same codebook is used for both basis vectors, and the length of the codebook 
vectors is 20 samples (half the subframe size). 

All rates (6.65, 5.8 and 4.55 kbps) use the same Gaussian codebook. The Gaussian 
codebook, CJJq^, has only 10 entries, and thus the storage requirement is 10-20 = 200 16- 
bit words. From the 10 entries, as many as 32 code vectors are generated. An index, idx s , to 
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The correlation in the numerator is given by: 



where m i is the position of the i th pulse and i9, is its amplitude. For the complexity reason, 
all the amplitudes { «9, } are set to +1 or -1 ; that is, 
S^SIGNii), i = n p =0 9 .... 9 N p -L . 
The energy in the denominator is given by: 

Ed= E *(«/.»»i) + 2 £ Y*3i&j<f>(m iy mj). 

1=0 i=0 

To simplify the search procedure, the pulse signs are preset by using the signal b(n ), 

which is a weighted sum of the normalized d(n) vector and the normalized target signal of 

X2(n) in the residual domain res2(n)\ 

res 2 (n) 2d(n) 
Uo + > n=0,l,....,39 



[35 + n^ - 

E*w 2 (i>w 2 (0 J2>('V(0 
I /=o V 1=0 



If the sign of the i th pulse located at m t is encoded, it is set to the sign of signal b(n ) 

at that position, i.e., SIGN(i)=sign[b(mt)]. 

In the present embodiment, the fixed codebook 261 has 2 or 3 subcodebooks for each 
of the encoding bit rates. Of course many more might be used in other embodiments. Even 
with several subcodebooks, however, the searching of the fixed codebook 261 is very fast 
using the following procedure. In a first searching turn, the encoder processing circuitry 
searches the pulse positions sequentially from the first pulse (n p =0) to the last pulse (n p =N p - 
7) by considering the influence of all the existing pulses. 

In a second searching turn, the encoder processing circuitry corrects each pulse 
position sequentially from the first pulse to the last pulse by checking the criterion value A k 
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respectively located into 10 segments. Since the position of each pulse is limited into one 
segment, the possible locations for the pulse numbered with n p are, {4n p }, {4n pt 4n p +2), or 
{41%^ 4n p +l f 4n p +2, 4n p +3 }. respectively for 0, 1, or 2 bits to code the pulse position. All 
the signs for all the 10 pulses are encoded. 

The fixed codebook 261 is searched by minimizing the mean square error between the 
weighted input speech and the weighted synthesized speech. The target signal used for the 
LTP excitation is updated by subtracting the adaptive codebook contribution. That is: 
x 2 (n)^x(n)-g p y(n), n=0,...,39, 

where y(n)=v(n)+h(n) is the filtered adaptive codebook vector and g p is the modified 
(reduced) LTP gain. 

If c k is the code vector at index k from the fixed codebook, then the pulse codebook 
is searched by maximizing the term: 



where d = H'x 2 is the correlation between the target signal xjjt) and the impulse response 

A(n), H is a the lower triangular Toepliz convolution matrix with diagonal A(0) and lower 

diagonals A(l) A(39), and <D = H'H is the matrix of correlations of h(n). The vector d 

(backward filtered target) and the matrix <D are computed prior to the codebook search. The 
elements of the vector d are computed by: 
39 

d (»)=^x 2 (i)h(i-n), «=0....39. 
i=n 

and the elements of the symmetric matrix <S> are computed by: 
39 

¥ *'X>= h ( n -0¥n-j), (jTzi). 
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where f=0,7, .... 7 or 15 (corresponding to 3 or 4 bits to code the position), is the possible 

position index, n p = 0, ...,Af p -7 (N p is the total number of pulses), distinguishes different 

pulses, m p ~0 or 7, defines two tracks, and phase jnode—0 or 7, specifies two phase modes. 

For 3 bits to code the pulse position, the two basic tracks are: 

{ TRACK(OJ) }={0, 4, 8. 12, 18, 24, 30, 36}, and 
{TRACK(IJ) }={0, 6, 12, 18, 22, 26, 30, 34). 

If the position of each pulse is coded with 4 bits, the basic tracks are: 

{ TRACK(OJ) MO, 2, 4, 6, 8, 10, 12, 14, 17, 20, 23, 26, 29, 32, 35, 38), and 
{TRACK(IJ) MO, 3, 6, 9, 12, 15, 18, 21, 23, 25, 27, 29, 31, 33, 35, 37). 

The initial phase of each pulse is fixed as: 

PHAS(n p ,0) = modulus(n p /MAXPHAS) 
PHAS(n p ,1) = PHAS(N p - 1 - n p , 0) 

where MAXPHAS is the maximum phase value. 

For any pulse subcodebook, at least the first sign for the first pulse, SIGN(rip) 9 n p =0 9 is 
encoded because the gain sign is embedded. Suppose is the number of pulses with 
encoded signs; that is, SIGN(n p ),for n p <N sign ><=N p , is encoded while SIGN(n p ), for 
n p >=N S ign 9 is not encoded. Generally, all the signs can be determined in the following way: 
SIGN(n p )= - SIGN(n p -l),for n p >=N sign , 

due to that the pulse positions are sequentially searched from n p —0 to n p =N p -l using an 
iteration approach. If two pulses are located in the same track while only the sign of the first 
pulse in the track is encoded, the sign of the second pulse depends on its position relative to 
the first pulse. If the position of the second pulse is smaller, then it has opposite sign, 
otherwise it has the same sign as the first pulse. 

In the second kind of pulse subcodebook, the innovation vector contains 10 signed 
pulses. Each pulse has 0, 1, or 2 bits to code the pulse position. One subframe with the size 
of 40 samples is divided into 10 small segments with the length of 4 samples. 10 pulses are 
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of the subcodebooks is populated with Gaussian noise. For the lower bit-rates (e.g., 6.65, 5.8, 
4.55 kbps), the speech classifier forces the encoder to choose from the Gaussian subcodebook 
in case of stationary noise-like subframes, excjnode = 0. For excjnode 1 all 
subcodebooks are searched using adaptive weighting. 

For the pulse subcodebooks, a fast searching approach is used to choose a 
subcodebook and select the code word for the current subframe. The same searching routine 
is used for all the bit rate modes with different input parameters. 

In particular, the long-term enhancement filter, F p (z), is used to filter through the 

selected pulse excitation. The filter is defined as F p (z) = V where Tis me 

/ \\ P z ) 

integer part of pitch lag at the center of the current subframe, and 0 is the pitch gain of 
previous subframe, bounded by [0.2, 1.0]. Prior to the codebook search, the impulsive 
response h(n) includes the filter F p (z). 

For the Gaussian subcodebooks, a special structure is used in order to bring down the 
storage requirement and the computational complexity. Furthermore, no pitch enhancement is 
applied to the Gaussian subcodebooks. 

There are two kinds of pulse subcodebooks in the present AMR coder embodiment. 
All pulses have the amplitudes of +1 or -1. Each pulse has 0, 1, 2, 3 or 4 bits to code the 
pulse position. The signs of some pulses are transmitted to the decoder with one bit coding 
one sign. The signs of other pulses are determined in a way related to the coded signs and 
their pulse positions. 

In the first kind of pulse subcodebook, each pulse has 3 or 4 bits to code the pulse 
position. The possible locations of individual pulses are defined by two basic non-regular 
tracks and initial phases: 

POS(n p ,i) = TRACK{m p ,/) + PHAS(n p ,phas_mode) , 
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G r = 0.6 R p +0.4; 

if (rate ==2) /*for 8.0kbps */ 
G r = 0.3 R p +0.7; 

if (rate==3) /* for 11. 0kbps V 
G r = 0.95; 

if (T op >L_SF & g p >0.5 & rate<=2) 
G r <=G r (0.3R p +0.7) ;and 

where normalized LTP gain, R p , is defined as: 



Another factor considered at the control block 275 in conducting the fixed codebook 
search and at the block 401 (Fig. 4) during gain normalization is the noise level + ")" which is 
given by: 



where E 5 is the energy of the current input signal including background noise, and E n is a 
running average energy of the background noise. E„ is updated only when the input signal is 
detected to be background noise as follows: 

if (first background noise frame is true) 

E n = 0.75E s ; 
else if (background noise frame is true) 

E„ = 0.75 E njn + 0.25 E s ; 

where E n m is the last estimation of the background noise energy. 

For each bit rate mode, the fixed codebook 261 (Fig. 2) consists of two or more 
subcodebooks which are constructed with different structure. For example, in the present 
embodiment at higher rates, all the subcodebooks only contain pulses. At lower bit rates, one 



39 




n=0 
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3. Update running mean of maximum: 
if(max_mes Ks2 < 0.5) 
if (consec < 51) 

consec + + 
endif 
else 

consec = 0 
endif 



if((exc_mode = 0 & (max^mes^ > 0.5 | consec > 50)) | 
(updates < 30 & ma cp < 0.6 & c/> < 0.65)) 

ma_max(n) = 0.9 • ma_max(n - 1) + 0. 1 - /wax^ («) 

if (updates < 30) 
updates + + 

endif 



When this process is completed, the final subframe based classification, exc_mode, and the 
smoothing parameter, p su b(n), are available. 

To enhance the quality of the search of the fixed codebook 261, the target signal, 
T g (n), is produced by temporally reducing the LTP contribution with a gain factor, G r : 
TgOO = T gs(n) - G r ♦ g p ♦ Y a (n), n=0,l,. . .,39 

where T gs (n) is the original target signal 253, Y a (n) is the filtered signal from the adaptive 

codebook, g p is the LTP gain for the selected adaptive codebook vector, and the gain factor is 

determined according to the normalized LTP gain, Rp, and the bit rate: 

if (rate <=0; /*for 4.45kbps and 5.8kbps*/ 
G r = 0.7R p +0.3; 

if (rate = 1) /* for 6.65kbps */ 
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Classify subframe and calculate smoothing: 

if (speech jnode - 1 1 max_mes ns2 > 1.75) 
excjnode = 1 /* class 2 * / 
Aub(") = 0 
N_modejsub(n ) = -4 

excjnode = 0 1* class 1 * / 
N_mode_sub(n) = N_modejsub{n - 1) + 1 
if(Njnode sub(ri) > 4) 
N_modejsub(n) - 4 

if (N_mode_sub(n) > 0) 

= ^ • (N_mode_sub(n) - if 

else 

Aub(*) = 0 
endif 
endif 
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where k x is the first reflection coefficient. 



6. Update running mean of maximum of class 2 segments, i.e. speech, music, tonal-like 
signals, non-stationary noise, etc, continued from above: 

elseif(ma_cp > update jnajzp_speech) 
if (updates ^speech < 80) 

« speech = 0.95 



"speech = 0*999 



else 
endif 

ma_maxjspeech(n) = a specch . majnax_speech{n - 1) + (1 - a spccch ) . max(n) 

if (updates _speech £ 80) 
updates ^speech + + 
endif 



The final classifier (exc_preseleci) provides the final class, excjnode, and the subframe 
based smoothing parameter, fi sub (n) . It has three steps: 

1. Calculate parameters: 

Maximum amplitude of ideal excitation in current subframe: 
max ft^ (*) = max \res2(i% i = 0, . . . , LJSF - 1} 



Measure of relative maximum: 



max^mes^ = 



majnax^ (n - 1) 
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Check for increase in level: 

if ((updates__noise >= 30 | lev reset = -1) & max_mes > 1.5 & ma cp < 0.70 & cp < 

0.85 

& kl < -0.4 & endmax2minmax < 50 & max2sum < 35 & slope > -100 & slope < 

120) 

if (consec_high < 15) 

consec__high++ 
endif 
else 

consec_high = 0 
endif 

if (consec_high = 15 & endmax2minmax < 6 & max2sum < 5)) 

updates_noise = 30 

lev_reset = 1 /* high level reset */ 
endif 

5. Update running mean of maximum of class 1 segments, i.e. stationary noise: 

n 

I * 1 . condition : regular update * / 

{maxjnes < update jnaxjnes & ma cp < 0.6 & cp < 0.65 & maxjnes > 0.3) | 

/ * 2. condition : VAD continued update * / 

{consec_yad_0 = 8) | 

/ * 3. condition : start - up/reset update * / 

{updates jioise < 30 & majop < 0.7 & cp < 0.75 & k } < -0.4 & endmax2minmax < 5 & 
{levj-eset * -1 1 {levjreset = -1 & maxjnes < 2))) 

) 

majnaxjioise(n ) = 0.9 • ma_max_noise(n - 1) + 0.1 • max(n) 

if {updates jioise < 30) 

updates jtoise + + 
else 

levjreset = 0 
endif 
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Maximum in groups of 3 subframes for past 15 subframes: 
max_group{n, k) = rnzx{max(n - 3 - (4 - k) - j\ j = 0, ... ,2}, k = 0, . . . ,4 



Group-maximum to minimum of previous 4 group-maxima: 

endmax2minmax = , »"x^™<P(»A) 

mm{max_group(n, k) 9 k = 0, „ 

Slope of 5 group maxima: 

4 

slope = 0. 1 • ]T (£ - 2) • max _group(n, k) 

3. Classify subframe: 

if(((max_mes < decijnaxjnes & ma_cp < deci_jnajcp) | (K4Z> = 0)) & 
{LTPJAODE = 1 1 S&kbitls 1 4.55kbit/s)) 
speech_mode = 0 /* c/asyl * / 

speech_mode = 1 /* class! * / 

4. Check for change in background noise level, i.e. reset required: 

Check for decrease in level: 

if (updates_noise = 31 & max_mes <= 0.3) 

if (consec_low< 15) 
consec_low-H- 

endif 
else 

consec_low = 0 
endif 

if (consecjow =15) 

updates_noise = 0 

lev_reset = -1 /* low level reset */ 
endif 
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2. Calculate parameters: 
Pitch correlation: 



L SF-\ 



cp = ,=0 



f L SF-\ \ (L 5F-\ \ 

£*(0*(0 • ^£?{i-lag) s{i-lag) 



Running mean of pitch correlation: 
ma_cp(n) = 0.9 • ma_cp(n — 1) + 0.1 • cp 



Maximum of signal amplitude in current pitch cycle: 
max{n ) = max \s(i )|, i = start, . . . , LJSF - 1 } 

where: 

start = min {Z,_£F - lag ft} 



Sum of signal amplitudes in current pitch cycle: 

L SF-l 

mean(n) = Xi^l 



Measure of relative maximum: 
max(n) 



max mes — 



ma_max_noise{n — 1) 



Maximum to long-term sum: 

- max(n) 
maxzsum = — ^-^ 

^Tmean{n -k) 



*=1 
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The initial classifier {speech jlassifwr) has adaptive thresholds and is performed in six 
steps: 

1 . Adapt thresholds: 

if {updates jioise > 30 & updates jspeech > 30) 

ejv , D . ( ma max speech ^S] 
SNRjnax = mini — = — ,32 

^ majmaxjnoise J 

else 

SNR_max = 3.5 
endif 



ij {SNRjnax <1.75) 
decijnaxjnes = 1.30 
decijnajzp = 0.70 
update jnaxjnes = 1.10 
update jna_cp jspeech = 0.72 

elseif {SNRjnax < 2.50) 
decijnaxjnes = 1.65 
decijnajp = 0.73 
update jnaxjnes = 1 .30 
update jnajpjpeech = 0.72 

decijnaxjnes = 1.75 
decijnajp = 0.77 
update jnaxjnes = 1 .30 
update jnajcpjspeech = 0.77 
endif 
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noise-like segments from segments of speech, music, tonal-like signals, non-stationary noise, 
etc. 

The speech classification is performed in two steps. An initial classification 
{speech jnode) is obtained based on the modified input signal. The final classification 
(exc_mode) is obtained from the initial classification and the residual signal after the pitch 
contribution has been removed. The two outputs from the speech classification are the 
excitation mode, excjnode, and the parameter fi mb {n) , used to control the subframe based 
smoothing of the gains. 

The speech classification is used to direct the encoder according to the characteristics 
of the input signal and need not be transmitted to the decoder. Thus, the bit allocation, 
codebooks, and decoding remain the same regardless of the classification. The encoder 
emphasizes the perceptually important features of the input signal on a subframe basis by 
adapting the encoding in response to such features. It is important to notice that 
misclassification will not result in disastrous speech quality degradations. Thus, as opposed 
to the VAD 235, the speech classifier identified within the block 279 (Fig. 2) is designed to 
be somewhat more aggressive for optimal perceptual quality. 
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interpolating the past excitation as previously described. The adaptive codebook gain, g p9 is 
temporally given then by: 



e = — 

Bp 39 

^y(n)y(n) 



n=0 



bounded by 0 < g p < 1 .2 , where y(n) = v(«) * h(n) is the filtered adaptive 
codebook vector (zero state response of H(z)W(z) to v(n) ). The adaptive codebook gain 
could be modified again due to joint optimization of the gains, gain normalization and 
smoothing. The term y(n) is also referred to herein as C p (n) . 

With conventional approaches, pitch lag maximizing correlation might result in two 
or more times the correct one. Thus, with such conventional approaches, the candidate of 
shorter pitch lag is favored by weighting the correlations of different candidates with constant 
weighting coefficients. At times this approach does not correct the double or treble pitch lag 
because the weighting coefficients are not aggressive enough or could result in halving the 
pitch lag due to the strong weighting coefficients. 

In the present embodiment, these weighting coefficients become adaptive by checking 
if the present candidate is in the neighborhood of the previous pitch lags (when the previous 
frames are voiced) and if the candidate of shorter lag is in the neighborhood of the value 
obtained by dividing the longer lag (which maximizes the correlation) with an integer. 

In order to improve the perceptual quality, a speech classifier is used to direct the 
searching procedure of the fixed codebook (as indicated by the blocks 275 and 279) and to- 
control gain normalization (as indicated in the block 401 of Fig. 4). The speech classifier 
serves to improve the background noise performance for the lower rate coders, and to get a 
quick start-up of the noise level estimation. The speech classifier distinguishes stationary 
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pitch resolution of 1/6 is always used for the rate 11.0 kbps in the range [7] - 5-,T x + 4-] , 

6 6 

where T x is the pitch lag of the previous (1 st or 3 rd ) subframe. 

The close-loop pitch search is performed by minimizing the mean-square weighted 
error between the original and synthesized speech. This is achieved by maximizing the term: 



excitation at delay k (past excitation convoluted with h(n) ). The convolution y k (n) is 
computed for the first delay in the search range, and for the other delays in the search 
range k = + 1,...,/^ , it is updated using the recursive relation: 
yt (*) = y k -x + u(-)H") > 

where «(«), n = -(143 + 1 1) to 39 is the excitation buffer. 

Note that in the search stage, the samples u(n),n = 0 to 39, are not available and are 
needed for pitch delays less than 40. To simplify the search, the LP residual is copied to 
u(ri) to make the relation in the calculations valid for all delays. Once the optimum integer 
pitch delay is determined, the fractions, as defined above, around that integor are tested. The 
fractional pitch search is performed by interpolating the normalized correlation and searching 
for its maximum. 

Once the fractional pitch lag is determined, the adaptive codebook vector, v(«) , is 
computed by interpolating the past excitation u(n) at the given phase (fraction). The 
interpolations are performed using two FIR filters (Hamming windowed sine functions), one 
for interpolating the term in the calculations to find the fractional pitch lag and the other for 



39 



R{k) = 




where T p («) is the target signal and y k («) is the past filtered 
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excitation (adaptive codebook) with the pitch lag contour, r f (n + m-L_SF), m = 0,1,2,3 . The 
interpolation is performed using an FIR filter (Hamming windowed sine functions): 

ext{MAX _LAG + n) = £ ext(MAX _LAG + n-T e (n) + i) • /, (i,T tc (»)), n = 0,1,..., L_SF- l- r 

<=-/, 

where Tc(n) and T lc (n) are calculated by 
T e (n)=trunc{r e (n + m -L_ SF)} , 

m is subframe number, {I s (i,T /c (n))} is a set of interpolation coefficients, f, is 10, 
MAXJAG is 145+1 1, and L_SF=40 is the subframe size. Note that the interpolated values 
{ ext(MAX_LA G+n), 0<=n<L_SF -17+11} might be used again to do the interpolation when 
the pitch lag is small. Once the interpolation is finished, the adaptive codevector 
Va={v a (n),n=0 to 39} is obtained by copying the interpolated values: 
v a (n)=ext(MAX_LAG+n), 0<=n<L_SF 

Adaptive codebook searching is performed on a subframe basis. It consists of 
performing closed-loop pitch lag search, and then computing the adaptive code vector by 
interpolating the past excitation at the selected fractional pitch lag. The LTP parameters (or 
the adaptive codebook parameters) are the pitch lag (or the delay) and gain of the pitch filter. 
In the search stage, the excitation is extended by the LP residual to simplify the closed-loop 
search. 

For the bit rate of 1 1.0 kbps, the pitch delay is encoded with 9 bits for the 1 st and 3 rd 

subframes and the relative delay of the other subframes is encoded with 6 bits. A fractional 

pitch delay is used in the first and third subframes with resolutions: 1/6 in the range 
4 

[17,93-], and integers only in the range [95,145]. For the second and fourth subframes, a 
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The target signal for the search of the adaptive codebook 257 is usually computed by 
subtracting the zero input response of the weighted synthesis filter H(z)W(z) from the 
weighted speech signal s w (n) . This operation is performed on a frame basis. An equivalent 
procedure for computing the target signal is the filtering of the LP residual signal r{n) 
through the combination of the synthesis filter 1/ A(z) and the weighting filter W(z) . 

After determining the excitation for the subframe, the initial states of these filters are 
updated by filtering the difference between the LP residual and the excitation. The LP 
residual is given by: 
10 

= s(n) + a jS (n - /), « = 0, L _ SF - 1 



The residual signal r(/i) which is needed for finding the target vector is also used in the 
adaptive codebook search to extend the past excitation buffer. This simplifies the adaptive 
codebook search procedure for delays less than the subframe size of 40 samples. 

In the present embodiment, there are two ways to produce an LTP contribution. One 
uses pitch preprocessing (PP) when the PP-mode is selected, and another is computed like the 
traditional LTP when the LTP-mode is chosen. With the PP-mode, there is no need to do the 
adaptive codebook search, and LTP excitation is directly computed according to past 
synthesized excitation because the interpolated pitch contour is set for each frame. When the 
AMR coder operates with LTP-mode, the pitch lag is constant within one subframe, and 
searched and coded on a subframe basis. 

Suppose the past synthesized excitation is memorized in { ext(MAX_LA G+n ), n<0), 
which is also called adaptive codebook. The LTP excitation codevector, temporally 
memorized in { ext (MAX_LA G+n), 0<=n<L_SF}, is calculated by interpolating the past 
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There are four different interpolation paths. For each path, a reference LSF set 
rq{ri) in cosine domain is obtained as follows: 
rq(n) = a(k)q 4 («) + (1 - a(k))q 4 (n - 1), k = 1 to 4 

a = {0.4,0.5,0.6, 0.7} for each path respectively. Then the following distance measure is 
computed for each path as: 

r> = |r/(«)-/(«)| r vv 

The path leading to the minimum distance D is chosen and the corresponding reference LSF 
set rq(n) is obtained as : 

rq (») = a opl q 4 («) + (1 - a opt )q 4 (n - 1) 

The interpolated LSF sets in the cosine domain are then given by: 

Qx in) = 0.5q 4 (n - 1) + 0.5rq (n) 
q 2 (n) = rq(n) 

q 3 (n) = 0.5rq(n) + 0.5q A (n) 

The impulse response, h(n) , of the weighted synthesis filter 
H{z)W{z) = A(z I y l )/[A(z)A(z/ y 2 )] is computed each subframe. This impulse response is 
needed for the search of adaptive and fixed codebooks 257 and 261 . The impulse response 
h(n) is computed by filtering the vector of coefficients of the filter A(z I y x ) extended by 
zeros through the two filters 1 / A(z) and 1 / A(z ly 2 ) • 
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quantized prediction error is added to the prediction states and the mean LSFs value to 
generate the quantized LSFs vector. 

For the 4.55 kbps coder, the number of order flips of the LSFs as the result of the 
quantization if counted, and if the number of flips is more than 1, the LSFs vector is replaced 
with 0.9 • (LSFs of previous frame) + 0.1 • (mean LSFs value) . For all the rates, the quantized 
LSFs are ordered and spaced with a minimal spacing of 50 Hz. 

The interpolation of the quantized LSF is performed in the cosine domain in two ways 
depending on the LTP_mode. If the LTP_mode is 0, a linear interpolation between the 
quantized LSF set of the current frame and the quantized LSF set of the previous frame is 
performed to get the LSF set for the first, second and third subframes as: 

q x (jt) = 0J5q A (n - 1) + 0.25? 4 (n) 
q 2 (*) = 0.5? 4 (ti - 1) + 0.5q A (n) 
q 3 (n) = Q.25q 4 (n - 1) 4- 0J5q 4 (n) 

where q 4 (n - 1) and q 4 (n) are the cosines of the quantized LSF sets of the previous and 
current frames, respectively, and q x (n) , q 2 (n) and q 3 (n) are the interpolated LSF sets in 
cosine domain for the first, second and third subframes respectively. 

If the LTP_mode is 1, a search of the best interpolation path is performed in order to 
get the interpolated LSF sets. The search is based on a weighted mean absolute difference 
between a reference LSF set rl (n) and the LSF set obtained from LP analysis_2 l{ti) . The 
weights w are computed as follows: 

m<0) = (1-/(0))(1-/(1) + /(0)) 
u<9) = (l-/(9))(l-/(9) + /(8)) 
^>r/ = lto9 

H</) = (1 - /(/))(! - Min(l(i + 1) - /(/),/(/) - /(/ - 1))) 
where Min(a y b) returns the smallest of a and b. 
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prediction 


1 M stage 


2" u stage 


3 ,u stage 


4"' stage 


5"' stage 


total 


4.55 kbps 


1 


6 


6 


6 






19 


5.8 kbps 


0 


6 


6 


6 


6 




24 


6.65 kbps 


0 


6 


6 


6 


6 




24 


8.0 kbps 


0 


6 


6 


6 


6 




24 


11.0 kbps 


0 


6 


6 


6 


6 


4 


28 



The number of surviving candidates for each stage is summarized in the following table. 



4.55 kbps 


1 prediction 
candidates 
into the 1 st 

1 stage 


Surviving 
candidates 
from the 
1 st stage 


surviving 
candidates 
from the 
2 nd stage 


surviving 
candidates 
from the 
3 rd stage 


surviving 
candidates 
from the 
4 th stage 




10 


6 


4 




5.8 kbps 


1 


8 


6 


4 




6.65 kbps 


1 


8 


8 


4 




8.0 kbps 


1 


8 


8 


4 




11.0 kbps 


1 


8 


6 


4 


4 



The quantization in each stage is done by rnirmnizing the weighted distortion measure 
given by: 



i'aO 

The code vector with index which rnimrnizes s k such that < e k for all k , is chosen 
to represent the prediction/quantization error (fe represents in this equation both the initial 
prediction error to the first stage and the successive quantization error from each stage to the 
next one). 

The final choice of vectors from all of the surviving candidates (and for the 4.55 kbps 
coder - also the predictor) is done at the end, after the last stage is searched, by choosing a 
combined set of vectors (and predictor) which minimizes the total error. The contribution 
from all of the stages is summed to form the quantized prediction error vector, and the 
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The LSFs are quantized once per 20 ms frame using a predictive multi-stage vector 
quantization. A minimal spacing of 50 Hz is ensured between each two neighboring LSFs 

before quantization. A set of weights is calculated from the LSFs, given by w i - K\p(J])f A 
where f t is the /'* LSF value and P{f t ) is the LPC power spectrum at /, ( K is an irrelevant 
multiplicative constant). The reciprocal of the power spectrum is obtained by (up to a 
multiplicative constant): 



p(j,r ~ 



(1 - cos(2<.)n[cos(2^;.) - cos(2^-)] 2 even i 

odd/ 

(1 + cos(2^) n[cos(2^.) - cos(2*f y )] 2 odd i 



even j 



and the power of - 0.4 is then calculated using a lookup table and cubic-spline interpolation 
between table entries. 

A vector of mean values is subtracted from the LSFs, and a vector of prediction error 
vector fe is calculated from the mean removed LSFs vector, using a full-matrix AR(2) 
predictor. A single predictor is used for the rates 5.8, 6.65, 8.0, and 1 1.0 kbps coders, and 
two sets of prediction coefficients are tested as possible predictors for the 4.55 kbps coder. 

The vector of prediction error is quantized using a multi-stage VQ, with multi- 
surviving candidates from each stage to the next stage. The two possible sets of prediction 
error vectors generated for the 4.55 kbps coder are considered as surviving candidates for the 
first stage. 

The first 4 stages have 64 entries each, and the fifth and last table have 16 entries. 
The first 3 stages are used for the 4.55 kbps coder, the first 4 stages are used for the 5.8, 6.65 
and 8.0 kbps coders, and all 5 stages are used for the 1 1 .0 kbps coder. The following table 
summarizes the number of bits used for the quantization of the LSFs for each rate. 
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The parameter J0(n) is controlled by the following logic: 

Stepl: 

if(Vad = 1 1 PastVad = 1 1 k x > 0.5) 
W„»de_ to («-1) = 0 
fi(n) = 0.0 

elseifiN^ ^in - 1) > 0 & (ASP > 0.0015 1 ASP*,, > 0.0024)) 

^ mode . ftm (»-i) = o 

>ff(n) = 0.0 

ekeifiN^ ^in - 1) > 1 & ASP > 0.0025) 
AW™ ("-D = l 

Step 2: 

//(Ha/ = 0 & PastVad = 0) 

^HMxk.fan (») = ^mode.^m (" ~ 0 + 1 

ifiN^^n) > 5) 

Aw)= M.( A ^ niodefiin(w) _ 1)2 

where £, is the first reflection coefficient. 

In step 1, the encoder processing circuitry checks the VAD and the evolution of the 
spectral envelope, and performs a full or partial reset of the smoothing if required. In step 2, 
the encoder processing circuitry updates the counter, N ntode ^(n) , and calculates the 
smoothing parameter, J3(n) . The parameter fi{n) varies between 0.0 and 0.9, being 0.0 for 
speech, music, tonal-like signals, and non-stationary background noise and ramping up 

towards 0.9 when stationary background noise occurs. 
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quantization. As an example, in stationary noise-like signals with constant spectral envelope 
introducing even very small variations in the spectral envelope is picked up easily by the 
human ear and perceived as an annoying modulation. 

The smoothing of the LSFs is done as a running mean according to: 
lsf.(n) = J3{n) - IsMn - 1) + (1 - P(n)) • Isf^est, («), / = 1, . _10 
where Is/jBst^n) is the i th estimated LSF of frame n , and Isf g (n) is the /* LSF for 
quantization of frame n . The parameter p(ri) controls the amount of smoothing, e.g. if 
P{n) is zero no smoothing is applied. 

Pin) is calculated from the VAD information (generated at the block 235) and two 
estimates of the evolution of the spectral envelope. The two estimates of the evolution are 
defined as: 

10 



ASP = - 1)) 2 



«=i 



majsfiin) = fi(n) • majsf^n - 1) + (1 - yff(w)) • Isfjest^n), i = 1,... ,10 



-31- 



WO 00/1 1659 PCT/US99/19136 



= J 0, if ***** T ° pt >M 
T ° pt \r optt otherwise 

The modified weighted speech of the current subframe, memorized in 
{$ w (mO + n),n = 0,1,..., L s - 1 } to update the buffer and produce the second target signal 253 
for searching the fixed codebook 261, is generated by warping the original weighted speech 
{s w (n)} from the original time region, 

[mO + r w , mO + r acc + L s + r opt ] , 
to the modified time region, 
[/wO, /wO+ZJ: 

// 

5 w (m0 + «)= 2^(^0 + /i + 7V(n) + 0 I 5 (UT lw (n)\ n = 0,1,...,Z, - 1, 

where Tw(n) and J/jpfo) are calculated by: 

T w (n) = truncir^ + « - r op/ / L s } , 
7V(w)= + w • r op , / L s - T w (n) , 

{/ x (i, T^in))} is a set of interpolation coefficients. 

After having completed the modification of the weighted speech for the current 
subframe, the modified target weighted speech buffer is updated as follows: 
s w (n) <^$ w (n + L s ) 9 « = 0,l,....,/i m -l. 

The accumulated delay at the end of the current subframe is renewed by: 

^acc ^ T acc + T opt * 

Prior to quantization the LSFs are smoothed in order to improve the perceptual 

quality. In principle, no smoothing is applied during speech and segments with rapid 

variations in the spectral envelope. During non-speech with slow variations in the spectral 

envelope, smoothing is applied to reduce unwanted spectral variations. Unwanted spectral 

variations could typically occur due to the estimation of the LPC parameters and LSF 
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Z.,/2-1 



2k(« + *0 + 4/2)| 



shl (L sr - L s 1 2)max{|s lv (* + »0 + 4 / 2)|, « = 0,1,..., 4 r - L 5 1 2 - 1} 
where n0 = trunc{mO+ r acc + 0.5} (here, ;?7 is subframe number and r acc is the previous 
accumulated delay). 

In order to find the best local delay, r oph at the end of the current processing subframe, 
a normalized correlation vector between the original weighted speech signal and the modified 
matching target is defined as: 



A best local delay in the integer domain, k oph is selected by maximizing Rj(k) in the range of 
k e [SR0,SRI] , which is corresponding to the real delay: 



liRi(k op t)<0.5 y k r is set to zero. 

In order to get a more precise local delay in the range { k r -0. 75+0Jj,j=0J, ...15} 
around k ry R s (k) is interpolated to obtain the fractional correlation vector, Rj(j), by: 



where {I/iJ)} is a set of interpolation coefficients. The optimal fractional delay index, j opi , 
is selected by maximizing R/j). Finally, the best local delay, r opi9 at the end of the current 
processing subframe, is given by, 




K =K P t + n0-m0-r c 



acc 



8 



R f U)= E +/,+') I/(iJ), 7 = 0,1,...,15, 



V=*r-0.75 + 01./ op , 



The local delay is then adjusted by: 
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The target for the modification process of the weighted speech temporally memorized 
in {$ w (m0+ n), n = 0,l,...,Z, 5r - 1 } is calculated by warping the past modified weighted 
speech buffer, s w (m0 +n\ n< 0, with the pitch lag contour, r^ + w-ij, m = 0,1,2 , 

// 

s w (m0 + n) = £s w (/wO + «-:T c (/i) + 0 I s (UT IC {n)\ n = 0,l,...,Z jr - 1, 

'=-// 

where Tc(n) and T IC (n) are calculated by: 
T c (n) =trunc{r c (n + m* L s )} , 

m is subframe number, T IC (n)) is a set of interpolation coefficients, and fi is 10. Then, 

the target for matching, s ( (n\ n = 0,1,..., L sr - 1, is calculated by weighting 

S w (m0 + n), n = 0,1,..., Z jr — 1, in the time domain: 

s t (n) - n s w (m0 + n)/L s , n = 0,1,..., X, -1, 
I («)=5 w (w0 + /!),/! = Z s „...,Z, r -1 

The local integer shifting range SR// for searching for the best local delay is 

computed as the following: 

if speech is unvoiced 
SR0=-l t 
SRJ=l t 

else 

SR0=round{-4min{].0, max{0.0, 1-0.4 (P sh '0.2)}}}, 
SRl=round{4min{L0 p max{0.0 t J-0.4 (P sh -0.2)}}}, 

where P sh =max{P shI , P sh2 J, P s h/ is the average to peak ratio (i.e., sharpness) from the target 
signal: 



/> _ . rr=0 
sh\ — 



I|i(w0 + n)\ 



L sr max$s w (mQ+n% n = 0,l,...,Z,, r - 1} 
and P sh2 is the sharpness from the weighted speech signal: 
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The possible candidates of the precise pitch lag are obtained from the table named as 

PitLagTab8b[i], i=0 t l 127. In the last step, the precise pitch lag P m = PitLagTab8b[I„J is 

possibly modified by checking the accumulated delay r acc due to the modification of the 
speech signal: 

if(r acc >5) l m eziain{I m +l 9 127}, and 
if{r acc <-5) 7 m c=max{/ m -1,0}. 

The precise pitch lag could be modified again: 

'/(^MO) / m <=min{/ m +1, 127}, and 
if(r acc <-10) 7 m <=max{/ m -1,0}. 

The obtained index I m will be sent to the decoder. 

The pitch lag contour, r c (w), is defined using both the current lag P m and the previous 

lagiV/: 

if ( \P m ~ P m -x\ < 0^ min{P m , P^} ) 

r c (") = Prn-i + n(P m -P m _ l )/L /9 n = 0,l,...,Z / -l 
T c (n) = P m ,w-Z/,...,170 

else 

r c (n) =P m .j ,n=0,l,...,39; 
r c (n) = P m ,n=40,...,170 

where Lf= 160 is the frame size. 

One frame is divided into 3 subframes for the long-term preprocessing. For the first 
two subframes, the subframe size, Z, 5 , is 53, and the subframe size for searching, L sr , is 70. 
For the last subframe, L s is 54 and L sr is: 
L sr = min{70, L^L^ -10- r acc ) , 
where Lum=25 is the look-ahead and the maximum of the accumulated delay r acc is limited 
to 14. 
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ejsf = ±-£abs(LSF(i) - LSF _m(i)) , 



if(abs(pit-lagl) < TH and abs(lagj[2> J-lag!) < lagl*02 ) 

if(Rp > 0.5 && pgain _past > 0.7 and ejsf < 05/30 ) LTP_ mode = 0; 
elseZ,7P_mode = l; 

where Rp is current frame normalized pitch correlation, pgain _ past is the quantized pitch 
gain from the fourth subframe of the past frame, TH = MIN(lagl*0\, 5 ) , and 
TH = MAX( 2.0, TH) . 

The estimation of the precise pitch lag at the end of the frame is based on the 
normalized correlation: 



where s w (n + ni) 9 n = 0,1,.... ,£ - 1, represents the last segment of the weighted speech signal 
including the look-ahead ( the look-ahead length is 25 samples), and the size L is defined 
according to the open-loop pitch lag T op with the corresponding normalized correlation C T : 

if (C Tep > 0.6) 
L = max{ 50, T op } 

L = min{ 80. L } 
else 
L = 80 

In the first step, one integer lag k is selected maximizing the Rk in the range 

k e[T op - 10, T op +10] bounded by [17, 145]. Then, the precise pitch lag P m and the 

corresponding index I m for the current frame is searched around the integer lag, [k-1, k+1], 
by up-sampling 



L 



^s w (n + nl)s w (n + nl — k) 




-26- 



WO 00/11659 



PCT/US99/19136 



ki/m+4] 9 m=2 t 3 t 4,5, and if k f > k f 0.95 A / < /, where D is 1.0, 0.85, or 0.65, depending 
on whether the previous frame is unvoiced, the previous frame is voiced and £, is in the 
neighborhood (specified by ± 8) of the previous pitch lag, or the previous two frames are 
voiced and k t is in the neighborhood of the previous two pitch lags. The final selected pitch 
lag is denoted by T op . 

A decision is made every frame to either operate the LTP (long-term prediction) as 
the traditional CELP approach (LTP_mode=l), or as a modified time warping approach 
(LTP_mode=0) herein referred to as PP (pitch preprocessing). For 4.55 and 5.8 kbps 
encoding bit rates, LTP_mode is set to 0 at all times. For 8.0 and 1 1 .0 kbps, LTPjmode is set 
to 1 all of the time. Whereas, for a 6.65 kbps encoding bit rate, the encoder decides whether 
to operate in the LTP or PP mode. During the PP mode, only one pitch lag is transmitted per 
coding frame. 

For 6.65 kbps, the decision algorithm is as follows. First, at the block 241, a 

prediction of the pitch lag pit for the current frame is determined as follows: 

if (LTP_MODE_m = 1 ) 

pit = lagll + 2A*(lagJf3 ]-lagl\ ); 

else 

pit = lagjllj + 2.75*(Iagjr3]-lagJ[lJ); 
where L7T_mode_m is previous frame L7T_mode, lag _f[l] 9 lag _/[3] arethepast 
closed loop pitch lags for second and fourth subframes respectively, lagl is the current frame 
open-loop pitch lag at the second half of the frame, and , lagll is the previous frame open- 
loop pitch lag at the first half of the frame. 

Second, a normalized spectrum difference between the Line Spectrum Frequencies 
(LSF) of current and previous frame is computed as: 
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10 

where Ipc _ gain = Y[ 0 ~ k f ) > where are the reflection coefficients obtained from LP 

m 

analysis_l. 

The voiced/unvoiced decision is derived if the following conditions are met: 

if P2_R\< 0.6 and PI _ SHP > 0.2 set mode = 2, 
ifP3_ZC> 0.4 and PI _ SHP > 0. 1 8 set mode = 2, 
ifPA_RE< 0.4 and PI _ S//P > 0.2 set mode = 2, 
'/ (P2 _ J?l < -L2 + 3.2P1 _ SHP) set VUV = -3 
if(P4_RE < -0.21 + 1A2S6PI_SHP) set VUV = -3 
ST (P3_ ZC > 0.8 - 0.6P1 _ SHP) set VUV = -3 
if {PA RE < 0.1) set VUV = -3 

Open loop pitch analysis is performed once or twice (each 10 ms) per frame 
depending on the coding rate in order to find estimates of the pitch lag at the block 241 (Fig. 
2). It is based on the weighted speech signal s w (n + n m \n = 0,1,....,79, in which n m defines 

the location of this signal on the first half frame or the last half frame. In the first step, four 
maxima of the correlation: 

79 

are found in the four ranges 17.. ..33, 34.. ..67, 68. ...135, 13 6.... 145, respectively. The 
retained maxima , / = 1,2,3,4, are normalized by dividing by: 

„ s w( n m + i = 1.....4, respectively. 

The normalized maxima and corresponding delays are denoted by (R i ,fc i ),i=l J 2,3,4. 

In the second step, a delay, k lt among the four candidates, is selected by maximizing 

the four normalized correlations. In the third step, kj is probably corrected to jfc, (/</) by 

favoring the lower ranges. That is, ifc, (/<7) is selected if*, is within [ki/m-4, 
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A voiced/unvoiced classification and mode decision within the block 279 using the 
input speech s(n) and the residual r w (n) is derived where: 

10 



r w (n) = + * 0, n = 0,L_SF- 1 



The classification is based on four measures: 1) speech sharpness P1_SHP; 2) normalized one 
delay correlation P2_R1; 3) normalized zero-crossing rate P3_ZC; and 4) normalized LP 
residual energy P4_RE. 

The speech sharpness is given by: 

PI £HP = -=2 

MaxL 

where Max is the maximum of abs(r w (n)) over the specified interval of length L . The 
normalized one delay correlation and normalized zero-crossing rate are given by: 

P2 Rl = , "~° 



P3_ZC = J-j? ) sgn[5(z)] - sgn[*(/ - 1)] |], 



where sgn is the sign function whose output is either 1 or -1 depending that the input sample 
is positive or negative. Finally, the normalized LP residual energy is given by: 

P4_RE = 1- ^jlpc _ gain 
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The modified autocorrelations r'(0) = l.OOOlr(O) and r(k) = r(k)w lag (£), k = 1,10 
are used to obtain the reflection coefficients k i and LP filter coefficients a g9 i = 1,10 using 
the Levinson-Durbin algorithm. Furthermore, the LP filter coefficients a f are used to obtain 
the Line Spectral Frequencies (LSFs). 

The interpolated unquantized LP parameters are obtained by interpolating the LSF 
coefficients obtained from the LP analysis_l and those from LP_analysis_2 as: 

q x in) = 0.5? 4 (n - 1) + 0.5q 2 («) 
q 3 (n) = Q.5q 2 (n) + 0.5q 4 (n) 

where q x (n) is the interpolated LSF for subframe 1, q 2 (n) is the LSF of subframe 2 obtained 
from LP_analysis_2 of current frame, q z (n) is the interpolated LSF for subframe 3, 
qA n - 1) is the LSF (cosine domain) from LP_analysis_l of previous frame, and q 4 (n) is the 
LSF for subframe 4 obtained from LP_analysis_l of current frame. The interpolation is 
carried out in the cosine domain. 

A VAD (Voice Activity Detection) algorithm is used to classify input speech frames 
into either active voice or inactive voice frame (background noise or silence) at a block 235 



(Fig. 2). 



The input speech s(n) is used to obtain a weighted speech signal s w (n) by passing 



s(n) through a filter: 




That is, in a subframe of size L_SF, the weighted speech is given by: 



s w (n) = s(n) + ^a.rlsin - i) - J^a^s^n - /), n = 0, L_SF - 1 . 



10 10 
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are performed twice per frame using two different windows. In the first LP analysis 
(LP_anaiysis_l), a hybrid window is used which has its weight concentrated at the fourth 
subframe. The hybrid window consists of two parts. The first part is half a Hamming 
window, and the second part is a quarter of a cosine cycle. The window is given by: 



*>,(«) = 



0.54- 0.46 cos^^ n = 0to214,L = 215 

f0.49(n-L)^ _ _ 

cosl K ' 1, n = 215 to 239 



In the second LP analysis (LP analysis_2), a symmetric Hamming window is used. 



w 2 (n) = 



0.54-0.46cos^y- 



n = 0toll9,L = 120 



0.54 + 0.46cosf (n " L);r \ n = 120 to 239 
V 120 J 



past frame current frame 



future frame 



55 160 25 (samples) 

In either LP analysis, the autocorrelations of the windowed speech s'(n\n = 0,239 are 
computed by: 



239 



r(fy = ^s'(n)s\n-k), k = 0,10. 



n=k 



A 60 Hz bandwidth expansion is used by lag windowing, the autocorrelations using the 
window: 



w ^(0 = exp 



1 ^ 2^60/ V 
2^ 8000 J 



,i = 1,10. 



Moreover, r(0)is multiplied by a white noise correction factor 1.0001 which is equivalent to 
adding a noise floor at -40 dB. 
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The LSF vectors are converted to the LP filter coefficients and interpolated to obtain 
LP filters at each subframe. At each subframe, the decoder processing circuitry constructs 
the excitation signal by: 1) identifying the adaptive and innovative code vectors from the 
codebooks 515 and 519; 2) scaling the contributions by their respective gains at the block 
521; 3) summing the scaled contributions; and 3) modifying and applying adaptive tilt 
compensation at the blocks 527 and 529. The speech signal is also reconstructed on a 
subframe basis by filtering the excitation through the LP synthesis at the block 53 L Finally, 
the speech signal is passed through an adaptive post filter at the block 535 to generate the 
reproduced speech signal 539. 

The AMR encoder will produce the speech modeling information in a unique 
sequence and format, and the AMR decoder receives the same information in the same way. 
The different parameters of the encoded speech and their individual bits have unequal 
importance with respect to subjective quality. Before being submitted to the channel 
encoding function the bits are rearranged in the sequence of importance. 

Two pre-processing functions are applied prior to the encoding process: high-pass 
filtering and signal down-scaling. Down-scaling consists of dividing the input by a factor of 
2 to reduce the possibility of overflows in the fixed point implementation. The high-pass 
filtering at the block 215 (Fig. 2) serves as a precaution against undesired low frequency 
components. A filter with cut off frequency of 80 Hz is used, and it is given by: 

H (z) = 0.92727435 - 1 .85449412'' + 0.92727435z" 2 
1 - 1.90594652 - 1 + 0.91 14024z' 2 

Down scaling and high-pass filtering are combined by dividing the coefficients of the 

numerator of H M (2) by 2. 

Short-term prediction, or linear prediction (LP) analysis is performed twice per speech 
frame using the autocorrelation approach with 30 ms windows. Specifically, two LP analyses 
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code vector) from x(n). The encoder processing circuitry uses the second target signal 253 in 
the fixed codebook search to find the optimum innovation. 

Fifth, for the 11.0 kbps bit rate mode, the gains of the adaptive and fixed codebook 
are scalar quantized with 4 and 5 bits respectively (with moving average prediction applied to 
the fixed codebook gain). For the other modes the gains of the adaptive and fixed codebook 
are vector quantized (with moving average prediction applied to the fixed codebook gain). 

Finally, the filter memories are updated using the determined excitation signal for 
finding the first target signal in the next subframe. 

The bit allocation of the AMR codec modes is shown in table 1 . For example, for 
each 20 ms speech frame, 220, 160, 133 , 1 16 or 91 bits are produced, corresponding to bit 
rates of 1 1.0, 8.0, 6.65, 5.8 or 4.55 kbps, respectively. 

Table 1 : Bit allocation of the AMR coding algorithm for 20 ms frame 



CODING RATE 


11 .0KB PS | 8 .0KB PS | 6.65KBPS | 5.80KBPS | 4. 5 5KB PS 


Frame size 


20ms 


Look ahead 


Sms 


LPC order 


KT-order 


Predictor for LSF 
Quantization 


1 predictor 
0 bit/frame 


2 predictors: 
1 bit/frame 


LSF Quantization 


28 bit/frame 


24 bit/frame 


18 


LPC interpolation 


2 bits/frame 


2bits/f | 0 


2bits/f | 0 | 0 


0 


Coding mode bit 


Obit 


Obit 


1 bit/frame 0 bit 


Obit 


Pitch mode 


LTP 


LTP | 


LTP | PP 1 PP 


PP 


Subframe size 




5ms 




Pitch Lag 


30 bits/frame (9696) 


8585 


8585 


0008 


0008 


0008 


Fixed excitation 


31 bits/subframe 


20 


13 


18 


14 bits/subframe 


10 bits/subframe 


Gain quantization 


9 bits (scalar) 


7 bits/subframe 


6 bits/subframe 






1 1 


1 


1 


Total 


220 bits/frame 


160 | 133 


133 | 116 


91 



With reference to Fig. 5, the decoder processing circuitry, pursuant to software 
control, reconstructs the speech signal using the transmitted modeling indices extracted from 
the received bit stream by the demultiplexor 511. The decoder processing circuitry decodes 
the indices to obtain the coder parameters at each transmission frame. These parameters are 
the LSF vectors, the fractional pitch lags, the innovative code vectors, and the two gains. 
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More specifically, LP analysis at the block 239 is performed twice per frame but only 
a single set of LP parameters is converted to line spectrum frequencies (LSF) and vector 
quantized using predictive multi-stage quantization (PMVQ). The speech frame is divided 
into subframes. Parameters from the adaptive and fixed codebooks 257 and 261 are 
transmitted every subframe. The quantized and unquantized LP parameters or their 
interpolated versions are used depending on the subframe. An open-loop pitch lag is 
estimated at the block 241 once or twice per frame for PP mode or LTP mode, respectively. 

Each subframe, at least the following operations are repeated. First, the encoder 
processing circuitry (operating pursuant to software instruction) computes x(n), the first 
target signal 229, by filtering the LP residual through the weighted synthesis filter 
W( z)H( z) with the initial states of the filters having been updated by filtering the error 
between LP residual and excitation. This is equivalent to an alternate approach of subtracting 
the zero input response of the weighted synthesis filter from the weighted speech signal. 

Second, the encoder processing circuitry computes the impulse response, h(n ) , of the 
weighted synthesis filter. Third, in the LTP mode, closed-loop pitch analysis is performed to 
find the pitch lag and gain, using the first target signal 229, x(n), and impulse response, 
h(n), by searching around the open-loop pitch lag. Fractional pitch with various sample 
resolutions are used. 

In the PP mode, the input original signal has been pitch-preprocessed to match the 
interpolated pitch contour, so no closed-loop search is needed. The LTP excitation vector is 
computed using the interpolated pitch contour and the past synthesized excitation. 

Fourth, the encoder processing circuitry generates a new target signal x 2 (n), the 
second target signal 253, by removing the adaptive codebook contribution (filtered adaptive 
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where T is the pitch delay and g p is the pitch gain. 

With reference to Fig. 2, the excitation signal at the input of the short-term LP 
synthesis filter at the block 249 is constructed by adding two excitation vectors from the 
adaptive and the fixed codebooks 257 and 261, respectively. The speech is synthesized by 
feeding the two properly chosen vectors from these codebooks through the short-term 
synthesis filter at the block 249 and 267, respectively. 

The optimum excitation sequence in a codebook is chosen using an 
analysis-by-synthesis search procedure in which the error between the original and 
synthesized speech is minimized according to a perceptually weighted distortion measure. 
The perceptual weighting filter, e.g., at the blocks 251 and 268, used in the 
analysis-by-synthesis search technique is given by: 

A(z/y 2 ) K } 

where A{z) is the unquantized LP filter and 0<y 2 <Yi<larethe perceptual weighting 
factors. The values y x = [0.9, 0.94] and y 2 = 0.6 are used. The weighting filter, e.g., at the 
blocks 251 and 268, uses the unquantized LP parameters while the formant synthesis filter, 
e.g., at the blocks 249 and 267, uses the quantized LP parameters. Both the unquantized and 
quantized LP parameters are generated at the block 239. 

The present encoder embodiment operates on 20 ms (millisecond) speech frames 
corresponding to 160 samples at the sampling frequency of 8000 samples per second. At 
each 160 speech samples, the speech signal is analyzed to extract the parameters of the CELP 
model, i.e., the LP filter coefficients, adaptive and fixed codebook indices and gains. These 
parameters are encoded and transmitted. At the decoder, these parameters are decoded and 
speech is synthesized by filtering the reconstructed excitation signal through the LP synthesis 
filter. 

-17- 



WO 00/11659 



PCT/US99/19136 



3) reconstruction filter including x/sin(x) correction; and 4) an output level adjustment 
device. 

In terminal equipment, the AID function may be achieved by direct conversion to 
13-bit uniform PCM format, or by conversion to 8-bit/A-law compounded format. For the 
D/A operation, the inverse operations take place. 

The encoder 1 17 receives data samples with a resolution of 13 bits left justified in a 
16-bit word. The three least significant bits are set to zero. The decoder 133 outputs data in 
the same format. Outside the speech codec, further processing can be applied to 
accommodate traffic data having a different representation. 

A specific embodiment of an AMR (adaptive multi-rate) codec with the operational 
functionality illustrated in Figs. 2-5 uses five source codecs with bit-rates 1 1.0, 8.0, 6.65, 5.8 
and 4.55 kbps. Four of the highest source coding bit-rates are used in the full rate channel 
and the four lowest bit-rates in the half rate channel. 

All five source codecs within the AMR codec are generally based on a code-excited 
linear predictive (CELP) coding model. A 10th order linear prediction (LP), or short-term, 
synthesis filter, e.g., used at the blocks 249, 267, 301, 407 and 531 (of Figs. 2-5), is used 
which is given by: 



where a,. , i = 1, . . . , m, are the (quantized) linear prediction (LP) parameters. 

A long-term filter, i.e., the pitch synthesis filter, is implemented using the either an 
adaptive codebook approach or a pitch pre-processing approach. The pitch synthesis filter is 
given by: 




0) 



1 1 



B(z) l-g p z 



(2) 
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speech encoder, the speech decoder, which comprises decoder processing circuitry, typically 
operates pursuant to software instruction carrying out the following functionality. 

A demultiplexer 51 1 receives a bit stream 513 of speech modeling indices from an 
often remote encoder via a channel decoder. As previously discussed, the encoder selected 
each index value during the multi-stage encoding process described above in reference to 
Figs. 2-4. The decoder processing circuitry utilizes indices, for example, to select excitation 
vectors from an adaptive codebook 515 and a fixed codebook 519, set the adaptive and fixed 
codebook gains at a block 521, and set the parameters for a synthesis filter 531. 

With such parameters and vectors selected or set, the decoder processing circuitry 
generates a reproduced speech signal 539. In particular, the codebooks 515 and 519 generate 
excitation vectors identified by the indices from the demultiplexer 51 1. The decoder 
processing circuitry applies the indexed gains at the block 521 to the vectors which are 
summed. At a block 527, the decoder processing circuitry modifies the gains to emphasize 
the contribution of vector from the adaptive codebook 515. At a block 529, adaptive tilt 
compensation is applied to the combined vectors with a goal of flattening the excitation 
spectrum. The decoder processing circuitry performs synthesis filtering at the block 531 
using the flattened excitation signal. Finally, to generate the reproduced speech signal 539, 
post filtering is applied at a block 535 deemphasizing the valley areas of the reproduced 
speech signal 539 to reduce the effect of distortion. 

In the exemplary cellular telephony embodiment of the present invention, the AID 
converter 115 (Fig. la) will generally involve analog to uniform digital PCM including: 1) an 
input level adjustment device; 2) an input anti-aliasing filter; 3) a sample-hold device 
sampling at 8 kHz; and 4) analog to uniform digital conversion to 13-bit representation. 

Similarly, the D/A converter 135 will generally involve uniform digital PCM to 
analog including: 1) conversion from 13-bit/8 kHz uniform PCM to analog; 2) a hold device; 
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The speech encoding circuitry searches for optimum gain values for the previously 
identified excitation vectors ( in the first stage) from both the adaptive and fixed codebooks 
257 and 261. As indicated by blocks 307 and 309, the speech encoding circuitry identifies 
the optimum gain by generating a synthesized and weighted signal, i.e., via a block 301 and 
303, that best matches the first target signal 229 (which minimizes the third error signal 311). 
Of course if processing capabilities permit, the first and second stages could be combined 
wherein joint optimization of both gain and adaptive and fixed codebook rector selection 
could be used. 

Fig. 4 is a functional block diagram depicting of a third stage of operations performed 
by the embodiment of the speech encoder illustrated in Figs. 2 and 3. The encoder processing 
circuitry applies gain normalization, smoothing and quantization, as represented by blocks 
401, 403 and 405, respectively, to the jointly optimized gains identified in the second stage of 
encoder processing. Again, the adaptive and fixed codebook vectors used are those identified 
in the first stage processing. 

With normalization, smoothing and quantization functionally applied, the encoder 
processing circuitry has completed the modeling process. Therefore, the modeling 
parameters identified are communicated to the decoder. In particular, the encoder processing 
circuitry delivers an index to the selected adaptive codebook vector to the channel encoder 
via a multiplexor 419. Similarly, the encoder processing circuitry delivers the index to the 
selected fixed codebook vector, resultant gains, synthesis filter parameters, etc., to the 
muliplexor 419. The multiplexor 419 generates a bit stream 421 of such information for 
delivery to the channel encoder for communication to the channel and speech decoder of 
receiving device. 

Fig. 5 is a block diagram of an embodiment illustrating functionality of speech 
decoder having corresponding functionality to that illustrated in Figs. 2-4. As with the 
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At blocks 247, 249 and 251, the resultant excitation vector is applied after adaptive 
gain reduction to both a synthesis and a weighting filter to generate a modeled signal that best 
matches the first target signal 229. The encoder processing circuitry uses LPC (linear 
predictive coding) analysis, as indicated by a block 239, to generate filter parameters for the 
synthesis and weighting filters. The weighting filters 2 1 9 and 25 1 are equivalent in 
functionality. 

Next, the encoder processing circuitry designates the first error signal 253 as a second 
target signal for matching using contributions from a fixed codebook 261 . The encoder 
processing circuitry searches through at least one of the plurality of subcodebooks within the 
fixed codebook 261 in an attempt to select a most appropriate contribution while generally 
attempting to match the second target signal. 

More specifically, the encoder processing circuitry selects an excitation vector, its 
corresponding subcodebook and gain based on a variety of factors. For example, the 
encoding bit rate, the degree of minimization, and characteristics of the speech itself as 
represented by a block 279 are considered by the encoder processing circuitry at control 
block 275. Although many other factors may be considered, exemplary characteristics 
include speech classification, noise level, sharpness, periodicity, etc. Thus, by considering 
other such factors, a first subcodebook with its best excitation vector may be selected rather 
than a second subcodebook's best excitation vector even though the second subcodebook's 
better minimizes the second target signal 265. 

Fig. 3 is a functional block diagram depicting of a second stage of operations 
performed by the embodiment of the speech encoder illustrated in Fig. 2. In the second stage, 
the speech encoding circuitry simultaneously uses both the adaptive the fixed codebook 
vectors found in the first stage of operations to minimize a third error signal 311. 
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Figs. 2-4 are functional block diagrams illustrating a multi-step encoding approach 
used by one embodiment of the speech encoder illustrated in Figs, la and lb. In particular, 
Fig. 2 is a functional block diagram illustrating of a first stage of operations performed by 
one embodiment of the speech encoder shown in Figs, la and lb. The speech encoder, which 
comprises encoder processing circuitry, typically operates pursuant to software instruction 
carrying out the following functionality. 

At a block 215, source encoder processing circuitry performs high pass filtering of a 
speech signal 211. The filter uses a cutoff frequency of around 80 Hz to remove, for 
example, 60 Hz power line noise and other lower frequency signals. After such filtering, the 
source encoder processing circuitry applies a perceptual weighting filter as represented by a 
block 219. The perceptual weighting filter operates to emphasize the valley areas of the 
filtered speech signal. 

If the encoder processing circuitry selects operation in a pitch preprocessing (PP) 
mode as indicated at a control block 245, a pitch preprocessing operation is performed on the 
weighted speech signal at a block 225. The pitch preprocessing operation involves warping 
the weighted speech signal to match interpolated pitch values that will be generated by the 
decoder processing circuitry. When pitch preprocessing is applied, the warped speech signal 
is designated a first target signal 229. If pitch preprocessing is not selected the control block 
245, the weighted speech signal passes through the block 225 without pitch preprocessing 
and is designated the first target signal 229. 

As represented by a block 255, the encoder processing circuitry applies a process 
wherein a contribution from an adaptive codebook 257 is selected along with a corresponding 
gain 257 which minimize a first error signal 253. The first error signal 253 comprises the 
difference between the first target signal 229 and a weighted, synthesized contribution from 
the adaptive codebook 257. 
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The encoding system 159 comprises both a speech processing circuit 185 that 
performs speech encoding, and a channel processing circuit 187 that performs channel 
encoding. Similarly, the decoding system 165 comprises a speech processing circuit 189 that 
performs speech decoding, and a channel processing circuit 191 that performs channel 
decoding. 

Although the speech processing circuit 185 and the channel processing circuit 187 are 
separately illustrated, they might be combined in part or in total into a single unit. For 
example, the speech processing circuit 185 and the channel processing circuitry 187 might 
share a single DSP (digital signal processor) and/or other processing circuitry. Similarly, the 
speech processing circuit 189 and the channel processing circuit 191 might be entirely 
separate or combined in part or in whole. Moreover, combinations in whole or in part might 
be applied to the speech processing circuits 185 and 189, the channel processing circuits 187 
and 191, the processing circuits 185, 187, 189 and 191, or otherwise. 

The encoding system 159 and the decoding system 165 both utilize a memory 161. 
The speech processing circuit 185 utilizes a fixed codebook 181 and an adaptive codebook 
183 of a speech memory 177 in the source encoding process. The channel processing circuit 
187 utilizes a channel memory 175 to perform channel encoding. Similarly, the speech 
processing circuit 189 utilizes the fixed codebook 181 and the adaptive codebook 183 in the 
source decoding process. The channel processing circuit 187 utilizes the channel memory 
175 to perform channel decoding. 

Although the speech memory 177 is shared as illustrated, separate copies thereof can 
be assigned for the processing circuits 185 and 189. Likewise, separate channel memory can 
be allocated to both the processing circuits 187 and 191. The memory 161 also contains 
software utilized by the processing circuits 185,187,189 and 191 to perform various 
functionality required in the source and channel encoding and decoding processes. 
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With lower bit rate encoding, the speech encoder 1 17 incorporates various techniques 
to generate better low bit rate speech reproduction. Many of the techniques applied are based 
on characteristics of the speech itself. For example, with lower bit rate encoding, the speech 
encoder 117 classifies noise, unvoiced speech, and voiced speech so that an appropriate 
modeling scheme corresponding to a particular classification can be selected and 
implemented. Thus, the speech encoder 117 adaptively selects from among a plurality of 
modeling schemes those most suited for the current speech. The speech encoder 117 also 
applies various other techniques to optimize the modeling as set forth in more detail below. 

Fig. lb is a schematic block diagram illustrating several variations of an exemplary 
communication device employing the functionality of Fig. la. A communication device 151 
comprises both a speech encoder and decoder for simultaneous capture and reproduction of 
speech. Typically within a single housing, the communication device 151 might, for 
example, comprise a cellular telephone, portable telephone, computing system, etc. 
Alternatively, with some modification to include for example a memory element to store 
encoded speech information the communication device 151 might comprise an answering 
machine, a recorder, voice mail system, etc. 

A microphone 155 and an A/D converter 157 coordinate to deliver a digital voice 
signal to an encoding system 159. The encoding system 159 performs speech and channel 
encoding and delivers resultant speech information to the channel. The delivered speech 
information may be destined for another communication device ( not shown) at a remote 
location. 

As speech information is received, a decoding system 165 performs channel and 
speech decoding then coordinates with a D/A converter 167 and a speaker 169 to reproduce 
something that sounds like the originally captured speech. 
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The channel encoder 1 19 coordinates with a channel decoder 131 to deliver the 
speech indices across the communication channel 103. The channel decoder 131 forwards 
the speech indices to a speech decoder 133. While operating in a mode that corresponds to 
that of the speech encoder 117, the speech decoder 133 attempts to recreate the original 
speech from the speech indices as accurately as possible at a speaker 137 via a D/A (digital to 
analog) converter 135. 

The speech encoder 1 1 7 adaptively selects one of the plurality of operating modes 
based on the data rate restrictions through the communication channel 103. The 
communication channel 103 comprises a bandwidth allocation between the channel encoder 
1 19 and the channel decoder 131. The allocation is established, for example, by telephone 
switching networks wherein many such channels are allocated and reallocated as need arises. 
In one such embodiment, either a 22.8 kbps (kilobits per second) channel bandwidth, i.e., a 
full rate channel, or a 1 1.4 kbps channel bandwidth, i.e., a half rate channel, may be allocated. 

With the full rate channel bandwidth allocation, the speech encoder 117 may 
adaptively select an encoding mode that supports a bit rate of 1 1.0, 8.0, 6.65 or 5.8 kbps. The 
speech encoder 117 adaptively selects an either 8.0, 6.65, 5.8 or 4.5 kbps encoding bit rate 
mode when only the half rate channel has been allocated. Of course these encoding bit rates 
and the aforementioned channel allocations are only representative of the present 
embodiment. Other variations to meet the goals of alternate embodiments are contemplated. 

With either the full or half rate allocation, the speech encoder 1 1 7 attempts to 
communicate using the highest encoding bit rate mode that the allocated channel will support. 
If the allocated channel is or becomes noisy or otherwise restrictive to the highest or higher 
encoding bit rates, the speech encoder 1 17 adapts by selecting a lower bit rate encoding 
mode. Similarly, when the communication channel 103 becomes more favorable, the speech 
encoder 117 adapts by switching to a higher bit rate encoding mode. 
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Detailed Description 
Fig. la is a schematic block diagram of a speech communication system illustrating 
the use of source encoding and decoding in accordance with the present invention. Therein, a 
speech communication system 100 supports communication and reproduction of speech 
across a communication channel 103. Although it may comprise for example a wire, fiber or 
optical link, the communication channel 103 typically comprises, at least in part, a radio 
frequency link that often must support multiple, simultaneous speech exchanges requiring 
shared bandwidth resources such as may be found with cellular telephony embodiments. 

Although not shown, a storage device may be coupled to the communication channel 
103 to temporarily store speech information for delayed reproduction or playback, e.g., to 
perform answering machine functionality, voiced email, etc. Likewise, the communication 
channel 103 might be replaced by such a storage device in a single device embodiment of the 
communication system 100 that, for example, merely records and stores speech for 
subsequent playback. 

In particular, a microphone 1 1 1 produces a speech signal in real time. The 
microphone 1 1 1 delivers the speech signal to an A/D (analog to digital) converter 1 15. The 
A/D converter 1 15 converts the speech signal to a digital form then delivers the digitized 
speech signal to a speech encoder 1 1 7. 

The speech encoder 117 encodes the digitized speech by using a selected one of a 
plurality of encoding modes. Each of the plurality of encoding modes utilizes particular 
techniques that attempt to optimize quality of resultant reproduced speech. While operating 
in any of the plurality of modes, the speech encoder 1 17 produces a series of modeling and 
parameter information (hereinafter "speech indices"), and delivers the speech indices to a 
channel encoder 119. 
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BRIEF DESCRIPTION OF THF DRAWINGS 

Fig. la is a schematic block diagram of a speech communication system illustrating 
the use of source encoding and decoding in accordance with the present invention. 

Fig. lb is a schematic block diagram illustrating an exemplary communication device 
utilizing the source encoding and decoding functionality of Fig. la. 

Figs. 2-4 are functional block diagrams illustrating a multi-step encoding approach 
used by one embodiment of the speech encoder illustrated in Figs, la and lb. In particular, 
Fig. 2 is a functional block diagram illustrating of a first stage of operations performed by 
one embodiment of the speech encoder of Figs, la and lb. Fig. 3 is a functional block 
diagram of a second stage of operations, while Fig. 4 illustrates a third stage. 

Fig. 5 is a block diagram of one embodiment of the speech decoder shown in Figs, la 
and lb having corresponding functionality to that illustrated in Figs. 2-4. 

Fig. 6 is a block diagram of an alternate embodiment of a speech encoder that is built 
in accordance with the present invention. 

Fig. 7 is a block diagram of an embodiment of a speech decoder having corresponding 
functionality to that of the speech encoder of Fig. 6. 

Fig. 8 is a flow diagram illustrating the functionality of gain normalization such as 
that represented in the block 401 of Fig. 4 by an encoder built in accordance with the present 
invention. 

Fig. 9 is a flow diagram providing a more detailed description of one embodiment of 
gain normalization functionality of Fig. 8. 
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specific way to accomplish this involves the use of the linear predictive coding gain as a 
weighting factor. 

When the speech signal comprises background noise, the encoder processing system 
may also select without combination either the open loop gain or the closed loop gain in 
generating the gain normalization factor. In some cases, this involves selection of the lesser 
of the open loop gain and the closed loop gain. 

Other aspects, advantages and novel features of the present invention will become 
apparent from the following detailed description of the invention when considered in 
conjunction with the accompanying drawings. 
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SUMMARY OF THE INVENTION 
Various aspects of the present invention can be found in a speech encoding system 
using an analysis by synthesis coding approach on a speech signal. Therein, the speech 
encoding system has an encoder processing circuit and a plurality of codebooks that generate 
excitation vectors. The encoder processing circuit calculates open loop gain and closed loop 
gain. The encoder processing circuit selectively applies the open and closed loop gains in 
gain normalization processing. 

The selective application of the open loop gain and the closed loop gain by the 
encoder processing circuit may further involve the use of a weighting factor of linear 
predictive coding gain. The encoder processing circuit may use weighting factor to linearly 
combine the open and closed loop gains. 

In certain embodiments, the selective application of the open and closed loop gains by 
the encoder processing circuit comprises applying the lesser of the open loop gain and the 
closed loop gain to the background noise. When using the linear predictive coding gain as a 
weighting factor, the encoder processing circuit may exclude such background noise from the 
application of the weighting factor. 

Additionally, the encoder processing circuit may apply a maximum limit, a minimum 
limit, or both in gain normalization processing. 

Further aspects of the present invention can also be found in a method used by a 
speech encoding system that applies an analysis by synthesis coding approach to a speech 
signal. The method involves the identification of open and closed loop gains for combination 
of contributions therefrom to generate a gain normalization factor. 

In some embodiments, such method may further involve the use of linear predictive 
coding gain to identify appropriate contributions of the open and closed loop gains. One 
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Further limitations and disadvantages of conventional systems will become apparent 
to one of skill in the art after reviewing the remainder of the present application with 
reference to the drawings. 
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Background 

1. Technical Field 

The present invention relates generally to speech encoding and decoding in voice 
communication systems; and, more particularly, it relates to various techniques used with 
code-excited linear prediction coding to obtain high quality speech reproduction through a 
limited bit rate communication channel. 

2. Related Art 

Signal modeling and parameter estimation play significant roles in communicating 
voice information with limited bandwidth constraints. To model basic speech sounds, speech 
signals are sampled as a discrete waveform to be digitally processed. In one type of signal 
coding technique called LPC (linear predictive coding), the signal value at any particular time 
index is modeled as a linear function of previous values. A subsequent signal is thus linearly 
predictable according to an earlier value. As a result, efficient signal representations can be 
determined by estimating and applying certain prediction parameters to represent the signal. 

Applying LPC techniques, a conventional source encoder operates on speech signals 
to extract modeling and parameter information for communication to a conventional source 
decoder via a communication channel. Once received, the decoder attempts to reconstruct a 
counterpart signal for playback that sounds to a human ear like the original speech. 

A certain amount of communication channel bandwidth is required to communicate 
the modeling and parameter information to the decoder. In embodiments, for example where 
the channel bandwidth is shared and real-time reconstruction is necessary, a reduction in the 
required bandwidth proves beneficial. However, using conventional modeling techniques, 
the quality requirements in the reproduced speech limit the reduction of such bandwidth 
below certain levels. 
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1 1 . The method of claim 1 0 further comprising: 
utilizing linear predictive coding gain to identify appropriate contributions of the open 
loop gain and the closed loop gain. 



12. The method of claim 10 further comprising : 

utilizing linear predictive coding gain as a weighting factor to identify the 
contributions of the open loop gain and the closed loop gain. 



13. The method of any of claims 10, 1 1 and 12 further comprising: 
selecting without combination either the open loop gain or the closed loop gain in 

generating the gain normalization factor when the speech signal comprises background noise. 

14. The method of any of claims 10, 1 1 and 12 further comprising: 
choosing the lesser of the open loop gain and the closed loop gain without 

combination in generating the gain normalization factor when the speech signal comprises 
background noise. 
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6. The speech encoding system of claim 2 wherein the smooth combination of 
the open loop gain and the closed loop gain by the encoder comprises utilizing linear 
predictive coding gain. 

7. The speech encoding system of claim 2 wherein the smooth combination of 
the open loop gain and the closed loop gain by the encoder comprises a linear combination of 
the open loop gain and the closed loop gain using linear predictive coding gain as a weighting 
factor. 

8. The speech encoding system of any of claims 2 and 7 wherein the encoder 
excludes background noise from the smooth combination of the open loop gain and the 
closed loop gain. 

9. The speech encoding system of claim 8 wherein the encoder applies the lesser 
of the open loop gain and the closed loop gain to the background noise. 

10. A method used by a speech encoding system that applies an analysis by 
synthesis coding approach to a speech signal, the method comprising: 

identifying open loop gain; 
identifying closed loop gain; and 

combining contributions of the open loop gain and the closed loop gain to generate a 
gain normalization factor. 
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CLAIMS 

We claim: 

1. A speech encoding system using an analysis by synthesis approach on a 
speech signal, the speech encoding system comprising: 

an encoder that calculates open loop gain and closed loop gain; 
a plurality of codebooks that generate excitation vectors; and 

the encoder selectively applying the open loop gain and the closed loop gain in gain 
normalization processing. 

2. The speech encoding system of claim 1 wherein the encoder smoothly 
combines the open loop gain and the closed loop gain. 

3. The speech encoding system of any of claims 1 and 2 wherein the plurality of 
codebooks comprise an adaptive codebook and a fixed codebook. 

4. The speech encoding system of claim 3 wherein the encoder generates an 
adaptive codebook gain and a fixed codebook gain, and the encoder applies the gain 
normalization factor to both the adaptive codebook gain and the fixed codebook gain. 

5. The speech encoding system of claim 4 wherein the smooth combination of 
the open loop gain and the closed loop gain by the encoder comprises utilizing linear 
predictive coding gain. 
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exc3-7 exc3(ltp) 


106 


exc3-8 


107 


cxc3-9 


108 


cxc3-10 


109 


exc3-l 1 


121 


cxc4-5 exc4<ltp) 


122 


exc4-6 exc4<ltp) 


123 


cxc4-7 exc4<ltp) 


124 


exc4-8 


125 


exc4-9 


126 


exc4~10 


127 


CXC4-I 1 


30 


gain 1-5 




gain 1-6 


37 


gam2-5 


35 


gain2-6 


44 


Rain3-5 




ftain3-6 


51 


ftun4-5 


52 




60 


pncn-o 


61 


pttcn- / 


74 


cxc 1-12 


75 


cxc 1-13 


76 


excl-14 


77 


cxcl-15 


92 


exc2-12 


93 


cxc2-l3 


94 


cxc2-14 


95 


exc2-I5 


1 10 


exc3-12 


111 


exc3-13 


112 


exc3-14 


113 


cxc3-15 


128 


exc4-12 


129 


exc4-13 


130 


exc4-l4 


131 


exc4-I5 


78 


cxc 1-16 


96 


exc2-16 


114 


exc3-16 
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bus. see table aaa 


Description 


53 


mode-O 


54 


pitch-0 


55 


pitch* I 


56 


pitch -2 


57 


pitch-3 


58 


piich-4 


59 


pitcb-5 


1 


lsfl-0 


2 


IsfM 


3 


lsfl-2 


4 


lsfl-3 


5 


lsfl-4 


6 


lsfl-5 


7 


un-o 


8 


lsf2-l 


9 


un-2 


to 


bf2-3 


11 


isn-4 


12 


UI2-5 


25 


gainl-0 


26 


gainl-1 


27 


gainl-2 


28 


gain 1-3 


32 


gain2-0 


33 


gain2-l 


34 


gam2-2 


35 


Katn2-3 


39 


gain3-0 


40 


gain3-l 


41 


gain3-2 


42 


Kain3-3 


46 


gain4-0 


47 


gain4-l 


48 


pain4-2 


49 


?ain4-3 


29 


gainl-4 


36 


gatn2-4 


43 


Rain3-4 


50 


gain4-4 


62 


excl-0 piich-0(Third subframe) 


63 


exc 1 - 1 pitch- 1 (Third subframe) 


64 


exc 1-2 piich-2(Third subframe) 


65 


exc 1 -3 pitch-3{Third subframe) 


Rn 
ou 


exc2-0 pitch-5(Third subframe) 


98 


exc3-0 pitcb-0( Second subframe) 


99 


exc3-l pitch- 1 (Second subframe) 


100 


exc3~2 pitch-2( Second subframe) 


116 


exc4-0 pitch-OfFourth subframe) 


117 


exc4-l pitch- KFourth subframe) 


118 


exc4-2 pitch-2(Fouith subframe) 


13 


lsf3-0 


14 


IsO-l 


15 


lsO-2 


16 


lsf3-3 


17 


lsf3-4 


18 


kO-5 


19 


lsf4-0 


20 


lsf4-l 


21 


lsf4-2 


22 


lsf4-3 


23 


lsf4-4 


24 


lsf4-5 


81 


exc2-l exc2(ltp) 



. HRTCH) 
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129 


cxc3-8 


130 


cxc3-9 


131 


exc3-IO 


132 


exc3-!l 


133 


exc3-12 


134 


exc3-!3 


135 


exc3-I4 


136 


exc3-15 


137 


cxc3-16 


138 


exc3-17 


139 


exc3-18 


140 


exc3-19 


141 


cxc4-0 


142 


CXC4-1 


143 


cxc4-2 


144 


exc4-3 


145 


exc4-4 


146 


cxc4-5 


147 


exc4-6 


148 


CXC4-7 


149 


cxc4-8 


150 


cxc4-9 


151 


exc4-10 


152 


cxc4-M 


153 


cxc4-12 


154 


exc4-13 


155 


exc4-l4 


156 


exc4-15 


157 


exc4-16 


158 


exc4-17 


159 


exc4-18 


160 


cxc4-19 
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24 


lsf4-5 


30 


gain 1-5 


31 


gain 1-6 


37 


gain2-5 


38 


gain2-6 


44 


Kain3-5 


45 


gain3-6 


51 


ftain4-5 


52 


gasn4-6 


59 


pitch 1-6 


67 


pitch3-6 


72 


pitch2-3 


77 


pitct>4-3 


60 


pitch 1-7 


68 


pitch3-7 


73 


pitch2«4 


78 


pitch4~4 


81 


excl-0 


82 


excl-l 


83 


excl-2 


84 


excl-3 


85 


excl-4 


86 


cxcl-5 


87 


excl-6 


88 


excl-7 


89 


excl-8 


90 


excl-9 


91 


excl-10 


92 


excUM 


93 


excl-12 


94 


excl-13 


95 


excl-U 


96 


excl-15 


97 


excl-l 6 


98 


cxcl-17 


99 


excl-18 


100 


excl-19 


101 


exc2-0 


102 


exc2-I 


103 


cxc2-2 


104 


cxc2-3 


105 


exc2-4 


106 


exc2-5 


107 


cxc2-6 


108 


cxc2-7 


109 


cxc2-8 


1 10 


cxc2-9 


111 


exc2-IO 


112 


exc2-U 


113 


cxc2-12 


114 


cxc2-13 


115 


cxc2-14 


116 


exc2-15 


117 


cxc2-16 


118 


exc2-17 


119 


exc2-18 


120 


cxc2-l9 


121 


exc3-0 


122 


exc3-l 


123 


cxc3-2 


124 


exc3-3 


125 


exc3-4 


126 


exc3-5 


127 


exc3-6 


128 


cxc3-7 
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Ordering of bits according to subjective 



Bits, see table XXX 


Description 


1 


lsfl-0 


2 


bfl-l 


3 


lsfl-2 


4 


lsfl-3 


5 


lsfl-4 


6 


lsfl-5 


25 


pain 1-0 


26 


gainl-l 


27 


gain 1-2 


28 


gain I -3 


32 


gain2-0 


33 


gain2-l 


34 


gain2-2 


35 


gain2-3 


39 


gain3-0 


40 


gain3-l 


41 


gain3-2 


42 


gain3-3 


46 


gain4-0 


47 


gain4-l 


48 


gain4-2 


49 


gain4-3 


53 


pitch 1-0 


54 


pitch 1*1 


55 


pitch 1-2 


56 


pitch 1-3 


57 


pitch 1-4 


58 


pitch 1-5 


61 


pttch3-0 


62 


Pitch3-1 


63 


pitch3-2 


64 


pitch3-3 


65 


pnch3-4 


66 


pitch3-5 


69 


pitch2-0 


70 


pitch2-l 


71 


pitch2-2 


74 


pitch4-0 


75 


pitch4-l 


76 


pttch4-2 


7 


lsf2-0 


8 


UI2-I 


9 


lsf2-2 


in 




11 


Uf7-4 


12 


lsf2-5 


29 


gain 1-4 


36 


gain2*4 


43 


gain3-4 


50 


gain4-4 


79 


tnterp-0 


80 


twerp- 1 


13 


hf3-0 


14 


IsO-l 


15 


lsf3-2 


16 


isfS-3 


17 


IsfM 


18 


lsO-5 


19 


lsf4-0 


20 


lsf4-l 


21 


lsf4-2 


22 


lsf4-3 


23 


Is f 4-4 



importance (8.0 kbit/s HRTCH). 
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J 1 


£31 n 1 *6 


38 


£9iin2~6 




£&in3~6 






Dl 


cac i -\t 


/ J 


CXCxHJ 




cxc3*0 




cxc4-0 




CXC J -I 


63 


cxc 1 -2 


64 


exc 1 -3 


65 


cxc 1-4 


66 


cxc 1-5 


67 


cxc 1-6 


68 


cxc 1-7 


69 


cxc 1 -8 


70 


excl-9 


71 


cxc 1-10 


72 


excl-1 1 


73 


excl-12 


74 


excl-13 


76 


cxc2-l 


77 


cxc2-2 


78 


exc2-3 


79 


exc2-4 


80 


exc2-5 


81 


exc2-6 


82 


cxc2-7 


83 


exc2-8 


84 


cxc2-9 


85 


cxc2-10 


86 


exc2-l 1 


87 


exc2-12 


88 


exc2-13 


90 


exc3-l 


91 


exc3-2 


92 


exc3-3 


93 


exc3-4 


94 


exc3-5 


95 


exc3-6 


96 


cxc3-7 




cxc 3-8 


98 


cxc3-9 


99 


cxc3-10 


100 


exc3-11 


101 


cxc3-12 


102 


cxc3-13 


104 


exc4-l 


105 


exc4-2 


106 


exc4-3 


107 


CXC4-4 


108 


cxc4-5 


109 


exc4-6 


110 


exc4-7 


111 


cxc4-8 


112 


cxc 4-9 


113 


cxc4-10 


114 


CXC4-U 


115 


cxc4-12 


116 


exc4-13 
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132 


exc4-l6 


79 


excl-17 


97 


CXC2-I7 


I IS 


exc3-l7 


133 


exc4-17 



Ordering of bits according to subjective importance (5.8 kbit/s FRTCH). 



Bits, see table XXX 


Description 


53 


pitch-0 


54 


pitch* 1 


55 


pitch- 2 


56 


pitch'3 


57 




58 


pitch-5 


] 


151 1 


2 


left 1 
151 1-1 


3 


1ST J 


4 


ic-n i 
ISIl'3 


< 

j 


(SI I-^ 




151 1 -J 


7 


151.* -v 


8 


IsHM 


9 


lsf2-2 


10 


tsf2-3 


1 1 


lsf2-4 


12 


lsl2-5 


25 




26 




27 




28 




29 




32 


sain2*0 


33 


nin2-l 


34 


gftin?-2 


35 


R8tn2-3 


36 


eain2-4 


39 


2&in3-0 


40 


*ain3-l 


41 


pain3-2 


42 


gain3-3 


43 


Rain3-4 


46 


gain4-0 


47 


gain4-l 


48 


jtain4-2 


49 


gain4-3 


50 


jrain4-4 


30 


sain 1-5 


37 


gain2-5 


44 


gain3-5 


51 


K«n4-5 


13 


lsD-0 


14 


lsf3-l 


15 


lsf3-2 


16 


U/3-3 


17 


lsf3-4 


18 


lsD-5 


59 


pitch -6 


60 


pitch-7 


19 


lsM-0 


20 


lsf4-l 


21 


bf4-2 


22 


lsf4-3 


23 


bf4-4 


24 


lsf4-5 
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2i 


lsf4-2 


22 


lsf4-3 


67 


CAC 1 "J CXC 1 I Kp) 


68 


aw a 1 mm** 1 

CAC I -O CXC 1 \ lip/ 


69 


mitt* 1 ^ AV**t/tBn\ 

CAC 1 - / CaC 1 \1%P) 


70 


CAC 1 -O CaC HUP 7 


"7 1 


cxc ■ -7 cxc i(itp) 


72 


CXC 1 - 1 \i 


8 1 
o I 


excz-i cxc^ucp) 


O J 


cxc2-2 exc2(ltp) 


Ol 
83 


exc2-3 exc2(Itp) 


84 


cxc2-4 exc2(ltp) 


85 


cxc2-5 exc2(ltp) 


86 


exc 2 -6 exc2(ltp) 


87 


cxc2-7 


88 


cxc2-8 


89 


cxc2-9 


90 


cxc2-10 


101 


cxc3-3 cxc3<»tp) 


102 


exc3-4 cxc3(ltp) 


103 


cxc3-5 exc3(ltp) 


104 


exc3-6 cxc3(tq» 


105 


exc3-7 exc3(ltp) 


106 


exc3«8 


107 


exc3-9 


108 


exc3-10 


119 


exc4-3 exc4(ltp) 


120 


exc4-4 exc4<hp) 


121 


cxc4-5 cxc4<ttp) 


122 


cxc4~6 exc40tp) 


123 


exc4-7 exc4(ltp) 


124 


exc4-8 


125 


cxc4-9 


126 


exc4-10 


73 


excl-11 


91 


cxc2-l I 


109 


exc3-l 1 


127 


CXC4-1 1 


74 


cxc 1-12 


92 


cxc2-12 


1 10 


cxc3- 1 2 


i *>o 


cxc4- 1 2 


OU 


pstch-6 


Al 
Ol 


pitch -7 


11 

Zj 


Is f 4-4 


24 


ISI40 


75 


exc 1-1 3 


93 


cxc2-13 


111 


exc3-13 


129 


exc4-13 


31 


gain 1-6 


38 


gain2-6 


45 


gain3-6 


52 


gain4-6 


76 


cxcl-14 


77 


cxcl-15 


94 


cxc2-14 


95 


cxc2-l5 


112 


cxc3-14 


113 


cxc3-15 


130 


exc4-!4 


131 


exc4-15 


78 


cxcl-16 


96 


exc2-16 


114 


cxc3-16 
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Ordering of bits according to subjective importance (6.6S kbi t/s FRTCH). 



Bits, see table XXX 


Description 


54 


pitcb-0 


55 


pitch- 1 


56 


pttch-2 


57 


pitch- 3 


58 


pitch-4 


59 


pitch- 5 


1 


lsfl-0 


2 


lsfl-1 


3 


lsfl-2 


4 


lsfl-3 


5 


lsfl-4 


6 


lsfl-5 


25 


gainl -0 


26 


gain! -1 


27 


gain 1 -2 


28 




32 


pain 2-0 


33 




34 


gain2-2 


35 




39 


gain3-0 


40 




41 




42 




46 




47 




48 


gauw-z 


49 


gain4-3 


29 


saints 


36 


gain2^ 


43 


gain3-4 


50 


gain4-4 


53 


mode-0 


98 


exc3-0 nitch-(tf SeennH ctthfrnrwl 


99 


CXC3-1 nitrh- 1 (^erruwi ciihfmnr^l 


7 


1sQ-0 


8 


lsf2-l 


9 


lsD-2 


10 


ls!2-3 


11 


UI2-4 


12 


lsf2-5 


30 


gain 1-5 


37 


gain2-5 


44 


gain3-5 


51 


gain4-S 


62 


excl-0 pitch-OCThird subframe) 


63 


excl-1 pitch- 1 (Third subframe) 


64 


exc I -2 pitch-2(Third subframe) 


65 


excl-3 pitch -3(Third subframe) 


66 


exc 1-4 pitcrwKThird subframe) 


80 


exc2-0 pitch«5(Third subframe) 


100 


exc3-2 pitch-2( Second subframe) 


116 


exc4-0 oitcb-0( Fourth subframe) 


117 


exc4-| pitch- 1 (Fourth subframe) 


118 


exc4-2 pitch -2(Founh subframe) 


13 


lsf3-0 


14 


isn-i 


15 


lsf3-2 


16 


lsf3-3 


17 


lsO-4 


18 


lsf3-5 


19 


lsf4-0 


20 | kf4-l 
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128 


exc3-7 


129 


exc3-8 


130 


exc3-9 


131 


exc3-10 


132 


exc3-ll 


133 


exc3-12 


134 


exc3-13 


135 


exc3-14 


136 


exc3-!5 


137 


exc3-16 


138 


exc3-!7 


139 


cxc3-18 


140 


exc3-19 


141 


exc4-0 


142 


exc4-l 


143 


exc4-2 


144 


exc4-3 


145 


CXC4-4 


146 


cxc4-5 


147 


cxc4-6 


148 


exc4-7 


149 


exc4-8 


150 


exc4-9 


151 


cxc4-10 


152 


exc4-l 1 


153 


cxc4-12 


154 


exc4-I3 


155 


exc4-14 


156 


exc4-15 


157 


exc4-16 


158 


cxc4-17 


159 


exc4-18 


160 


cxc4-19 
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lag is less than the subframe size 40 and the chosen excitation is pulse type, the pitch 
sharpening is applied. This translates into modifying c(n) as c(n) = c(n) + fic{n-T) + 
where p is the decoded pitch gain g p from the previous subframe bounded by [0.2,1.0]. 
The excitation at the input of the synthesis filter is given by 
u ( n ) = g p v( n ) + g c c ( n )> n - °>39 . Before the speech synthesis, a post-processing of the 
excitation elements is performed. This means that the total excitation is modified by 
emphasizing the contribution of the adaptive codebook vector: 

\u(n) + 0.25 J3g p v(n), g p >0.5 
u(ri) = < p p 

1 u{n\ g p <=0.5 

Adaptive gain control (AGC) is used to compensate for the gain difference between the 
unemphasized excitation u(ri) and emphasized excitation u{n) . The gain scaling factor rj 
for the emphasized excitation is computed by: 



7 = 



39 




I> 








39 






2 («) 



f,>0.5 



B«0 

10 f p <=0.5 
The gain-scaled emphasized excitation u(ri) is given by: 

The reconstructed speech is given by: 
s(n) = u \ri) - ^ a { s(n - /)> w = 0 to 39 , 

where a i are the interpolated LP filter coefficients. The synthesized speech s(n) is then 
passed through an adaptive postfilter. 

Post-processing consists of two functions: adaptive postfiltering and signal up-scaling. 

The adaptive postfilter is the cascade of three filters: a formant postfilter and two tilt 
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compensation filters. The postfilter is updated every subframe of 5 ms. The forxnant 
postfilter is given by: 



where A(z) is the received quantized and interpolated LP inverse filter and y n and y d control 
the amount of the formant postfiltering. 

The first tilt compensation filter H n (z) compensates for the tilt in the formant 
postfilter H f (z) and is given by: 

where // = y tX k x is a tilt factor, with k x being the first reflection coefficient calculated on the 
truncated impulse response h f (n) , of the formant postfilter k x = with: 

r h (i) = A,C/)*/O' + 0 , (L h = 22) . 

The postfiltering process is performed as follows. First, the synthesized speech s(n) 
is inverse filtered through A( z / ) to produce the residual signal r(ri) . The signal r(n) is 

filtered by the synthesis filter \/A(zfy d ) is passed to the first tilt compensation filter 
h n (z) resulting in the postfiltered speech signal s f (n) . 

Adaptive gain control (AGC) is used to compensate for the gain difference between 
the synthesized speech signal s(n)and the postfiltered signal s f (w) . The gain scaling factor 
y for the present subframe is computed by: 
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r = 



f> 2 (") 



/»=0 



The gain-scaled postfiltered signal s (n) is given by: 
s(n) = j3(rt)s f (n) 

where J3(n) is updated in sample by sample basis and given by: 

where a is an AGC factor with value 0.9. Finally, up-scaling consists of multiplying the 
postfiltered speech by a factor 2 to undo the down scaling by 2 which is applied to the input 
signal. 

Figs. 6 and 7 are drawings of an alternate embodiment of a 4 kbps speech codec that 
also illustrates various aspects of the present invention. In particular, Fig. 6 is a block 
diagram of a speech encoder 601 that is built in accordance with the present invention. The 
speech encoder 601 is based on the analysis-by-synthesis principle. To achieve toll quality at 
4 kbps, the speech encoder 601 departs from the strict waveform-matching criterion of 
regular CELP coders and strives to catch the perceptual important features of the input signal. 

The speech encoder 601 operates on a frame size of 20 ms with three subframes (two 
of 6.625 ms and one of 6.75 ms). A look-ahead of 15 ms is used. The one-way coding delay 
of the codec adds up to 55 ms. 

At a block 615, the spectral envelope is represented by a 10 th order LPC analysis for 
each frame. The prediction coefficients are transformed to the Line Spectrum Frequencies 
(LSFs) for quantization. The input signal is modified to better fit the coding model without 
loss of quality. This processing is denoted "signal modification" as indicated by a block 621. 
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In order to improve the quality of the reconstructed signal, perceptual important features are 
estimated and emphasized during encoding. 

The excitation signal for an LPC synthesis filter 625 is build from the two traditional 
components: 1) the pitch contribution; and 2) the innovation contribution. The pitch 
contribution is provided through use of an adaptive codebook 627. An innovation codebook 
629 has several subcodebooks in order to provide robustness against a wide range of input 
signals. To each of the two contributions a gain is applied which, multiplied with their 
respective codebook vectors and summed, provide the excitation signal. 

The LSFs and pitch lag are coded on a frame basis, and the remaining parameters (the 
innovation codebook index, the pitch gain, and the innovation codebook gain) are coded for 
every subframe. The LSF vector is coded using predictive vector quantization. The pitch lag 
has an integer part and a fractional part constituting the pitch period. The quantized pitch 
period has a non-uniform resolution with higher density of quantized values at lower delays. 
The bit allocation for the parameters is shown in the following table. 

Table of Bit Allocation 



Parameter 


Bits per 20 ms 


LSFs 


21 


Pitch lag (adaptive codebook) 


8 


Gains 


12 


Innovation codebook 


3x13 = 39 


Total 


80 



When the quantization of all parameters for a frame is complete the indices are multiplexed 

to form the 80 bits for the serial bit-stream. 

Fig. 7 is a block diagram of a decoder 701 with corresponding functionality to that of 

the encoder of Fig. 6. The decoder 701 receives the 80 bits on a frame basis from a 

demultiplexer 71 1. Upon receipt of the bits, the decoder 701 checks the sync-word for a bad 

frame indication, and decides whether the entire 80 bits should be disregarded and frame 
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erasure concealment applied. If the frame is not declared a frame erasure, the 80 bits are 
mapped to the parameter indices of the codec, and the parameters are decoded from the 
indices using the inverse quantization schemes of the encoder of Fig. 6. 

When the LSFs, pitch lag, pitch gains, innovation vectors, and gains for the 
innovation vectors are decoded, the excitation signal is reconstructed via a block 715. The 
output signal is synthesized by passing the reconstructed excitation signal through an LPC 
synthesis filter 721 . To enhance the perceptual quality of the reconstructed signal both short- 
term and long-term post-processing are applied at a block 731. 

Regarding the bit allocation of the 4 kbps codec (as shown in the prior table), the 
LSFs and pitch lag are quantized with 21 and 8 bits per 20 ms, respectively. Although the 
three subframes are of different size the remaining bits are allocated evenly among them. 
Thus, the innovation vector is quantized with 13 bits per subframe. This adds up to a total of 
80 bits per 20 ms, equivalent to 4 kbps. 

The estimated complexity numbers for the proposed 4 kbps codec are listed in the 
following table. All numbers are under the assumption that the codec is implemented on 
commercially available 16-bit fixed point DSPs in full duplex mode. All storage numbers are 
under the assumption of 16-bit words, and the complexity estimates are based on the floating 
point C-source code of the codec. 

Table of Complexity Estimates 



Computational complexity 


30 MIPS 


Program and data ROM 


18 kwords 


RAM 


3 kwords 



The decoder 701 comprises decode processing circuitry that generally operates 
pursuant to software control. Similarly, the encoder 601 (Fig. 6) comprises encoder 
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processing circuitry also operating pursuant to software control. Such processing circuitry 
may coexists, at least in part, within a single processing unit such as a single DSP. 
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Fig. 8 is a flow diagram illustrating the functionality of gain normalization such as 
that represented in the block 401 of Fig. 4 by an encoder built in accordance with the present 
invention. In particular, at a block 807, encoding processing circuitry operating pursuant to 
software direction begins normalization processing. At a block 811, closed loop gain (CL_G) 
is calculated. Similarly, open loop gain (OL_G) is calculated at a block 815. 

For background noise, as determined at the block 819, a gain normalization factor (G) 
is not limited or smoothed. Instead, for best perceptual quality, the encoding processing 
circuitry selects the lesser of open and closed loop gains to be the normalization factor, at a 
block 823. Thereafter, the adaptive and fixed codebook gains, represented by the blocks 259 
and 263 (of Fig. 4), are normalized as follows: 

gp' = g p x G, 
gc' = gc x G, 

where the resultant g p ' and g c ' comprise the normalized gains applied to adaptive codebook 

vectors and fixed codebook vectors, respectively. 

If the speech signal does not constitute background noise, the encoder processing 
circuit performs a smoothing algorithm at a block 835. The smoothing algorithm is one in 
which contributions from the closed and open loop gains are combined so as to prevent 
discontinuity that might otherwise accompany switching between such gains. Having 
smoothed the contributions of both the closed and open loop gains into a gain normalization 
factor (G), g p * and g c * are calculated as previously with reference to a block 827. 

The smoothing could be linear or otherwise, and might be based on any of a variety of 
factors. In particular embodiments, the smoothing comprises a linear combination based on a 
weighting coefficient of linear predictive coding (LPC) gain. The lower the LPC gain, the 
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greater the open loop gain contribution to a gain normalization factor. The greater the LPC 
gain, the greater the closed loop gain contribution. 

Fig. 9 is a flow diagram providing a more detailed description of one embodiment of 
gain normalization functionality of Fig. 8. Upon beginning the gain normalization procedure, 
the encoder processing circuitry calculates the open and closed loop gains and the LPC gain 
(^LPC). a * a block 91 1 . Next the encoder processing circuitry determines whether the 

present speech signal comprising background noise as indicated at a block 915. If the signal 
comprises background noise, the gain normalization factor (G) is set to the lesser of the open 
and close loop gains at a block 919. Thereafter, at a block 923, the adaptive and fixed 
codebook gains are normalized as follows: 

gp' - gp x G, 
gc'=g c xG, 

where the resultant g p ' and g c ' comprise the normalized gains applied to the adaptive and 
fixed codebook vectors, respectively. 

If background noise is not indicated at the block 915, the encoder processing circuitry 
checks to see if the LPC gain is between zero and one at a block 927. If so, the gain 
normalization factor is set to the value of the open loop gain, at a block 931, and applied as 
described above in reference to the block 923. Otherwise, the open and closed loop gains are 
smoothly combined at a block 935. 

Because the LPC gain changes slowly, a linear combination of the open and closed 
loop gains based on LPC gain weighting yields a slowly changing gain normalization factor. 
In particular, the generation of this gain normalization factor is given by: 
G = (QClpc) CL_G + (1 - cClpc) OL_G. 
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Thus, when the LPC gain is low for example with unvoiced speech, the open loop gain 
contribution is emphasized. Likewise, when the LPC gain is high, the closed loop gain 
contribution is emphasized. 

After calculating the gain normalization factor at the block 935, the encoder 
processing circuitry limits the gain normalization factor to a predetermined range. 
Specifically, if the gain normalization factor is not greater or equal to one as tested at a block 
939, the gain normalization factor is set to one at a block 943. With the value of one, no 
adjustment is made to the adaptive and fixed codebook vector gains at the block 923. 

If the gain normalization factor is greater or equal to one but less than one plus the 
LPC gain (as indicated by the blocks 939 and 947), the gain normalization factor need not be 
limited and is therefore applied in the gain normalization process at the block 923. However, 
if the gain normalization factor is not less than one plus the LPC gain, the gain normalization 
factor is set to one plus the LPC gain and delivered to the gain normalization process at the 
block 923. Such limitation is used mainly because the large memory energy of the synthesis 
filter due to the large normalization factor could influence the coding quality of the following 
subframe when the LPC gain is high. 

By applying such smoothed gain normalization, especially at low bit rates for noisy 
speech, the perceptual quality of the synthesized speech is significantly improved. Of course 
many other ways of smoothly combining the closed and open loop gains and other limits (or 
none at all) could be alternatively or additionally applied. 

Of course, many other modifications and variations are also possible. In view of the 
above detailed description of the present invention and associated drawings, such other 
modifications and variations will now become apparent to those skilled in the art. It should 
also be apparent that such other modifications and variations may be effected without 
departing from the spirit and scope of the present invention. 
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In addition, the following Appendix A provides a list of many of the definitions, 
symbols and abbreviations used in this application. Appendices B and C respectively provide 
source and channel bit ordering information at various encoding bit rates used in one 
embodiment of the present invention. Appendices A, B and C comprise part of the detailed 
description of the present application, and, otherwise, are hereby incorporated herein by 
reference in its entirety. 
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APPENDIX A 



For purposes of this application, the following symbols, definitions and abbreviations 



apply. 



adaptive postfilter: 



adaptive codebook: The adaptive codebook contains excitation vectors that are 

adapted for every subframe. The adaptive codebook is derived 
from the long term filter state. The pitch lag value can be 
viewed as an index into the adaptive codebook. 

The adaptive postfilter is applied to the output of the short term 
synthesis filter to enhance the perceptual quality of the 
reconstructed speech. In the adaptive multi-rate codec (AMR), 
the adaptive postfilter is a cascade of two filters: a formant 
postfilter and a tilt compensation filter. 

Adaptive Multi Rate codec: The adaptive multi-rate code (AMR) is a speech and channel 

codec capable of operating at gross bit-rates of 1 1 .4 kbps 
("half-rate") and 22.8 kbs ("full-rate"). In addition, the codec 
may operate at various combinations of speech and channel 
coding (codec mode) bit-rates for each channel mode. 



AMR handover: 

channel mode: 

channel mode adaptation: 

channel repacking: 

closed-loop pitch analysis: 



codec mode: 



codec mode adaptation: 



direct form coefficients: 



Handover between the full rate and half rate channel modes to 
optimize AMR operation. 

Half-rate (HR) or full-rate (FR) operation. 

The control and selection of the (FR or HR) channel mode. 

Repacking of HR (and FR) radio channels of a given radio cell 
to achieve higher capacity within the cell. 

This is the adaptive codebook search, i.e., a process of 
estimating the pitch (lag) value from the weighted input speech 
and the long term filter state. In the closed-loop search, the lag 
is searched using error minimization loop 
(analysis-by-synthesis). In the adaptive multi rate codec, 
closed-loop pitch search is performed for every subframe. 

For a given channel mode, the bit partitioning between the 
speech and channel codecs. 

The control and selection of the codec mode bit-rates. 
Normally, implies no change to the channel mode. 

One of the formats for storing the short term filter parameters. 
In the adaptive multi rate codec, all filters used to modify 
speech samples use direct form coefficients. 
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fixed codebook: 



fractional lags: 

full-rate (FR): 
frame: 

gross bit-rate: 

half-rate (HR); 
in-band signaling: 

integer lags: 
interpolating filter: 

inverse filter: 

lag: 

Line Spectral Frequencies: 
Line Spectral Pair: 



The fixed codebook contains excitation vectors for speech 
synthesis filters. The contents of the codebook are non-adaptive 
(i.e., fixed). In the adaptive multi rate codec, the fixed 
codebook for a specific rate is implemented using a multi- 
function codebook. 

A set of lag values having sub-sample resolution. In the 
adaptive multi rate codec a sub-sample resolution between l/6 th 
and 1.0 of a sample is used. 

Full-rate channel or channel mode. 

A time interval equal to 20 ms (160 samples at an 8 kHz 
sampling rate). 

The bit-rate of the channel mode selected (22.8 kbps or 1 1.4 
kbps). 

Half-rate channel or channel mode. 

Signaling for DTX, Link Control, Channel and codec mode 
modification, etc. carried within the traffic. 

A set of lag values having whole sample resolution. 

An FIR filter used to produce an estimate of sub-sample 
resolution samples, given an input sampled with integer sample 
resolution. 

This filter removes the short term correlation from the speech 
signal. The filter models an inverse frequency response of the 
vocal tract. 

The long term filter delay. This is typically the true pitch 
period, or its multiple or sub-multiple. 

(see Line Spectral Pair) 

Transformation of LPC parameters. Line Spectral Pairs are 
obtained by decomposing the inverse filter transfer function 
A(z) to a set of two transfer functions, one having even 
symmetry and the other having odd symmetry. The Line 
Spectral Pairs (also called as Line Spectral Frequencies) are the 
roots of these polynomials on the z-unit circle). 
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LP analysis window: 



LP coefficients: 

LTP Mode: 
mode: 

multi-function codebook: 



open-loop pitch search: 



out-of-band signaling: 
PP Mode: 
residual: 

short term synthesis filter: 
perceptual weighting filter: 



For each frame, the short term filter coefficients are computed 
using the high pass filtered speech samples within the analysis 
window. In the adaptive multi rate codec, the length of the 
analysis window is always 240 samples. For each frame, two 
asymmetric windows are used to generate two sets of LP 
coefficient coefficients which are interpolated in the LSF 
domain to construct the perceptual weighting filter. Only a 
single set of LP coefficients per frame is quantized and 
transmitted to the decoder to obtain the synthesis filter. A 
lookahead of 25 samples is used for both HR and FR 

Linear Prediction (LP) coefficients (also referred as Linear 
Predictive Coding (LPC) coefficients) is a generic descriptive 
term for describing the short term filter coefficients. 

Codec works with traditional LTP. 

When used alone, refers to the source codec mode, i.e., to one 
of the source codecs employed in the AMR codec. (See also 
codec mode and channel mode.) 

A fixed codebook consisting of several subcodebooks 
constructed with different kinds of pulse innovation vector 
structures and noise innovation vectors, where codeword from 
the codebook is used to synthesize the excitation vectors. 

A process of estimating the near optimal pitch lag directly from 
the weighted input speech. This is done to simplify the pitch 
analysis and confine the closed-loop pitch search to a small 
number of lags around the open-loop estimated lags. In the 
adaptive multi rate codec, open-loop pitch search is performed 
once per frame for PP mode and twice per frame for LTP 
mode. 

Signaling on the GSM control channels to support link control. 
Codec works with pitch preprocessing. 

The output signal resulting from an inverse filtering operation. 

This filter introduces, into the excitation signal, short term 
correlation which models the impulse response of the vocal 
tract. 

This filter is employed in the analysis-by-synthesis search of 
the codebooks. The filter exploits the noise masking properties 
of the formants (vocal tract resonances) by weighting the error 
less in regions near the formant frequencies and more in 
regions away from them. 



-79- 



f ^ *-i ni™^ ~rrF*«-F* 



WO 00/11659 



PCT/US99/I9136 



subframe: 

vector quantization: 
zero input response: 
zero state response: 

A(z) 
A(z) 



H(z) = 



1 

B(z) 
W{z) 

r t >r 2 

F E {z) 
T 



f A{z/y d ) 



H t (z) 



A time interval equal to 5-10 ms (40-80 samples at an 8 kHz 
sampling rate). 

A method of grouping several parameters into a vector and 
quantizing them simultaneously. 

The output of a filter due to past inputs, i.e. due to the present 
state of the filter, given that an input of zeros is applied. 

The output of a filter due to the present input, given that no past 
inputs have been applied, i.e., given the state information in the 
filter is all zeroes. 

The inverse filter with unquantized coefficients 
The inverse filter with quantized coefficients 

The speech synthesis filter with quantized coefficients 

The unquantized linear prediction parameters (direct form 
coefficients) 

The quantized linear prediction parameters 
The long-term synthesis filter 

The perceptual weighting filter (unquantized coefficients) 
The perceptual weighting factors 
Adaptive pre-filter 

The nearest integer pitch lag to the closed-loop fractional pitch 
lag of the subframe 

The adaptive pre-filter coefficient (the quantized pitch gain) 
The formant postfilter 

Control coefficient for the amount of the formant post-filtering 
Control coefficient for the amount of the formant post-filtering 
Tilt compensation filter 
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Y/ Control coefficient for the amount of the tilt compensation 

filtering 

f* = Yt k \ A tilt factor, with k{ being the first reflection coefficient 

hf( n ) The truncated impulse response of the formant postfilter 

Lh The length of hf(n) 

r h(0 The auto-correlations of hf(n) 

M z hn) The inverse filter (numerator) part of the formant postfilter 

\/A(z/y d ) The synthesis filter (denominator) part of the formant postfilter 

rfa) The residual signal of the inverse filter A(z /y n ) 

Impulse response of the tilt compensation filter 

Psc( n ) The AGC-controlled gain scaling factor of the adaptive 

postfilter 

<x The AGC factor of the adaptive postfilter 

H hx {z) Pre-processing high-pass filter 

w 7 (/2) , w TI (ri) LP analysis windows 

Length of the first part of the LP analysis window w /( n ) 
^ (/) Length of the second part of the LP analysis window w i( n ) 

t(JI) 

1 Length of the first part of the LP analysis window w in n ) 

^ </7) Length of the second part of the LP analysis window 

r ac W The auto-correlations of the windowed speech s f (n) 

w tag(0 Lag window for the auto-correlations (60 Hz bandwidth 

expansion) 

f° The bandwidth expansion in Hz 

f s The sampling frequency in Hz 
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ac The modified (bandwidth expanded) auto-correlations 

e ld(*) The prediction error in the zth iteration of the Levinson 

algorithm 

k i The zth reflection coefficient 

^ The yth direct form coefficient in the ith iteration of the 

Levinson algorithm 

F \ ( z ) Symmetric LSF polynomial 

Fl ^ Antisymmetric LSF polynomial 

F \ ( z ) Polynomial F x ' (z) with root z = -1 eliminated 

F 2 ( 2 ) Polynomial F i 00 with root z = 1 eliminated 

9i The line spectral pairs (LSFs) in the cosine domain 

q An LSF vector in the cosine domain 

<l/ n) The quantized LSF vector at the xth subframe of the frame n 

mi The line spectral frequencies (LSFs) 

T m (*) A m th order Chebyshev polynomial 

f\ (0» fi (0 The coefficients of the polynomials F x (z) and F 2 (z) 

f\ (0> h (0 The coefficients of the polynomials F x ' (z) and F 2 ' (z) 

f(0 The coefficients of either F } (z) or F 2 (z ) 

C( x ) Sum polynomial of the Chebyshev polynomials 

■* Cosine of angular frequency co 

*>k Recursion coefficients for the Chebyshev polynomial 

evaluation 

fi The line spectral frequencies (LSFs) in Hz 

=[fl fl • • • -/io] The vector representation of the LSFs in Hz 
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z (lJ (n) 9 z w (n) The mean-removed LSF vectors at frame n 

n ) , r^ 2 V " ) The LSF prediction residual vectors at frame n 

p(/i) The predicted LSF vector at frame n 

n-l) The quantized second residual vector at the past frame 

f * The quantized LSF vector at quantization index k 

The LSF quantization error 
w. , i = 1, . . . ,1 0 f LSF-quantization weighting factors 

d i The distance between the line spectral frequencies f M and 

fi-X 

h(n) The impulse response of the weighted synthesis filter 

O k The correlation maximum of open-loop pitch analysis at delay 

k 

O tj , f=l, ... ,3 The correlation maxima at delays t i ? , / = 1, . . . ,3 

(A/,,*,), i = 1,..., 3 The normalized correlation maxima M i and the corresponding 

delays t h i = 1,...,3 

H(z)W(z) = . A ( z 'y^ The weighted synthesis filter 
A(z)A(z/r 2 ) 

A(z/y\) The numerator of the perceptual weighting filter 

l/A(z/y 2 ) The denominator of the perceptual weighting filter 

7} The nearest integer to the fractional pitch lag of the previous 

(1st or 3rd) subframe 

s f (n) The windowed speech signal 

s w ( n ) The weighted speech signal 

s(n) Reconstructed speech signal 

s '(n) The gain-scaled post-filtered signal 

Sf(n) Post-filtered speech signal (before scaling) 
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x ( n ) The target signal for adaptive codebook search 

x 2 (n) ^ %2 The target signal for Fixed codebook search 

res^n) The LP residual signal 

c{n) The fixed codebook vector 

v ( n) The adaptive codebook vector 

y(n) = v(w)* h{n) The filtered adaptive codebook vector 

The filtered fixed codebook vector 

y k ( n ) The past filtered excitation 

u(n) The excitation signal 

£( n ) The fully quantized excitation signal 

" f ( w ) The gain-scaled emphasized excitation signal 

T op The best open-loop lag 

tmin Minimum lag search value 

t max Maximum lag search value 

Correlation term to be maximized in the adaptive codebook 
search 

R {k) t The interpolated value of R(k) for the integer delay k and 

fraction/ 

Correlation term to be maximized in the algebraic codebook 
search at index k 

Q The correlation in the numerator of A k at index k 

^Dk The energy in the denominator of A k at index A: 

d = H' x 2 The correlation between the target signal x z {ti) and the impulse 

response h(n) 9 i.e., backward filtered target 

H The lower triangular Toepliz convolution matrix with diagonal 

h(0) and lower diagonals /r(l), . . . , h(39) 

-84- 



WO 00/11659 



PCT/US99/I9136 



^ = H*H The matrix of correlations of h(n) 

d («) The elements of the vector d 

J) The elements of the symmetric matrix <X> 

c i The innovation vector 

C The correlation in the numerator of A k 

n*i The position of the i th pulse 

«9| The amplitude of the i th pulse 

N p The number of pulses in the fixed codebook excitation 

E D The energy in the denominator of A k 

res LTP {n) The normalized long-term prediction residual 

b ( n ) The sum of the normalized d(n) vector and normalized 

long-term prediction residual res LTP (n) 

s bi n ) The sign signal for the algebraic codebook search 

> The fixed codebook vector convolved with h(n ) 

E( n ) The mean-removed innovation energy (in dB) 

E The mean of the innovation energy 

E( n) The predicted energy 

[^1 ^2 h ^4] The MA prediction coefficients 

R(k) The quantized prediction error at subframe k 

E/ The mean innovation energy 

R( n ) The prediction error of the fixed-codebook gain quantization 

Eq The quantization error of the fixed-codebook gain quantization 

e(n) The states of the synthesis filter \/A(z) 
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e w( n ) The perceptually weighted error of the analysis-by-synthesis 

search 

V The gain scaling factor for the emphasized excitation 

g c The fixed-codebook gain 

g c The predicted fixed-codebook gain 

g c The quantized fixed codebook gain 

g p The adaptive codebook gain 

g p The quantized adaptive codebook gain 

« 

y gc-Sc 1 ' S c A correction factor between the gain g c and the estimated one 

Sc 

Y gc The optimum value for y gc 

y sc Gain scaling factor 

AGC Adaptive Gain Control 

AMR Adaptive Multi Rate 

CELP Code Excited Linear Prediction 

C/l Carrier-to-Interferer ratio 

DTX Discontinuous Transmission 

EFR Enhanced Full Rate 

FIR Finite Impulse Response 

FR Full Rate 

HR HalfRate 

LP Linear Prediction 

LPC Linear Predictive Coding 

LSF Line Spectral Frequency 

LSF Line Spectral Pair 
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LT P Long Term Predictor (or Long Term Prediction) 

MA Moving Average 

TFO Tandem Free Operation 

VAD Voice Activity Detection 
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APPENDIX B 

Bit ordering (source coding) 



Bit ordering of output bits from source encoder (1 1 kbit/s). 



Bits 


Description 


1-6 


Index of 1"LSF stage 


7-12 


Index of 2 nd LSF stage 


13-18 


Index of 3'° LSF stage 


19-24 


Index of 4" LSF stage 


25-28 


Index of 5™ LSF stage 


29-32 


Index of adaptive codebook gain, 1 u subframe 


33-37 


Index of fixed codebook gain, 1 " subframe 


38-41 


Index of adaptive codebook gain, 2 W subframe 


42^6 


Index of fixed codebook gam, 2°° subframe 


47-50 


Index of adaptive codebook gain, 3 ,B subframe 


51-55 


Index of fixed codebook gain, 3™ subframe 


56-59 


Index of adaptive codebook gain, 4™ subframe 


60-64 


Index of fixed codebook gain, 4 U subframe 


65-73 


Index of adaptive codebook, 1" subframe 


74-82 


Index of adaptive codebook, 3 ,a subframe 


83-88 


Index of adaptive codebook (relative), 2 M subframe 


89-94 


Index of adaptive codebook (relative), 4" subframe 


95-94 


Index for LSF interpolation 


97-127 


Index for fixed codebook, 1" subframe 


128-158 


Index for fixed codebook, 2"° subframe 


159-189 


Index for fixed codebook, 3'° subframe 


190-220 


Index for fixed codebook, 4 a subframe 



Bit ordering of output bits from source encoder (8 kbit/s). 



Bits 


Description 


1-6 


Index of 1" LSF stage 


7-12 


Index of2" u LSF stage 


13-18 


Index of 3^ LSF stage 


19-24 


Index of 4" LSF stage 


25-31 


Index of fixed and adaptive codebook gains, 1 M subframe 


32-38 


Index of fixed and adaptive codebook gains, 2 00 subframe 


39-45 


Index of fixed and adaptive codebook gains, 3'° subframe 


46-52 


Index of fixed and adaptive codebook gains, 4" subframe 


53-60 


Index of adaptive codebook, 1" subframe 


61-68 


Index of adaptive codebook, 3 re subframe 


69-73 


Index of adaptive codebook (relative), 2™ 1 subframe 


74-78 


Index of adaptive codebook (relative), 4 U subframe 


79-80 


Index for LSF interpolation 


81-100 


Index for fixed codebook, I s * subframe 


101-120 


Index for fixed codebook, 2°° subframe 


121-140 


Index for fixed codebook, 3™ subframe 


141-160 


Index for fixed codebook, 4 m subframe 
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Bit ordering of output bits from source encoder (6.65 kbit/s). 



Bits 




1-6 


Index of 1"LSF stage " 


7-12 


Index of 2°° LSF stage ~ 


13-18 


Index of3' B LSF stage ' 


19-24 


Index of 4*" LSF stage — 


25-31 


Index of fixed and adaptive codebook gains, 1* subframe 


32-38 


Index of fixed and adaptive codebook gains, 2°° subframe 


39-45 


Index of fixed and adaptive codebook gains, 3 ,Q subframe 


46-52 


Index of fixed and adaptive codebook gains, 4" subframe ~" 1 — " — ~ 


53 


Index for mode (LTP or PP) " — " 


LTPmode 




PP mode 


54-61 


Index of adaptive codebook, 1" subframe 




Index of pitch 


62-69 


Index of adaptive codebook, 3 W subframe 






70-74 


Index of adaptive codebook (relative). 2°° subframe 






75-79 


Index of adaptive codebook (relative), 4 U subframe 






80-81 


index for LSF interpolation 




Index for LSF interpolation 


82-94 


Index for fixed codebook, 1 " subframe 




Index for fixed codebook, 1 " subframe 


95-107 


Index for fixed codebook, 2 M subframe 




Index for fixed codebook, 2 00 subframe 


108-120 


Index for fixed codebook, 3 ,u subframe 




Index for fixed codebook, 3" subframe 


121-133 


Index for fixed codebook, 4" subframe 




Index for fixed codebook, 4" subframe 



Bit ordering of output bits from source encoder (5.8 kbit/s). 



Bits 


Description 


1-6 


Index of 1 M LSF stage 


7-12 


Index of2 nB LSF stage 


13-18 


Index of 3" LSF stage 


19-24 


Index of4 w LSF stage 


25-31 


Index of fixed and adaptive codebook gains, 1" subframe 


32-38 


Index of fixed and adaptive codebook gains, 2°° subframe 


39-45 


Index of fixed and adaptive codebook gams, 3'° subframe 


46-52 


Index of fixed and adaptive codebook gains, 4*° subframe 


53-60 


Index of pitch 


61-74 


Index for fixed codebook, 1" subframe 


75-88 


Index for fixed codebook, 2°° subframe 


89-102 


Index for fixed codebook, 3 ra subframe 


93-116 


Index for fixed codebook, 4 U subframe 



Bit ordering of output bits from source encoder (4.55 kbit/s). 



Bits 


Description 


1-6 


Index of 1" LSF stage 


7-12 


Index of2 M LSF stage 


13-18 


Index of 3" LSF stage 


19 


Index of predictor 


20-25 


Index of fixed and adaptive codebook gains, 1 N subframe 


26-31 


Index of fixed and adaptive codebook gains, 2 00 subframe 


32-37 


Index of fixed and adaptive codebook gains, 3 ,u subframe 


38-43 


Index of fixed and adaptive codebook gains, 4" subframe 


44-51 


Index of pitch 


52-61 


Index for fixed codebook, 1 H subframe 


62-71 


Index for fixed codebook, 2™ subframe 


72-81 


index for fixed codebook, 3" subframe 


82-91 


index for fixed codebook, 4 W subframe 
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APPENDIX C 



Bit ordering (channel coding) 



Ordering of bits according to subjective importance (1 1 kbit/s FRTCH). 



Bits, see table XXX 


i^vak n p uun 


1 


lai I ~v 


2 


Icfl.l 




151 1 -ii 


4 


J5I i -J 


5 


lotfl A 
ISIl-H 




IS! I -3 


7 


ISIZ-v 


g 


leO 1 


9 




10 




1 1 




12 


lsf2-5 


65 


pilvlil "V 


66 


pilCIU-l 


67 


nitrh 1 


68 


ni***h 1 -1 


69 


nit«*h 1 —A 

piicm*** 


70 


piicn 1 -J 


74 


pilLflj-V 


?5 


piUflU"! 


76 


pitch3*2 


77 


nitrh 1.1 


78 


nit*»K1_A 

piiciM*** 


79 


pitch3-5 


29 


gpl-0 


30 


gpl-I 


38 


gp2-0 


39 


BP*"* 


47 


gpj-o 


"48 ~~ 


gPJ-I 


"56* 


gp4-U 


57"" 


gp4-l 


33 


gci-u 


34 " 


gci-i 


35 


gci-^ 


42 




43 


gcz-i 


44 


gcx-x 


51 


gc3-0 


52 " 


gCJ-i 


& " 








61 


gc4-l 


62 


gc4-2 


71 


pitchI-6 


72 


pitch! -7 


73 


pitch 1-8 


80 ' 


pitch3-6 


81 


pitch3-7 


82 


pitch3-8 


83 


pttch2-0 


84 


pitch2-l 


85 


pttch2-2 


86 


pitch2-3 


87 


pitch2-4 


88 


pitch2-5 


89 


pitch4-0 


90 


pitch4-1 


91 


pitch4-2 


92 


pitch4-3 
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93 


pitch4-4 


94 


pitch4-5 


13 


lsf3-0 


14 


lsD-1 


15 


lsO-2 


16 


lsl3-3 


17 


lsD-4 


18 


lsO-5 


19 


lsf4-0 


20 


lsf4-l 


21 


lsf4-2 


22 


Isf4-3 


23 


lsf4-4 


24 


lsf4-5 


25 


lsf5-0 


26 


lsf5-l 


27 


lsf5-2 


28 


lsf5-3 


31 


gpl-2 


32 


gpl-3 


40 


gp2-2 


41 


gp2-3 


49 


gp3-2 


50 


gp3-3 


58 


gp4-2 


59 


gp4-3 


36 


gcl-3 


45 


gc2-3 


54 


gc3-3 


63 


gc4-3 


97 


excl-0 


98 


excl-1 


99 


excl-2 


100 


cxcl-3 


101 


cxcl-4 


102 


cxcl-5 


103 


cxcl-6 


104 


cxcl-7 


105 


cxci-8 


106 


cxcl-9 


107 


excl-10 


108 


cxcl-U 


109 


excl-12 


110 


cxcl-13 


111 


excl-14 


112 


excl-15 


113 


cxcl-16 


114 


cxcl-17 


115 


excl-1 8 


U6 


excl-1 9 


117 


excl-20 


118 


excl-21 


119 


excl-22 


120 


excl-23 


121 


cxcl-24 


122 


cxcl-25 


123 


excl-26 


124 


excl-27 


125 


excl-28 


128 


cxc2-0 


129 


cxc2-l 


130 


cxc2-2 


131 


cxc2-3 


132 


cxc2-4 


133 


cxc2-5 


134 


cxc2-6 


135 


CXC2-7 


136 


cxc2-8 


137 


CXC2-9 
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138 


exc2-10 


139 


cxc2-l 1 


140 


cxc2- 1 2 


141 


exc2-13 


142 


exc2-14 


143 


exc2-15 


144 


exc2-l6 


145 


exc2-17 


146 


exc2-18 


147 


exc2-l9 


148 


exc2-20 


149 


exc2-21 


150 


exc2-22 


151 


cxc2-23 


152 


exc2-24 


153 


cxc2-25 


154 


exc2-26 


155 


exc2-27 


156 


exc2-28 


159 


exc3-0 


160 


cxc3-I 


161 


exc3-2 


162 


exc3-3 


163 


cxc3^ 


164 


exc3-5 


165" 


cxc3-6 


166 


cxc3-7 


167 


cxc3-8 


168 


exc3-9 


169 


exc3-10 


170 


cxc3-l 1 


171 


cxc3-12 


172 


exc3-13 


173 


exc3-14 


174 


exc3-15 


175 


exc3-16 


176 


exc3-17 


177 


exc3-l 8 


178 


exc3-19 


179 


exc3-20 


180 


exc3-2 1 


181 


exc3-22 


182 


cxc3-23 


183 


exc3-24 


184 


exc3-25 


185 


exc3-26 


186 


cxc3-27 


187 


exc3-28 


190 


exc4-0 


191 


exc4-l 


192 


cxc4-2 




exc4-3 


194 


cxc4^ 


195 " " 


exc4-5 


196 


cxc4*6 


197 


exc4-7 


198 


exc4-8 


199 


cxc4-9 


200 


cxc4-10 


201 


exc4.ll 


102 


cxc4-12 


203 


cxc4-13 


264 


exc4-l4 


205 


cxc4-I5 


206 


exc4-16 


267 - 


CXC4-17 


208 


cxc4-18 


209 


cxc4-I9 


210 


cxc4-20 
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211 


cxc4-21 


212 


exc4-22 


213 


cxc4-23 


214 


exc4-24 


215 


cxc4-25 


216 


cxc4-26 


217 


cxc4-27 


218 


exc4-28 


37 


gcM 


46 


gc2~4 


55 


gc3-4 


64 


gc4-4 


126 


excI-29 


127 


cxcl-30 


157 


cxc2-29 


158 


cxc2-30 


188 


cxc3-29 


189 


cxc3-30 


219 


cxc4-29 


220 


cxc4-30 


95 


interp-0 


96 


interp-1 
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Ordering of bits according to subjective importance (8.0 kbit/s FRTCH). 



Bits, see table XXX 


Description 


1 


lsfl-0 


2 


lsn-1 


3 


lsn-2 


4 


!sfl-3 


5 


lsfl-4 


6 


lsfl-5 


7 


lsf2-0 


8 


lsf2-l 


9 


lsf2-2 


10 


lsf2-3 


11 


lsf2-4 


12 


Isf2-5 


25 


gain 1-0 


26 


gain 1-1 


27 


gam 1-2 


28 


gain 1-3 


29 


gam 1-4 


32 


gain2-0 


33 


gain2-l 


34 


gain2-2 


35 


gain2-3 


36 


gam2-4 


39 


gain3-0 


40 


gain3-l 


41 


gain3-2 


42 


gain3-3 


43 


gain3-4 


46 


gain4-0 


47 


gain4-l 


48 


gain4-2 


49 


gain4-3 


50 


gam4-4 


53 


pitch 1-0 


54 


pUchl-1 


55 


pitch 1-2 


56 


pitch 1-3 


57 


pttchl-4 


58 


pitch) -5 


61 


pitch3-0 


62 


pitch3-l 


63 


pitch3-2 


64 


pitch3-3 


65 


pitch3~4 


66 


pitch3-5 


69 


pitch2-0 


70 


pitch2-l 


71 


pitch2-2 


74 


pitch4-0 


75 


pitch4-l 


76 


pttch4-2 


13 


IsO-O 


14 


lsb-1 


15 


ls*3-2 


16 


lsO-3 


17 


lsO-4 


IS 


UO-5 


30 


gam 1-5 


J7 


gam2-5 


44 


gain3-5 


51 


gam4-5 


59 


pitch 1-6 


67 


pitch3-6 


72 ■ 


pitch2-3 


77 


pttch4-3 


79 


interp-0 


80 


mterp-1 1 
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31 


gain 1 -6 


38 


gain2-6 


45 


gain3-6 


52 


gain4-6 


19 


lsf4-0 


20 




21 


lsf4-2 


22 


!sf4-3 


23 


lsf4-4 


24 


lsf4-5 


60 


pitch 1-7 


68 


pitch3*7 


73 


pitch2-4 


78 


pitch4-4 


81 


cxcl-0 


82 


excl-1 


83 


excl-2 


84 


exc 1 -3 


85 


excl-4 


86 


exc 1-5 


87 


exc 1-6 


88 


exc 1-7 


89 


exc 1-8 


90 


exc 1-9 


91 


exc 1-10 


92 


excl-1 1 


93 


excl-1 2 


94 


excl-13 


95 


exc 1-14 


96 


excl-15 


97 


excl-1 6 


98 


excl-1 7 


99 


excl-1 8 


100 


excl-1 9 


101 


exc2-0 


102 


exc2-l 


103 


cxc2-2 


104 


exc2-3 


105 


exc2-4 


106 


exc2-5 


107 


exc2-6 


108 


exc2-7 


109 


exc2-8 


110 


exc2-9 


111 


exc2-10 


112 


exc2-1 1 


113 


exc2-12 


114 


exc2-13 


115 


cxc2-14 


116 


exc2-15 


117 


exc2-16 


118 


exc2-17 


119 


exc2-18 


120 


exc2-19 


121 


exc3-0 


122 


cxc3-l 


123 


exc3-2 


124 


exc3-3 


125 


exc3-4 


126 


exc3-5 


127 


exc3-6 


128 


exc3-7 


129 


exc3-8 


130 


cxc3-9 


131 


exc3-IO 


132 


exc3.Il 


133 


exc3-12 


134 


exc3-l3 


135 


exc3-l4 
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1 1f\ 

1 JO 


excj-i 3 


1 17 " 

i j t 


exc3- 1 6 


1 IS 


avj*1 1 *7 


139 


^v/*l 1 ft 


140 


CaCJ* t zf 


141 




142 


cxc4-l 


143 


cxc4*2 


144 


cxc4*3 


145 




146 


cxc4o 


147 


cxc4"6 


148 




149 


exc4-8 


150 




151 


CXC4-10 


152 


CXC4-1 1 


153 


cxc4-12 


154 


exc4-13 


155 


cxc4-14 


156 


cxc4-15 


157 


exc4-!6 


158 


cxc4-17 


159 


exc4-18 


160 


cxc4-19 
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Ordering of bits according to subjective importance (6.65 kbit/s FRTCH). 



Bits, see table XXX 


Description 


54 


pitch-0 


55 


pitch - 1 


5<i 


pitch -2 


57 


pitch-3 


58 


pitch -4 


59 


pitch-5 


1 


lsfl-0 


2 


Isfl-1 


3 


lsfl-2 


4 


lsfl-3 


5 


lsfl-4 


6 


lsfl-5 


25 


gain 1*0 


26 


gain 1-1 


27 


gain 1 -2 


28 


gain 1-3 


32 


gain2-0 


33 


gain2-l 


34 


gain2-2 


35 


gain2-3 


39 


gain3-0 


40 


gain3-l 


41 


gain3-2 


42 


gain3-3 


46 


gain4-0 


47 


gain4-l 


48 


gain4-2 


49 


gain4-3 


29 


gain 1-4 


36 




43 


gain3-4 


50 


gain4-4 


53 


mode-0 


98 


exc3-0 niteh-OfSeermfi Qiihframpl 


99 


exc_l-l nitr:K-l ( ^rrrmH ciiKtfram^t 


7 


ls£2-0 


8 


lsf2-l 


9 


lsf2-2 


10 


lsf2-3 


1 1 


lsQ-4 


12 


ls!2-5 


30 


gain 1 -5 


37 


gain2-5 


44 


gain3-5 


51 


gatn4-5 


62 


excl-0 pitch-OfThird subrrame) 


63 


exc 1 - 1 pitch- 1 (Third subrrame) 


64 


exc 1-2 pitch-2(Third subrrame) 


65 


exc 1-3 pitch-3(Third subframe) 


66 


exc 1-4 pitch-4(Third subrrame) 


80 


exc2-0 pitch-5(Third subframe) 


100 


exc3-2 pi tch-2( Second subrrame) 


116 


exc4-0 pitch-0( Fourth subrrame) 


117 


exc4-l pitch-1 (Fourth subframe) 


118 


exc4-2 Ditch-2f Fourth subframe 1 


13 


tsf3-0 


14 


IsO-l 


15 


1SO-2 


16 


lsD-3 


17 


lsO-4 


18 


Isf3-5 


19 


Isf4-0 


20 


lsf4-l 


21 


lsf4-2 


22 


lsf4-3 


67 


exc 1-5 exc I dtp) 


68 


exc 1-6 exc 1 (Itp) 
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69 


excl-7 excl(ltp) 


70 


excl-8 excl(ltp) 


71 


excl-9 cxcl(ltp) 


72 


excl-10 


81 


exc2-l exc2(ltp) 


82 


exc2-2 exc2(ltp) 


83 


cxc2-3 cxc2(ltp) 


84 


cxc2-4 exc2(ltp) 


85 


exc2-5 exc2(ltp) 


86 


exc2-6 exc2(ltp) 


87 


cxc2-7 


88 


exc2-8 


89 


cxc2-9 


90 


exc2-10 


101 


exc3-3 exc3(ltp) 


102 


exc3-4 exc3(ltp) 


103 


exc3-5 exc3<ltp) 


104 


exc3-6 exc3(ltp) 


105 


cxc3-7 cxc3(ltp) 


106 


cxc3-8 


107 


cxc3-9 


108 


exc3-10 


119 


exc4-3 exc4(ltp) 


120 


exc4-4 cxc4(ltp) 


121 


exc4-5 exc4(ltp) 


122 


exc4-6 exc4(ltp) 


123 


exc4-7 exc4(ltp) 


124 


exc4-8 


125 


exc4~9 


126 


exc4-10 


73 


excl-11 


91 


exc2-l 1 


109 


cxc3-ll 


127 


cxc4-l 1 


74 


cxcl-12 


92 


cxc2-12 


110 


exc3-12 


128 


cxc4-12 


66 


pitcb-6 


61 


pitch-7 


23 


lsf4-4 


24 


Uf4-5 


75 


excl-13 


93 


exc2-13 


in 


cxc3-13 


129 


cxc4-13 


31 


gain 1-6 


38 


gam2-6 


45 


jgani3-6 


52 


gain4-6 


16 


excl-14 


77 


excl-15 


94 


exc2-14 


95 


cxc2-15 


112 


cxc3-14 


113 


exc3-l5 


130 


CXC4-I4 


131 


cxc4-15 


78 


cxcl-16 


96 


cxc2-16 


114 


exc3-16 


\ii 


cxc4-16 


79 


excl-17 


97 


cxc2-17 


115 


cxc3-17 


133 


CXC4-I7 
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Ordering of bits according to subjective importance (5.8 kbit/s FRTCH). 



Bits, see table XXX 


Description 


53 


pitch -0 


54 


pitch- 1 


55 


pttch-2 


56 


pitch-3 


57 


pitch-4 


58 


pitch-5 


1 


lsH-0 


2 


isn-i 


3 


lsfl-2 


4 


isn-3 


5 


lsfl-4 


6 


lsfl-5 


7 


ls£2-0 


8 


lsf2-l 


9 


lsG-2 


10 


lsG-3 


11 


lsf2-4 


12 


lsC-5 


25 


gainl-0 


26 


gain 1-1 


27 


gain 1-2 


28 


gain 1-3 


29 


gain 1-4 


32 


gain2-0 


33 


gain2-l 


34 


gam2-2 


35 


gain2-3 


36 


gain2-4 


39 


gain3-0 


40 


gatn3-l 


41 


gain3-2 


42 


gain3-3 


43 


gain3-4 


46 


gain4-0 


47 


gain4-l 


48 


gatn4-2 


49 


gain4-3 


50 


gain4-4 


30 


gain 1-5 


37 


gain2-5 


44 


gam3-5 


51 


gain4-5 


13 


lsO-0 


14 


lsl3-l 


15 


sO-2 


16 


lsD-3 


17 


sO-4 


18 


sO-5 


59 


pitch-6 


60 


pitch-7 


19 


lSf4-0 


20 


Isf4*l 


21 


SI4-2 


22 


IS14-3 


23 


lSf4-4 


24 


lsf4-5 


31 


gainl-6 


38 


gain2-6 


45 


gam3-6 


52 


gam4-6 


61 


excl-0 


75 


exc2-0 


89 


exc3-0 


103 


exc4-0 


62 


excl-l 


63 


excl-2 
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64 


cxc 1 "3 


65 


cxc 1 ^4 


66 


excl-5 


67 


excl-6 


68 


cxc 1 -7 


69 


cxc 1-8 


70 


cxc 1 -9 


71 


excl-10 


72 


cxc 1*1 1 


73 


cxc I -12 


74 


excl-13 


76 


exc2-l 


77 


exc2-2 


78 


cxc2-3 


79 


exc2-4 


80 


exc2-5 


81 


exc2-6 


82 


cxc2-7 


83 1 


cxc2-8 


84 


exc2*9 


85 


exc2-l0 


86 


exc2-l 1 


87 


exc2-12 


88 


cxc2-13 


90 


exc3-l 


91 


exc3-2 


92 ' 


cxc3-3 


93 


cxc3*4 


94 


cxc3-5 


95 


cxc3-6 


96 


exc3-7 


97 


cxc3-& 


98 


exc3-9 


99 


exc3-10 


100 


exc3-l 1 


101 


cxc3-12 


102 


exc3-13 


104 


CXC4-1 


105 


exc4-2 


106 


exc4-3 


107 


exc4-4 


108 


exc4-5 


109 


exc4-6 


110 


cxc4-7 


111 


cxc4-8 


112 - 


cxc4-9 


113 


cxc4-10 


114 


CXC4-1 1 


115 


cxc4-12 


116 


cxc4-13 
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Ordering of bits according to subjective importance (8.0 kbil/s HRTCH) 



Bits, see table XXX 


Description 


I 


lsfl-0 


2 


lsfl-1 


3 


lsfl-2 


4 


lsfl-3 


5 


lsfl-4 


6 


lsfl-5 


25 


gain 1 -0 


26 


gain 1 - 1 


27 


gain 1-2 


28 


fraif) l -3 


32 


gain2-0 


33 


oaip2- 1 


34 


g3in2~2 


35 


gain2-3 


39 


gatfi3**0 


40 




41 


gain3-2 


42 


gain3-3 


46 


gain4-0 


47 


i?ain4-l 


48 


gain4-2 


49 


oain4-3 


53 


pitch 1-0 


54 


pitch I -1 


55 


pitch 1 -2 


56 


pitch I -3 


57 


pitch 1-4 


58 


pitch 1-5 


61 


flitch*) -fl 


62 


pitch3-l 


63 


pitch3-2 


64 


pitch3-3 


65 


pitch3*4 


& 


fit IMIJ J 


69 




70 


pitch2-l 


71 


pitch 2-2 


74 


pitch4-0 


75 


pi^h4~l 


76 


pitch4-2 


7 


lsf2-0 


8 


Isf2-1 


9 


lsf2-2 


10 


lsf2-3 


11 


lsO-4 


12 


lsf2-5 


29 


gain 1-4 


36 


gain2-4 


43 


gain3-4 


50 


gain4-4 


79 


interp-0 


SO 


tntcrp-l 


13 


lsf3-0 


14 


IsO-l 


15 


lsO-2 


16" 


lsO-3 


17 


lsO-4 


18 


lsG-5 


19 


lsf4-0 


20 


lsf4-l 


21 


lsf4-2 


22 


lsf4-3 


23 


lsf4-4 


24 


lsf4-5 


30 


gain 1-5 


31 


gain 1-6 


37 


gain2-5 
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38 


gain2-6 


44 


gain3-5 


45 


gain3-6 


51 


gain4-5 


52 


gairv4-6 


59 


pitch 1 -6 


67 


pitch3-6 


72 


pitch2-3 


77 


pitch4-3 


60 


pitch 1-7 


68 


pitch3-7 


73 


pitch2-4 


78 


pitch4-4 


81 


excl-0 


82 


excl-1 


83 


excl-2 


84 


excl-3 


85 


cxcl -4 


86 


excl-5 


87 


cxc 1 "6 


88 


cxcl -7 


89 


cxc 1-8 


90 


cxcl -9 


91 


cxcl -10 


92 


cxc 1-1 1 


93 


excl-1 2 


94 


excl-1 3 


95 


excl-14 


96 


excl-1 5 


97 


excl-1 6 


98 


excl-1 7 


99 


excl-1 8 


100 


excl-1 9 


101 


exc2-0 


102 


exc2-l 


103 


exc2-2 


104 


exc2-3 


105 


cxc2-^ 


106 


exc2-5 


107 


exc2-6 


108 


exc2-7 


109 


cxc2-8 


110 


exc2-9 


111 


exc2-10 


112 


exc2-l 1 


113 


cxc2-12 


114 


exc2-13 


115 


cxc2-14 


116 


exc2-15 


117 


exc2-16 


118 


exc2-l 7 


119 


cxc2-18 


120 


exc2-19 


121 


exc3^> 


122 


exc3-l 


123 


exc3-2 


124 


cxc3-3 


125 


exc3-4 


126 


exc3-5 


127 


exc3-6 


128 


CXC3-7 


129 


exc3-8 


130 


exc3-9 


131 


exc3-10 


132 


exc3-ll 


133 


exc3-12 


134 


exc3-13 


135 


exc3-14 


136 


cxc3-l5 
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137 


exc3-16 


138 


exc3-17 


139 


exc3- 18 


1 Ai\ 

140 


exc3-19 


1 A 1 

141 


exc4-0 




exc4- 1 


\A1 


exc4-2 


| A A 


CXC4-J 


i 4 * j 


exc4^ 


\AA 
I HO 


€xc4-5 


1*9 t 


CXC4-0 


\ASt 
l4o 


CXC4-/ 




exc4-8 




cxc4-9 


151 




152 


cxc4-U 


153 


exc4-12 


154 


exc4-13 


155 


exc4-14 


156 


exc4-15 


157 


cxc4-16 


158 


cxc4-17 


159 


cxc4-I8 


160 


exc4-l9 



RWQIYYMfV ^Uki/-v 
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Ordering of bits according to subjective importance (6.65 kbtt/s HRTCH). 



dii»» see laoie aaa 


Description 


j j 


mode-0 
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CLAIMS 

We claim: 

1 . A speech encoding system using an analysis by synthesis approach on a 
speech signal, the speech encoding system comprising: 

an encoder that calculates open loop gain and closed loop gain; 
a plurality of codebooks that generate excitation vectors; and 

the encoder selectively applying the open loop gain and the closed loop gain in gain 
normalization processing. 

2. The speech encoding system of claim 1 wherein the encoder smoothly 
combines the open loop gain and the closed loop gain. 

3. The speech encoding system of any of claims 1 and 2 wherein the plurality of 
codebooks comprise an adaptive codebook and a fixed codebook. 

4. The speech encoding system of claim 3 wherein the encoder generates an 
adaptive codebook gain and a fixed codebook gain, and the encoder applies the gain 
normalization factor to both the adaptive codebook gain and the fixed codebook gain. 

5. The speech encoding system of claim 4 wherein the smooth combination of 
the open loop gain and the closed loop gain by the encoder comprises utilizing linear 
predictive coding gain. 
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6. The speech encoding system of claim 2 wherein the smooth combination of 
the open loop gain and the closed loop gain by the encoder comprises utilizing linear 
predictive coding gain. 

7. The speech encoding system of claim 2 wherein the smooth combination of 
the open loop gain and the closed loop gain by the encoder comprises a linear combination of 
the open loop gain and the closed loop gain using linear predictive coding gain as a weighting 
factor. 

8. The speech encoding system of any of claims 2 and 7 wherein the encoder 
excludes background noise from the smooth combination of the open loop gain and the 
closed loop gain. 

9. The speech encoding system of claim 8 wherein the encoder applies the lesser 
of the open loop gain and the closed loop gain to the background noise. 

10. A method used by a speech encoding system that applies an analysis by 
synthesis coding approach to a speech signal, the method comprising: 

identifying open loop gain; 
identifying closed loop gain; and 

combining contributions of the open loop gain and the closed loop gain to generate a 
gain normalization factor. 
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1 1 . The method of claim 1 0 further comprising: 

utilizing linear predictive coding gain to identify appropriate contributions of 
the open loop gain and the closed loop gain. 

12. The method of claim 10 further comprising: 

utilizing linear predictive coding gain as a weighting factor to identify the 
contributions of the open loop gain and the closed loop gain. 

13. The method of any of claims 10, 1 1 and 12 further comprising: 

selecting without combination either the open loop gain or the closed loop gain in 
generating the gain noimalization factor when the speech signal comprises background noise. 

14. The method of any of claims 10, 1 1 and 12 further comprising: 
choosing the lesser of the open loop gain and the closed loop gain without 

combination in generating the gain normalization factor when the speech signal comprises 
background noise. 
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1 1 . The method of claim 10 further comprising: 
utilizing linear predictive coding gain to identify appropriate contributions of the open 
loop gain and the closed loop gain. 



12. The method of claim 10 further comprising: 

utilizing linear predictive coding gain as a weighting factor to identify the 
contributions of the open loop gain and the closed loop gain. 



13. The method of any of claims 10, 1 1 and 12 further comprising: 
selecting without combination either the open loop gain or the closed loop gain in 
generating the gain normalization factor when the speech signal comprises background noise. 



14. The method of any of claims 10, 1 1 and 12 further comprising: 
choosing the lesser of the open loop gain and the closed loop gain without 

combination in generating the gain normalization factor when the speech signal comprises 

background noise. 
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